Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,348 A→G 9.6% intergenic (‑144/‑161) parE ← / → RCAP_rcc02080 DNA topoisomerase IV subunit B/HAD superfamily hydrolase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,249,3480A→G 9.6% 29.8 ‑2.3 ‑1.6145intergenic (‑144/‑161)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (4/10):  ref base (69/62):  total (73/72)
Fisher's exact test strand distribution p-value = 9.94e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.43e-02

Only 100 of 161 total aligned reads displayed.

GCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.G.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATT  >  NC_014034/2249249‑2249447
                                                                                                           |                                                                                                         
gCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGa                                                                                                           <  2:2066088/101‑1 (MQ=255)
 cTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaa                                                                                                          <  1:2927752/101‑1 (MQ=255)
 cTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaa                                                                                                          <  1:758475/101‑1 (MQ=255)
  ttGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaa                                                                                                         >  2:1538005/1‑101 (MQ=255)
   tGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTGATCTTGCGGaaaa                                                                                                        >  2:439733/1‑101 (MQ=38)
   tGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaaa                                                                                                        <  2:1298623/101‑1 (MQ=255)
   tGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaaaa                                                                                                        <  1:986743/101‑1 (MQ=255)
    gggCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                                       >  2:625976/1‑101 (MQ=255)
    gggCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                                       >  1:1337162/1‑101 (MQ=255)
    gggCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                                       <  2:4468/101‑1 (MQ=255)
     ggCACGCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCTGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAata                                                                                                      <  2:1455475/101‑1 (MQ=38)
      gCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAatat                                                                                                     <  2:2484105/101‑1 (MQ=255)
       cACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATc                                                                                                    <  2:2180251/101‑1 (MQ=255)
          ccGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCg                                                                                                 >  1:2346695/1‑101 (MQ=255)
           cGGGGCAAGAAGGGGGCTCT.GCCCCCGTGCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgg                                                                                                <  2:4371622/101‑1 (MQ=255)
           cGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgg                                                                                                >  1:1909028/1‑101 (MQ=255)
           cGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgg                                                                                                >  2:3468556/1‑101 (MQ=255)
             gggCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggg                                                                                              <  1:4246974/101‑1 (MQ=255)
               gCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggggg                                                                                            >  2:1158809/1‑101 (MQ=255)
               gCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggggg                                                                                            <  1:2587984/101‑1 (MQ=255)
               gCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggggg                                                                                            <  2:3245345/101‑1 (MQ=255)
                  agaggGGGCCTCT.CCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGa                                                                                         <  1:2338540/101‑1 (MQ=37)
                    aagGGGGCTCTCC.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.t                                                                                      <  2:951015/101‑1 (MQ=255)
                    aagGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.t                                                                                      >  2:1626727/1‑101 (MQ=255)
                    aagGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.t                                                                                      >  1:3570651/1‑101 (MQ=255)
                     agGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.Tg                                                                                     <  2:1526650/101‑1 (MQ=255)
                      gggggCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.Tgg                                                                                    <  1:513595/101‑1 (MQ=255)
                          gCTCT.GCCCCCGTCCGGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..TCTTTTCATGTTGCCGAAACTATCCCGGTCGCTGAA.TGgcggc                                                                                >  2:295579/1‑99 (MQ=21)
                          gCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCa                                                                                <  1:1006474/101‑1 (MQ=255)
                          gCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCa                                                                                <  2:3436257/101‑1 (MQ=255)
                             tctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATTAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gc                                                                            >  1:2175671/6‑101 (MQ=255)
                              ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcc                                                                           >  1:258922/5‑101 (MQ=255)
                              ctacccccGTCCCGTGCCGGGGCT.CCCCCGGGATATTTGATGCAAGATGAA...GGGG..GGGTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcc                                                                           <  1:51823/97‑1 (MQ=37)
                                 cccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCaga                                                                        <  1:4090986/100‑1 (MQ=255)
                                  ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGAA..GCAagag                                                                       >  2:1764339/1‑101 (MQ=255)
                                    ccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggg                                                                     >  2:929917/1‑101 (MQ=255)
                                       tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGGTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                  >  1:3163885/1‑101 (MQ=255)
                                       tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                  >  1:320558/1‑101 (MQ=255)
                                        cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                 >  1:2453951/1‑101 (MQ=255)
                                         ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                <  2:4291910/101‑1 (MQ=255)
                                         ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                <  1:3924492/101‑1 (MQ=255)
                                          cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                   >  2:1957836/8‑101 (MQ=21)
                                          cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                   <  1:1337573/94‑1 (MQ=18)
                                          cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                   <  2:3398244/94‑1 (MQ=18)
                                          cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                   <  2:272454/94‑1 (MQ=18)
                                            caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                                 <  1:2900124/96‑1 (MQ=18)
                                            caagcggACTCCCCCCGGGATATTGGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                                 <  1:3882924/96‑1 (MQ=18)
                                             gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c.g.c                                                         >  2:2978001/1‑101 (MQ=255)
                                              ccGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cc                                                        <  1:499750/101‑1 (MQ=255)
                                                gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccc                                                      >  2:508026/1‑101 (MQ=255)
                                                gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccc                                                      >  2:630121/1‑101 (MQ=255)
                                                gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccc                                                      >  1:421189/1‑101 (MQ=255)
                                                cggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.c.g                                                           <  2:2752959/100‑1 (MQ=18)
                                                 ggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.ccccc                                                     >  1:4337870/1‑101 (MQ=255)
                                                  gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccccc                                                    <  2:2566420/101‑2 (MQ=255)
                                                   gcT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccccct                                                   <  2:2654579/100‑3 (MQ=255)
                                                   actccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAG.C.G.ccc                                                       <  1:2070380/97‑1 (MQ=17)
                                                   aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccccct                                                   >  1:3225093/1‑99 (MQ=255)
                                                   aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccccct                                                   >  1:9319/1‑99 (MQ=255)
                                                       cccccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.C.G.CCCCCt                                                    <  1:3122854/101‑1 (MQ=0)
                                                       cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.ccccccttcg                                                <  1:591924/101‑6 (MQ=255)
                                                        ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTcgc                                                <  2:1700602/101‑1 (MQ=255)
                                                          ccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCgcgg                                              >  2:863530/1‑101 (MQ=255)
                                                          ccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCgcgg                                              >  1:1603494/1‑101 (MQ=255)
                                                           cGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCgcggg                                             >  2:1215378/1‑101 (MQ=255)
                                                            gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCgcgggc                                            >  1:3009782/1‑101 (MQ=255)
                                                             ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCgcgggcg                                           >  1:3601854/1‑101 (MQ=255)
                                                                 atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCcg                                       <  1:1026188/101‑1 (MQ=255)
                                                                   ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGTAG.C.GCCCCCCTCCGCGCGCGGCcgcc                                     >  2:1098390/1‑100 (MQ=25)
                                                                       ggCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGAcgc                                 <  2:144875/101‑1 (MQ=255)
                                                                          aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGCGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGAcgcgcg                              >  2:671280/1‑101 (MQ=255)
                                                                          aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGAcgcgcg                              >  2:434029/1‑101 (MQ=255)
                                                                            gaggAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGcgcc                            <  2:1825631/98‑1 (MQ=255)
                                                                             aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGTCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGcgccg                           <  2:421738/101‑1 (MQ=255)
                                                                             aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGcggcg                           >  2:539962/1‑101 (MQ=255)
                                                                             aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGcgccg                           <  2:4337601/101‑1 (MQ=255)
                                                                               gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CAC.CCCCCTTCGCGGGCGACCGCGACGCGCGCCGCa                         >  2:4016059/1‑101 (MQ=37)
                                                                               gAA...GGGG..GGCTTTCATATTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa                         >  2:485481/1‑101 (MQ=255)
                                                                                aggaggggg..gcCTTTCATCTGGCCGAAAATTTCCGGGGGGGGGAA.TGGCGCA..GCCAGGGGGGGGCAGCC.C.CCCCCTTCGCGGGCGGCCGCGACGCGcgccg                           <  1:1626833/96‑1 (MQ=255)
                                                                                 a...GGGG..GGCTTTCATCTTGCGGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                       <  2:1976457/101‑1 (MQ=255)
                                                                                 a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                       <  2:588497/101‑1 (MQ=255)
                                                                                 a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                       <  2:2410089/101‑1 (MQ=255)
                                                                                 a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                       <  1:1887955/101‑1 (MQ=255)
                                                                                 a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                       >  2:1113979/1‑101 (MQ=255)
                                                                                      ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                     >  1:1937042/1‑101 (MQ=255)
                                                                                           ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc                  <  2:4231711/101‑1 (MQ=255)
                                                                                           ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.CGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc                  <  2:110667/101‑1 (MQ=255)
                                                                                           ggCTTTCATCTTGCCAAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCC.C.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc                  <  1:83749/101‑1 (MQ=38)
                                                                                            ggtttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCg                 <  2:2138991/99‑1 (MQ=255)
                                                                                              tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg               >  2:3201805/1‑101 (MQ=255)
                                                                                              tttCATCTTGCCAAAAATGGCCGGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg               <  1:956221/101‑1 (MQ=37)
                                                                                               ttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCGGCAACAGGGCGTGc              >  2:1858893/1‑101 (MQ=25)
                                                                                                  aTCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcccc           >  1:991653/1‑101 (MQ=255)
                                                                                                    cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt         <  2:239874/101‑1 (MQ=255)
                                                                                                     ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt        >  2:477158/1‑101 (MQ=255)
                                                                                                     ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt        <  1:2801147/101‑1 (MQ=255)
                                                                                                         cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa    <  2:4034098/101‑1 (MQ=255)
                                                                                                          gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt   >  1:1560187/1‑101 (MQ=255)
                                                                                                          gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt   >  2:2391509/1‑101 (MQ=255)
                                                                                                           aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt  >  2:3879374/1‑101 (MQ=255)
                                                                                                           |                                                                                                         
GCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.G.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATT  >  NC_014034/2249249‑2249447

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