Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,333 G→C 6.4% intergenic (‑129/‑176) 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,3330G→C 6.4% 14.8 ‑2.9 ‑0.0160intergenic (‑129/‑176)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (2/9):  ref base (75/73):  total (77/83)
Fisher's exact test strand distribution p-value = 5.78e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.66e-01

Only 100 of 176 total aligned reads displayed.

TGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.G.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTC  >  NC_014034/2249234‑2249431
                                                                                                           |                                                                                                        
tGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...gggg..gg                                                                                                          <  1:3768299/101‑1 (MQ=255)
tGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...gggg..gg                                                                                                          >  1:4236405/1‑101 (MQ=255)
 gCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGc                                                                                                         <  2:2203615/101‑1 (MQ=255)
 gCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGc                                                                                                         <  1:431239/101‑1 (MQ=255)
    caGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCttt                                                                                                      <  2:3655832/101‑1 (MQ=255)
     aGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTc                                                                                                     <  2:778834/101‑1 (MQ=255)
     aGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTc                                                                                                     <  1:1594615/101‑1 (MQ=255)
     aGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTc                                                                                                     <  1:4216614/101‑1 (MQ=255)
      gATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCa                                                                                                    >  2:2338390/1‑101 (MQ=255)
       aTACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCACGATGAA...GGGG..GGCTTTCAt                                                                                                   >  1:2654751/1‑101 (MQ=255)
         aCCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCt                                                                                                 >  1:907763/1‑101 (MQ=255)
         aCCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCt                                                                                                 >  1:588565/1‑101 (MQ=255)
          ccGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCtt                                                                                                <  1:1894832/101‑1 (MQ=255)
          ccGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCtt                                                                                                <  1:1806711/101‑1 (MQ=255)
           cGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTg                                                                                               <  2:1451728/101‑1 (MQ=255)
            gCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGc                                                                                              <  1:1359721/101‑1 (MQ=255)
              aGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCg                                                                                            <  2:4228665/101‑1 (MQ=255)
              aGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCg                                                                                            >  1:239901/1‑101 (MQ=255)
              aGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCg                                                                                            <  2:905584/101‑1 (MQ=255)
              aGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCg                                                                                            >  1:3377750/1‑101 (MQ=255)
                cTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaa                                                                                          <  1:758475/101‑1 (MQ=255)
                cTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGaa                                                                                          <  1:2927752/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)
                   gggCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                       >  1:1337162/1‑101 (MQ=255)
                   gggCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAat                                                                                       >  2:625976/1‑101 (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:3232114/1‑101 (MQ=255)
                          cGGGGCAAGAAGGGGGCTCT.GCCCCCGTGCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgg                                                                                <  2:4371622/101‑1 (MQ=255)
                             ggCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCggggg                                                                             >  2:4366761/1‑101 (MQ=255)
                              gCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggggg                                                                            <  1:2587984/101‑1 (MQ=255)
                              gCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCgggggg                                                                            >  2:1158809/1‑101 (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)
                                     gggggCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.Tgg                                                                    <  1:513595/101‑1 (MQ=255)
                                     gggggCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.Tgg                                                                    >  2:3739028/1‑101 (MQ=255)
                                      ggggCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgc                                                                   <  1:69665/101‑1 (MQ=255)
                                         gCTCT.GCCCCCGTCCGGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..TCTTTTCATGTTGCCGAAACTATCCCGGTCGCTGAA.TGgcggc                                                                >  2:295579/1‑99 (MQ=21)
                                         gCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCa                                                                <  2:3436257/101‑1 (MQ=255)
                                         gCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCa                                                                <  1:1006474/101‑1 (MQ=255)
                                             ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcc                                                           <  2:452480/97‑1 (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)
                                                  cccggcccgtGCCGGGACT.CCCCGGGGATATTTGAGGCAAGATGAA...GGGG..GGGTTTCATCTTGCCGAAAATATCGCGGGGGGTGAA.TGGCGCA..GCCAGAgg                                                      <  1:2938630/101‑1 (MQ=38)
                                                      tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGGTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                  >  1:3163885/1‑101 (MQ=255)
                                                      tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                  >  1:320558/1‑101 (MQ=255)
                                                      tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                  >  2:2576331/1‑101 (MQ=255)
                                                       cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                 >  1:2453951/1‑101 (MQ=255)
                                                       cccggGCCGGGCCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                 <  1:1570400/101‑1 (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)
                                                         cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                               <  2:1010498/101‑1 (MQ=255)
                                                         cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                   >  2:1957836/8‑101 (MQ=21)
                                                         cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                   <  2:3398244/94‑1 (MQ=18)
                                                           caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                 <  1:2900124/96‑1 (MQ=18)
                                                           caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                 <  1:2327243/96‑1 (MQ=18)
                                                           caagcggACTCCCCCCGGGATATTGGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                 <  1:3882924/96‑1 (MQ=18)
                                                             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:110774/1‑101 (MQ=255)
                                                               gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccc                                      >  2:508026/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)
                                                                  actccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAG.C.G.ccc                                       <  1:2070380/97‑1 (MQ=17)
                                                                   cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccccctt                                  >  2:1921007/1‑98 (MQ=255)
                                                                   actacgccGGGATATTTGGGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATCTCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.T.cccccct                                   <  2:1125168/95‑3 (MQ=25)
                                                                     ccccccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.C.G.ccccc                                     <  2:2095384/100‑1 (MQ=0)
                                                                      cccccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.C.G.CCCCCt                                    <  1:3122854/101‑1 (MQ=0)
                                                                      cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.ccccccttcg                                <  2:3163756/101‑6 (MQ=255)
                                                                       ccccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.C.G.CCCCCtg                                   >  2:670123/1‑100 (MQ=2)
                                                                          cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.C.G.CCCCCtgccc                                >  2:2044605/1‑97 (MQ=2)
                                                                           gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCgcgggc                            >  2:2718678/1‑101 (MQ=255)
                                                                           gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCgcgggc                            >  1:954626/1‑101 (MQ=255)
                                                                            ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCgcgggcg                           >  1:3601854/1‑101 (MQ=255)
                                                                             gATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCACGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCgccggcgg                          >  1:3619563/1‑101 (MQ=25)
                                                                              atgtttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGggcggc                         <  2:1725741/98‑1 (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)
                                                                                        cAAGATGAA...GGTG..GTCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGAcgcgc               <  2:488969/101‑1 (MQ=25)
                                                                                         aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGCGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGAcgcgcg              >  2:671280/1‑101 (MQ=255)
                                                                                         aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGAcgcgcg              >  2:1480789/1‑101 (MQ=255)
                                                                                         aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGAcgcgcg              >  2:1578778/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.GCCCCCCTTCGCGGGCGGCCGCGACGCGcgccg           >  1:981528/1‑101 (MQ=255)
                                                                                            aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGcgccg           <  2:4337601/101‑1 (MQ=255)
                                                                                            aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGcgccg           >  1:4210006/1‑101 (MQ=255)
                                                                                              gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CAC.CCCCCTTCGCGGGCGACCGCGACGCGCGCCGCa         >  2:4016059/1‑101 (MQ=37)
                                                                                               aggaggggg..gcCTTTCATCTGGCCGAAAATTTCCGGGGGGGGGAA.TGGCGCA..GCCAGGGGGGGGCAGCC.C.CCCCCTTCGCGGGCGGCCGCGACGCGcgccg           <  1:1626833/96‑1 (MQ=255)
                                                                                               aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCaa        <  2:1254047/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:588497/101‑1 (MQ=255)
                                                                                                a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc       <  2:2410089/101‑1 (MQ=255)
                                                                                                     ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg     >  1:1937042/1‑101 (MQ=255)
                                                                                                     ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg     <  1:3468796/101‑1 (MQ=255)
                                                                                                          ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc  <  2:4231711/101‑1 (MQ=255)
                                                                                                          ggCTTTCATCTTGCCAAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAGCC.C.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc  <  1:83749/101‑1 (MQ=38)
                                                                                                           |                                                                                                        
TGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.G.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTC  >  NC_014034/2249234‑2249431

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