Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,374 +G 2.5% intergenic (‑170/‑135) 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,3742.→G 2.5% 10.8 ‑0.4 ‑0.1237intergenic (‑170/‑135)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (4/2):  ref base (129/102):  total (133/104)
Fisher's exact test strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.91e-01

Only 100 of 263 total aligned reads displayed.

CTCTGCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGC  >  NC_014034/2249276‑2249473
                                                                                                          |                                                                                                        
tctgcccccGTCCCGTGCCGCGACT.CTCCCGGGATATATGAAGCAAGAAGAA...GGGG..GCCTTACATCTTGCCGAAAATATCCCGGGGCGTGAA.TGGCGCA..Gc                                                                                                       <  2:5020550/96‑1 (MQ=25)
 ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..Gcc                                                                                                      >  1:3524313/5‑101 (MQ=255)
   gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGAAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATACCGGGGGGGGAA.GGGCGCA..GCCag                                                                                                    >  1:2491384/3‑101 (MQ=38)
         gtcccgtGCCGGGACT.CCCCCGGGTTATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                              >  2:3670845/1‑101 (MQ=255)
         gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                              >  1:3213276/1‑101 (MQ=255)
         gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                              >  1:4881073/1‑101 (MQ=255)
          tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                             >  2:4829050/1‑101 (MQ=255)
          tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                             <  2:4118106/101‑1 (MQ=255)
           cccgtGCCGGGACT.CCCCGGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCGCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                            <  2:107116/101‑1 (MQ=255)
           cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGCGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                            >  2:2443535/1‑101 (MQ=255)
            ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                           <  2:608220/101‑1 (MQ=255)
            ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                           <  2:3319237/101‑1 (MQ=255)
            ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GACAGAGGGGGGCa                                                                                           <  1:2728576/101‑1 (MQ=255)
             cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                          <  1:3312132/101‑1 (MQ=255)
             cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                          <  1:3498506/101‑1 (MQ=255)
             cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                          <  1:675225/101‑1 (MQ=255)
              gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                                                                        >  1:2566735/1‑101 (MQ=255)
              gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                                                                        >  1:4454983/1‑101 (MQ=255)
               tGCCTGGACT.CCCCCGGGATATTTGAGGCAAGATCAG...GGGG..GGTTTTCATCTTTCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                                                       <  1:2240441/101‑1 (MQ=37)
               caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                                                            <  2:902873/96‑1 (MQ=18)
                gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCAg.ag.c                                                                                     >  2:1737242/1‑101 (MQ=255)
                  gcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.c                                                                                        <  1:2286656/99‑1 (MQ=18)
                  cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc                                                                                   >  1:1052733/1‑101 (MQ=255)
                     gactccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAG.CG.cc                                                                                    <  1:4232038/96‑1 (MQ=17)
                     gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.ccc                                                                               <  1:2807368/101‑2 (MQ=255)
                     gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.ccc                                                                               <  2:270951/101‑2 (MQ=255)
                      aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.ccct                                                                              >  1:485248/1‑99 (MQ=255)
                      aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.ccct                                                                              >  1:2140057/1‑99 (MQ=255)
                        t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.cccttc                                                                            <  2:1645830/101‑5 (MQ=255)
                        t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.cccttc                                                                            <  1:2181228/101‑5 (MQ=255)
                        t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc.cccttc                                                                            >  1:4658426/1‑97 (MQ=255)
                           ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTcgc                                                                          <  2:4298079/101‑1 (MQ=255)
                            cccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCC.CC.tgc                                                                            <  2:1204802/101‑3 (MQ=0)
                              cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCC.CC.tgccc                                                                          >  2:3195667/1‑97 (MQ=2)
                               gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCgcgggc                                                                      >  1:568043/1‑101 (MQ=255)
                                ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCgcggggg                                                                     >  1:3003233/1‑99 (MQ=255)
                                 gATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCC..GCCAGAGGGGGGGAG.CGCCCC.CC.TTCgccggcgc                                                                    >  2:1464777/1‑100 (MQ=25)
                                    atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCcg                                                                 >  2:1375132/1‑101 (MQ=255)
                                     tttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCcgc                                                                <  2:4880746/101‑1 (MQ=255)
                                      ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGCGGGGGGCAG.CGCCCC.CC.TTAGCGGGCGGCcgcg                                                               >  1:4345644/1‑101 (MQ=37)
                                        gAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGAc                                                             >  2:4379902/1‑101 (MQ=255)
                                        gAGGAAAGATGAA...GGGG..GGCTTTCATTTTGCCGAAAATTTCGGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGAc                                                             <  2:1789784/101‑1 (MQ=37)
                                         aggtaggaggAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGAcg                                                            <  2:1582857/92‑1 (MQ=25)
                                            cAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGAcgcgc                                                         <  2:341805/101‑1 (MQ=255)
                                            cAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCACC.TTCGCGGGCGGCCGCGAcgcgc                                                         >  2:1280495/1‑101 (MQ=255)
                                               gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCC.CCATTCGCGGGCGGCCGCGACGCGcgcc                                                      >  2:225834/1‑101 (MQ=255)
                                                aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGcgccg                                                     >  1:2085237/1‑101 (MQ=255)
                                                 tGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGcgcagc                                                    >  2:2760495/1‑101 (MQ=255)
                                                  gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                   >  2:2560414/1‑101 (MQ=255)
                                                  gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                   <  2:4682148/101‑1 (MQ=255)
                                                  gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                   >  2:515019/1‑101 (MQ=255)
                                                  gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                   <  2:540454/101‑1 (MQ=255)
                                                  gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                   >  2:2673307/1‑101 (MQ=255)
                                                  gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                   >  1:3437886/1‑101 (MQ=255)
                                                  gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                   >  1:1613799/1‑101 (MQ=255)
                                                   aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                  >  1:1143316/1‑101 (MQ=255)
                                                    a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                 >  2:1840746/1‑101 (MQ=255)
                                                    a...GGGG..GGCTTTCATCTTGCCAAAAATGTCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                 <  2:396183/101‑1 (MQ=37)
                                                        gggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                <  2:2897380/101‑1 (MQ=255)
                                                        gggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                <  1:4806444/101‑1 (MQ=255)
                                                               gCTTTCATCATGCCGAAAATATCCCGGGGGGTGAA.TGGCCCA..GCCAGGGGGGGGAAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGCCg                                           >  2:1677257/1‑101 (MQ=37)
                                                                cTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGt                                          <  2:3826913/101‑1 (MQ=255)
                                                                 tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                         >  1:2365250/1‑101 (MQ=255)
                                                                 tttCATCTTGCCGAAAATATCCCGGGGGCTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                         >  2:398831/1‑101 (MQ=255)
                                                                   tCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                       >  2:1030871/1‑101 (MQ=255)
                                                                   tCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                       >  1:2075832/1‑101 (MQ=255)
                                                                    cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                      <  1:3892216/101‑1 (MQ=255)
                                                                     aTCTTGCCGAAAATATCCCGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcccc                                     <  1:3981074/101‑1 (MQ=255)
                                                                      tCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                    <  1:2859550/101‑1 (MQ=255)
                                                                           ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                               <  2:4216837/101‑1 (MQ=255)
                                                                              aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                            >  1:3784659/1‑101 (MQ=255)
                                                                               aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGCCGCGCGCCGCAACAGGCCGTGCCCCTTTGTCATTc                           >  1:429738/1‑101 (MQ=37)
                                                                               aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                           >  1:4295762/1‑101 (MQ=255)
                                                                               aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                           <  2:308107/101‑1 (MQ=255)
                                                                               aaaTATCCCGGGGGGTGAA.TGGCGCA..GACAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                           <  2:2399495/101‑1 (MQ=255)
                                                                                  tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                        >  2:3526695/1‑101 (MQ=255)
                                                                                   atCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                       >  2:2879455/1‑101 (MQ=255)
                                                                                   atCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                       <  1:4945884/101‑1 (MQ=255)
                                                                                   acccgggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                       <  1:1874073/96‑1 (MQ=255)
                                                                                     cccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGc                     >  1:600989/1‑101 (MQ=255)
                                                                                       gggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                   <  2:1009809/100‑1 (MQ=255)
                                                                                        ggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                  >  2:3737446/1‑101 (MQ=255)
                                                                                           gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTGGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt               >  1:1550806/1‑101 (MQ=255)
                                                                                           gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt               >  1:4000052/1‑101 (MQ=255)
                                                                                           gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt               >  1:3919276/1‑101 (MQ=255)
                                                                                           gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt               >  2:4470104/1‑101 (MQ=255)
                                                                                            ggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg              >  1:707479/1‑101 (MQ=255)
                                                                                              ggAA.TGGCGCA..GCCAGAGGGGGGCAGCCC.CCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc            <  2:3478431/100‑1 (MQ=255)
                                                                                               gAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa           >  2:2122268/1‑101 (MQ=255)
                                                                                               gAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGACCCTTTGTCATTCTGTCGCACCCATGTCa           >  2:824026/1‑101 (MQ=255)
                                                                                                aa.TGGCGCA..GCCAGGGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt          <  1:4564322/101‑1 (MQ=255)
                                                                                                aa.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt          <  2:2274322/101‑1 (MQ=255)
                                                                                                aa.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGCCATTCTGTCGCACCCATGTCAt          >  1:2320318/1‑101 (MQ=255)
                                                                                                 a.TGGCGAA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGGCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa         <  2:2491210/101‑1 (MQ=255)
                                                                                                   tGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa        >  1:4413198/1‑101 (MQ=255)
                                                                                                   tGGCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa        <  2:2415479/101‑1 (MQ=255)
                                                                                                    ggCGCA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg       <  1:758587/101‑1 (MQ=255)
                                                                                                       gcA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg    <  1:2936855/101‑1 (MQ=255)
                                                                                                        cA..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg   >  1:3135105/1‑101 (MQ=255)
                                                                                                         a..GCCAGAGGGGGGCAG.CGCCCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGc  >  2:4603970/1‑101 (MQ=255)
                                                                                                          |                                                                                                        
CTCTGCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCC.CC.TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGC  >  NC_014034/2249276‑2249473

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG