Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,177 +C 91.2% intergenic (‑105/+120) RCAP_rcc01564 ← / ← mmsA isobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,696,1771.→C91.2% 64.7 ‑2.5 ‑1.0262intergenic (‑105/+120)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (106/133):  ref base (15/8):  total (121/141)
Fisher's exact test strand distribution p-value = 7.83e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.11e-01

Only 100 of 290 total aligned reads displayed.

GCTCCTCCCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCG.G.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGG  >  NC_014034/1696075‑1696275
                                                                                                        |                                                                                                  
gcTCCTCCCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCg.g                                                                                                     >  1:5454871/1‑101 (MQ=255)
 ctcctcCCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.c                                                                                                   <  1:5193728/101‑1 (MQ=255)
  tcctccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.cc                                                                                                  >  1:1624267/1‑100 (MQ=255)
   cctccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.cca                                                                                                 <  1:2767065/101‑3 (MQ=255)
   cctccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.cca                                                                                                 <  2:1634982/101‑3 (MQ=255)
    ctccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.ccaa                                                                                                >  2:4809386/1‑98 (MQ=255)
     tccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.ccaat                                                                                               >  1:4489797/1‑97 (MQ=255)
        cataGTGCCTGAATTGAACGTCTGTTGAATTTCGAATCCGTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTca                                                                                             <  1:83130/98‑2 (MQ=255)
         aaaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTcac                                                                                            >  2:2106293/1‑99 (MQ=255)
           aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAccc                                                                                         >  2:2900889/1‑101 (MQ=255)
           aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAccc                                                                                         >  1:2403507/1‑101 (MQ=255)
           aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAccc                                                                                         >  1:1724361/1‑101 (MQ=255)
           aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAccc                                                                                         >  1:5414966/1‑101 (MQ=255)
           aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCCTTCG.G.cccattcaccc                                                                                         >  2:1110872/1‑91 (MQ=255)
            gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAcccc                                                                                        >  1:5060482/1‑100 (MQ=255)
            gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAcccc                                                                                        <  1:935325/101‑2 (MQ=37)
             tGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACTAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAccccg                                                                                       <  1:5768081/101‑3 (MQ=255)
             tGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAccccg                                                                                       <  1:1151131/101‑3 (MQ=255)
              gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCAcccccg                                                                                      >  1:3108476/1‑98 (MQ=37)
              gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCACCCcgg                                                                                      >  2:4121638/1‑101 (MQ=37)
              gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.G.CCAATTCACCCcgg                                                                                      <  2:799158/101‑1 (MQ=255)
                cTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGa                                                                                     <  1:1496748/101‑1 (MQ=255)
                 tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGAC.CAATTCACCCCGGGat                                                                                    <  2:5244090/101‑2 (MQ=255)
                 tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCggcat                                                                                    <  1:3298167/101‑4 (MQ=255)
                  gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGatg                                                                                   >  2:3510539/1‑99 (MQ=255)
                  gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGatg                                                                                   <  1:2706380/101‑3 (MQ=25)
                  gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGatg                                                                                   >  2:4697478/1‑99 (MQ=255)
                   aaTTTATCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAG.AGGGTCGC.TCCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGatgt                                                                                  <  2:4616024/101‑4 (MQ=255)
                   aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGatgt                                                                                  >  1:649354/1‑98 (MQ=255)
                   aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGatgt                                                                                  <  1:1552590/101‑4 (MQ=255)
                    aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGatgtt                                                                                 >  1:425845/1‑97 (MQ=25)
                     ttGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGatgttt                                                                                <  2:3651499/101‑6 (MQ=255)
                      tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTg                                                                               >  1:1759181/1‑101 (MQ=255)
                       gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCgggctgttagg                                                                              >  2:4877616/1‑93 (MQ=25)
                        aaCGTCTGGTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGGCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGCTGTTTggg                                                                             >  2:5763773/1‑101 (MQ=255)
                          cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGcc                                                                           >  2:1045481/1‑101 (MQ=25)
                           gtgtGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCa                                                                          <  1:76336/98‑1 (MQ=19)
                           gTCTTTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCa                                                                          <  2:4231027/101‑1 (MQ=19)
                           gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCa                                                                          >  2:1848175/1‑101 (MQ=255)
                           gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCa                                                                          >  2:9278/1‑101 (MQ=255)
                            tCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaa                                                                         >  2:1934612/1‑101 (MQ=255)
                            tCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaa                                                                         >  2:4469657/1‑101 (MQ=255)
                             cTGTTTAATTTCGAATCCTTGCACGACTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaag                                                                        <  1:3478717/101‑2 (MQ=255)
                             cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaag                                                                        >  2:2626880/1‑100 (MQ=21)
                              tGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaagg                                                                       <  2:3344486/101‑3 (MQ=19)
                              tGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaagg                                                                       <  1:251494/101‑3 (MQ=19)
                                tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGCGG..TGCCCCCCCCGGCCCTTCG.GCCCATTCACCCCCGGGTGTTTGGGCCAAGggg                                                                     >  2:2063718/1‑101 (MQ=255)
                                tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                     <  1:3778367/101‑5 (MQ=255)
                                 ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                    <  1:1214201/101‑6 (MQ=16)
                                  tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtgaa                                                                   <  1:5644703/101‑7 (MQ=255)
                                   aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                  <  2:5387312/101‑8 (MQ=255)
                                   aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                  >  2:3585925/1‑94 (MQ=16)
                                     tttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAaca                                                              >  1:5313002/1‑100 (MQ=255)
                                      ttCGAATCCGTGCGCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacac                                                             <  1:1574542/101‑3 (MQ=255)
                                       tCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacc                                                            <  1:974318/101‑4 (MQ=12)
                                         gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                          >  2:2250233/1‑96 (MQ=12)
                                         gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                          <  2:725613/101‑6 (MQ=255)
                                         gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                          <  2:2394705/101‑6 (MQ=255)
                                           aTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgc                                                        <  2:545941/101‑8 (MQ=255)
                                            tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                       >  2:3157585/1‑93 (MQ=12)
                                             ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                      >  2:2540259/1‑92 (MQ=255)
                                              cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                     >  1:5080158/1‑91 (MQ=255)
                                                 gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                <  2:1598833/101‑1 (MQ=255)
                                                 gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                <  1:4745554/101‑1 (MQ=255)
                                                  cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                               <  1:2889633/101‑1 (MQ=255)
                                                    cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                             >  1:101160/1‑101 (MQ=255)
                                                          ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                       >  2:4568820/1‑101 (MQ=255)
                                                          ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                       >  1:4806610/1‑101 (MQ=255)
                                                           tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                      <  1:4084682/101‑1 (MQ=255)
                                                           tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                      <  1:2915783/101‑1 (MQ=255)
                                                           tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                      >  1:4960620/1‑101 (MQ=9)
                                                             cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                    <  2:5337219/101‑1 (MQ=255)
                                                             cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                    <  1:2215698/101‑1 (MQ=255)
                                                               aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                  <  1:4121887/101‑1 (MQ=9)
                                                                 gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                <  1:4582439/101‑1 (MQ=255)
                                                                 gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                >  1:390964/1‑101 (MQ=9)
                                                                  aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                               <  2:1464180/101‑1 (MQ=9)
                                                                  aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                               <  1:3157774/101‑1 (MQ=9)
                                                                    gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                             <  2:1376288/101‑1 (MQ=255)
                                                                    gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                             <  1:3586142/101‑1 (MQ=255)
                                                                    gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                             <  2:4290626/101‑1 (MQ=255)
                                                                         agaggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGc                          <  2:286209/91‑1 (MQ=255)
                                                                         agaggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGc                          >  2:389819/11‑101 (MQ=9)
                                                                            ggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                       <  1:5708334/94‑1 (MQ=255)
                                                                              ggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGaa                     <  1:2459546/96‑1 (MQ=255)
                                                                                cgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                   <  2:2452049/98‑1 (MQ=16)
                                                                                 gctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAggg                  >  2:433399/3‑101 (MQ=16)
                                                                                  ctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                 <  1:1727542/100‑1 (MQ=255)
                                                                                    gCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCGTCTGGCCCGAAGGGCaa               >  2:5195249/1‑101 (MQ=255)
                                                                                     ccccTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg              >  1:1362150/1‑101 (MQ=16)
                                                                                      cccTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAACGGCAAGc             <  2:1178357/101‑1 (MQ=16)
                                                                                        cTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg           >  1:2585213/1‑101 (MQ=16)
                                                                                            tGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAgatcggaa       >  2:5706634/1‑94 (MQ=16)
                                                                                            tGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc       <  2:3108288/101‑1 (MQ=19)
                                                                                             ggCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg      <  2:4529224/101‑1 (MQ=19)
                                                                                                cTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg   <  1:3521573/101‑1 (MQ=19)
                                                                                                cTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg   >  2:272152/1‑101 (MQ=255)
                                                                                                cTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTCTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg   <  2:5010192/101‑1 (MQ=255)
                                                                                                 tttCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg  <  1:213777/101‑1 (MQ=19)
                                                                                                 tttCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg  >  1:1449903/1‑101 (MQ=255)
                                                                                                        |                                                                                                  
GCTCCTCCCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCG.G.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGG  >  NC_014034/1696075‑1696275

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG