Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,303,840 Δ1 bp 3.0% intergenic (+27/+130) urtE → / ← RCAP_rcc01231 urea ABC transporter ATP‑binding protein UrtE/UbiH/UbiF/VisC/COQ6 family ubiquinone biosynthesis hydroxylase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,303,8400C→. 3.0% 4.7 ‑4.0 ‑0.1132intergenic (+27/+130)urtE/RCAP_rcc01231urea ABC transporter ATP‑binding protein UrtE/UbiH/UbiF/VisC/COQ6 family ubiquinone biosynthesis hydroxylase
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (46/82):  total (50/82)
Fisher's exact test strand distribution p-value = 1.91e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.08e-01

Only 100 of 149 total aligned reads displayed.

GGCGCGGTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCT.C..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGCC  >  NC_014034/1303742‑1303921
                                                                                                       |                                                                                  
ggcgcggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGc                                                                                    <  1:2996112/101‑1 (MQ=255)
   gcggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTc                                                                                 <  2:3003291/101‑1 (MQ=255)
   gcggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTc                                                                                 <  1:1495276/101‑1 (MQ=255)
    cggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGGGCTGCCCCC.C..CGGCTTCCCCTcc                                                                                >  1:4393033/1‑100 (MQ=255)
    cggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCa                                                                                <  2:1952793/101‑1 (MQ=255)
    cggTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCATGGTCAGCGTTTAGGCAAGGGGTCGGTG..CCC.CCTCGGCTTCGCCTcc                                                                                >  2:3427209/1‑100 (MQ=255)
       tGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCAccc                                                                             <  2:3237410/101‑1 (MQ=255)
        ggCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCAcccc                                                                            <  1:656928/101‑1 (MQ=255)
        ggCGATGGAGGGCACAAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCAcccc                                                                            <  1:1096884/101‑1 (MQ=255)
         gCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCCccccc                                                                           >  1:1302840/1‑101 (MQ=255)
         gCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCAccccc                                                                           <  2:3536889/101‑1 (MQ=39)
         gCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGATGCCCCC.C..CGGCTTCCCCTCGccccc                                                                           >  2:4407804/1‑101 (MQ=255)
          cGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCCCCTCACCCCCg                                                                          >  2:2532750/1‑101 (MQ=255)
           gATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCgg                                                                         >  1:3152089/1‑101 (MQ=39)
            aTGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCggg                                                                        <  2:1341999/101‑1 (MQ=39)
            aTGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCggg                                                                        >  1:3513776/1‑101 (MQ=39)
            aTGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCggg                                                                        <  1:3127328/101‑1 (MQ=39)
             tgtaGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGa                                                                       <  2:471683/98‑1 (MQ=38)
             tGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCgggg                                                                       >  1:3063966/1‑100 (MQ=39)
              ggAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGat                                                                      <  1:248800/101‑1 (MQ=38)
               gAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGGTG..CCC.CCTCGGCTTGGCCTCACCCCCGGGata                                                                     >  2:2988052/1‑101 (MQ=38)
               gAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGata                                                                     <  1:2018953/101‑1 (MQ=38)
               gAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGata                                                                     <  1:1209026/101‑1 (MQ=38)
                 gggCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATAtt                                                                   <  1:618758/101‑1 (MQ=38)
                  ggCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATAttt                                                                  <  2:1246470/101‑1 (MQ=38)
                   gCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTg                                                                 <  2:1886666/101‑1 (MQ=38)
                    cACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTgg                                                                >  1:350301/1‑101 (MQ=38)
                    cACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTgg                                                                <  1:2164282/101‑1 (MQ=38)
                     aCGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGCGCTTTggg                                                               >  2:4438322/1‑101 (MQ=38)
                       gAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGc                                                             <  1:1036236/101‑1 (MQ=38)
                          ggACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGGTATTTGGGGGAAg                                                          >  2:2348740/1‑101 (MQ=37)
                          ggACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAg                                                          <  2:2915276/101‑1 (MQ=37)
                           gACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGa                                                         >  1:2789282/1‑101 (MQ=37)
                             cAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.a                                                      <  1:2940153/101‑1 (MQ=37)
                             cAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.a                                                      <  2:434002/101‑1 (MQ=37)
                               gCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.aaa                                                    <  2:3607043/101‑1 (MQ=37)
                               gCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.aaa                                                    <  1:3667202/101‑1 (MQ=37)
                                 cTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGCGAAATTTGGGGCAAGAT.aaaac                                                  >  1:2367926/1‑99 (MQ=25)
                                  tGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAccc                                                 >  1:3005976/1‑101 (MQ=37)
                                  tGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAccc                                                 >  2:3876525/1‑101 (MQ=37)
                                   gTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAAcccc                                                >  2:447535/1‑101 (MQ=37)
                                    tCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCa                                               <  1:4521256/101‑1 (MQ=37)
                                    tCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGCGATATTTGGGGCAAGAT.AAACCCCa                                               >  1:519772/1‑101 (MQ=37)
                                     cccGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCAt                                              >  1:2116037/1‑101 (MQ=37)
                                        gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                           <  1:3385630/101‑1 (MQ=37)
                                        gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                           >  1:324403/1‑101 (MQ=37)
                                        gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                           <  1:1345128/101‑1 (MQ=37)
                                        gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                           <  2:3374429/101‑1 (MQ=37)
                                         cgACACCCTGCGCGCAGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                          <  2:2820085/101‑1 (MQ=25)
                                         cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                          <  2:40047/101‑1 (MQ=37)
                                         cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                          >  2:3809545/1‑101 (MQ=37)
                                         cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                          <  2:3753589/101‑1 (MQ=37)
                                           aaaCCCTGCGCGGCGTGGTCAGCGTTTAGGCAAGGGGGCGCTGCCCCC.C..CGGCTTCGCCTCACCCCCGGGCGATTTGGGGCAAGAT.AAACCCCATCCTTTc                                        >  1:2520620/1‑101 (MQ=37)
                                           aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                        >  2:2915590/1‑101 (MQ=37)
                                           aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                        >  2:672659/1‑101 (MQ=37)
                                           aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCCCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                        >  2:1331851/1‑101 (MQ=37)
                                            aaCCCTGCGCGCCGTTGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCGCCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCa                                       <  1:4612286/101‑1 (MQ=25)
                                             aCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCAc                                      <  2:350269/101‑1 (MQ=37)
                                              cccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCAcc                                     >  1:3069401/1‑101 (MQ=37)
                                              cccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCAcc                                     >  2:762776/1‑101 (MQ=37)
                                               ccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCt                                    >  2:1356017/1‑101 (MQ=37)
                                                tgcagcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCtt                                   >  2:2187067/5‑101 (MQ=25)
                                                cTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCtt                                   <  1:3809830/101‑1 (MQ=37)
                                                 tGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTg                                  >  2:2185049/1‑101 (MQ=37)
                                                  gcgcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGc                                 <  2:4118282/101‑1 (MQ=37)
                                                   cgcggCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGcc                                <  2:2964949/101‑1 (MQ=25)
                                                    gcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTGCCCCC.C..CGGCTTCGCCTCACCCCCGGGGTATTTTTGGCAAGAT.AAACCCCATCCTTTCTCCTTTCCg                               >  1:1575536/1‑101 (MQ=25)
                                                    gcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCg                               >  1:3830234/1‑101 (MQ=37)
                                                     cgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGc                              <  1:489274/101‑1 (MQ=37)
                                                      gcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCa                             <  2:3567447/101‑1 (MQ=37)
                                                       ccGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaa                            <  1:2874406/101‑1 (MQ=37)
                                                       ccGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaa                            >  1:1585085/1‑101 (MQ=37)
                                                        cGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCCCCCCCCGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaaa                           >  1:777301/1‑101 (MQ=37)
                                                         gTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGGTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAat                          >  2:2545248/1‑101 (MQ=37)
                                                         gTGGTCAGCGTTTAGGCAAGGGAGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAat                          <  2:4474152/101‑1 (MQ=37)
                                                          tGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                         <  1:2986212/101‑1 (MQ=37)
                                                          tGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                         <  2:4252110/101‑1 (MQ=37)
                                                          gggTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCGCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                         <  2:4362535/100‑1 (MQ=25)
                                                           ggTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAatat                        <  1:4116336/101‑1 (MQ=37)
                                                           ggTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAatat                        <  1:1129550/101‑1 (MQ=37)
                                                            gTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATc                       <  1:2185199/101‑1 (MQ=37)
                                                             tCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATcc                      >  1:4187963/1‑101 (MQ=37)
                                                               aGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCg                    >  1:3606527/1‑101 (MQ=25)
                                                               aGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCg                    <  1:3876810/101‑1 (MQ=37)
                                                                gCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                   <  1:442255/101‑1 (MQ=37)
                                                                gCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                   <  2:2272683/101‑1 (MQ=37)
                                                                  gTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgggg                 <  1:447572/101‑1 (MQ=37)
                                                                  ctgacggacggggggCTCTGC.CCC.C..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCaggggg               >  1:3167541/11‑101 (MQ=2)
                                                                  ctgacggacggggggCTCTGC.CCC.C..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCaggggg               >  2:2585515/11‑101 (MQ=2)
                                                                     tAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCCAATATCCCggggggc              >  2:567525/1‑100 (MQ=25)
                                                                      aGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgggggggg             >  1:541675/1‑99 (MQ=25)
                                                                      aGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTg             <  1:457214/101‑1 (MQ=37)
                                                                      aGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTg             <  1:2882259/101‑1 (MQ=37)
                                                                          aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt         <  1:762851/101‑1 (MQ=37)
                                                                          aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt         >  2:3190409/1‑101 (MQ=25)
                                                                          aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCGATCACCTTGCCGCAAATATCCCGGGGGGTGAAtt         >  1:4005449/1‑101 (MQ=25)
                                                                           aGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTg        <  1:3012741/101‑1 (MQ=37)
                                                                              gggCTCTG..CCCCC..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGG.GAATTGgccaaa  >  1:3204317/1‑98 (MQ=2)
                                                                              gggCTCTG..CCCCC..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGgccaaa  >  1:2983703/1‑98 (MQ=2)
                                                                               ggCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGgccg    <  2:3069157/101‑1 (MQ=37)
                                                                                                       |                                                                                  
GGCGCGGTGGCGATGGAGGGCACGAAGGACAGCCTGTCCCGCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCT.C..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGCC  >  NC_014034/1303742‑1303921

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG