Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,373 C→A 5.9% intergenic (‑169/‑136) 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,3730C→A 5.9% 15.4 ‑1.1 ‑1.7170intergenic (‑169/‑136)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (7/3):  ref base (80/80):  total (87/83)
Fisher's exact test strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.78e-02

Only 100 of 195 total aligned reads displayed.

GGCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGG  >  NC_014034/2249274‑2249472
                                                                                                            |                                                                                                        
ggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGgcgc                                                                                                          <  2:1721042/101‑1 (MQ=255)
  ctctGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..g                                                                                                      <  1:106933/101‑1 (MQ=255)
   tctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..Gc                                                                                                     >  1:4224026/6‑101 (MQ=255)
   tctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGCGCA..Gc                                                                                                     >  1:811589/6‑101 (MQ=255)
    ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..Gcc                                                                                                    <  1:1853884/97‑1 (MQ=255)
      gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCag                                                                                                  <  2:1149245/99‑1 (MQ=255)
      gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCag                                                                                                  <  1:2604459/99‑1 (MQ=255)
      gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGCGCA..GCCag                                                                                                  >  1:2243123/3‑101 (MQ=255)
         cccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAgg                                                                                               >  1:4129767/1‑101 (MQ=255)
         cccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAgg                                                                                               <  1:368489/101‑1 (MQ=255)
          ccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAggg                                                                                              <  2:2269323/101‑1 (MQ=255)
              ccgatctcGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                          <  2:425147/94‑1 (MQ=255)
              cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                          >  2:3872364/1‑101 (MQ=255)
               ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                         <  1:560713/101‑1 (MQ=255)
                cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.g                                                                                       <  2:3415444/101‑1 (MQ=255)
                cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.g                                                                                       >  1:84139/1‑101 (MQ=255)
                cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.g                                                                                       <  2:24645/101‑1 (MQ=255)
                cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.g                                                                                       >  2:4102191/1‑101 (MQ=255)
                cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCC.GGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                                            >  2:2990221/8‑101 (MQ=18)
                 gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.g.c                                                                                     <  1:1001720/101‑1 (MQ=255)
                   gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.g.cg.c                                                                                  >  1:3537406/1‑101 (MQ=255)
                     cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.gatcgg                                                                                  >  1:425204/1‑96 (MQ=255)
                     cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.ccc                                                                                >  1:1918090/1‑101 (MQ=255)
                     cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.ccc                                                                                >  1:2850463/1‑101 (MQ=255)
                       ggACT.CCCCCGGGATATTTGAGGCAAGATGAC...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.ccccc                                                                              >  1:3808942/1‑101 (MQ=255)
                        gactccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCC.GGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCA.G.CG.cc                                                                                 <  2:3320767/96‑1 (MQ=17)
                         aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGTGGGGCA.G.CG.cccccct                                                                            <  2:4129505/101‑3 (MQ=255)
                           t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.ccccccttc                                                                          >  1:1710117/1‑97 (MQ=255)
                           t.CCCCCGAGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.ccccccttc                                                                          >  2:1624239/1‑97 (MQ=255)
                            ccccccGGGTTTTTTGGGGGAAGAGGAA...GGGG..GGTTTTCATTTTGCCAAAAAAATCCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.cccccctt                                                                           <  1:1422384/100‑4 (MQ=21)
                             cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.ccccccttcg                                                                         <  2:3959378/101‑6 (MQ=255)
                             cccccGGGATATTTGAGGCAAGAGGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.ccccccttcg                                                                         <  2:450745/101‑6 (MQ=255)
                              ccccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC..AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCA.G.CG.CCCCCtg                                                                            >  2:1635175/1‑100 (MQ=0)
                              ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTcgc                                                                         <  1:1547977/101‑1 (MQ=255)
                              ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTcgc                                                                         >  1:2354331/1‑101 (MQ=255)
                               cccGGGATATTTGCAGCAAGGTGAAACCGGTT..GCCTTTCATCTGGCCCGAAATATCCC.GGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCA.G.CG.ccctcctgc                                                                          <  1:291216/101‑7 (MQ=14)
                               cccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTcgcg                                                                        <  2:3190971/101‑1 (MQ=255)
                               cccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAAGATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTcgcg                                                                        >  1:3579042/1‑101 (MQ=255)
                                ccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCgcgg                                                                       >  1:82437/1‑101 (MQ=255)
                                 cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGAGAAT.TGGCCAAAGGCCAAGAGGGGGCA.G.CG.CCCCCtgacc                                                                         >  2:1879142/1‑97 (MQ=2)
                                 cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC..AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCA.G.CG.CCCCCtgccc                                                                         >  2:1996119/1‑97 (MQ=2)
                                 cGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCgcggg                                                                      <  2:3162182/101‑1 (MQ=255)
                                  gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCG.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCgcgggc                                                                     <  2:2916758/101‑1 (MQ=255)
                                  gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCgcgggc                                                                     >  1:1763632/1‑101 (MQ=255)
                                  gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCgcgggc                                                                     >  2:3783320/1‑101 (MQ=255)
                                   gaatatTGGAAGAAAGATGAG...GGGG..TTTTTTCCTCGTGCAAAAAATCTCCG.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.GCCC.CCCCCTTCgcgggcg                                                                    <  1:762483/99‑1 (MQ=255)
                                      tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGcc                                                                 >  2:1689121/1‑101 (MQ=255)
                                       atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTCCGCGGGCGGCcg                                                                >  1:991321/1‑101 (MQ=255)
                                             ggCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGAcgc                                                          >  2:370932/1‑101 (MQ=255)
                                              gCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGAcgcg                                                         <  2:1112260/101‑1 (MQ=255)
                                               cAAGATGAA...GGGG..GGCTTTCATCTTGCCGAGAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgc                                                        >  1:3765924/1‑101 (MQ=255)
                                                aaGATGAG...GGGG..GGCTTTCATCTTGCCGAAAATAGCGC.GGGGGGTGAA.TGGCGCA..GCGAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg                                                       <  2:1066977/101‑1 (MQ=37)
                                                aaGAGGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg                                                       <  1:2440432/101‑1 (MQ=255)
                                                   aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgccg                                                    >  2:3618452/1‑101 (MQ=255)
                                                     gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                  <  1:1461113/101‑1 (MQ=255)
                                                       a...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                <  1:1899413/101‑1 (MQ=255)
                                                       a...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                <  2:1817901/101‑1 (MQ=255)
                                                          ggggg..ggTTTTCTTCTTGCGAAAAATTCCGG.GGGGGGTAAA.TGGCGCG..GCCGGGGGGGGGCG.GCCC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                <  1:3132242/100‑1 (MQ=255)
                                                           gggg..ggTTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                               <  2:899651/101‑1 (MQ=255)
                                                           gggg..ggCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                               >  2:1269296/1‑101 (MQ=255)
                                                                 ggCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc                                           >  2:2568943/1‑101 (MQ=255)
                                                                  gCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCg                                          <  1:3207496/101‑1 (MQ=255)
                                                                  gCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCg                                          <  1:3527858/101‑1 (MQ=255)
                                                                    tttCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                        >  1:1186225/1‑101 (MQ=255)
                                                                      tCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTTGTGcc                                      >  2:2107527/1‑101 (MQ=255)
                                                                       cATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                     <  2:136097/101‑1 (MQ=255)
                                                                         tCTTGCCGAAAATATCCG.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                   <  1:2631767/101‑1 (MQ=255)
                                                                          cTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGCCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                  >  1:2265057/1‑101 (MQ=255)
                                                                           ttGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCGGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCTCCttt                                 <  2:1677557/101‑1 (MQ=37)
                                                                           ttGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                 <  1:1265657/101‑1 (MQ=255)
                                                                            tggcGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                >  1:2420108/4‑101 (MQ=255)
                                                                            tGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                >  1:1651336/1‑101 (MQ=255)
                                                                            tGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                >  2:2885210/1‑101 (MQ=255)
                                                                            tGCCAAAAATATCCG.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                <  2:2924812/101‑1 (MQ=38)
                                                                            ggcagAAAATACCCG.GGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                <  2:811511/97‑1 (MQ=37)
                                                                              ccGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                              >  1:3153813/1‑101 (MQ=255)
                                                                               cGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                             <  2:1009263/101‑1 (MQ=255)
                                                                                 aaaaTATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGGGCCCCTTTGTCAtt                           >  1:1053140/1‑101 (MQ=255)
                                                                                         cc.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa                   <  1:4053584/101‑1 (MQ=255)
                                                                                          c.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                  <  2:2787310/101‑1 (MQ=255)
                                                                                            ggggggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                 <  1:208656/101‑1 (MQ=255)
                                                                                            ggggggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                 >  2:4010255/1‑101 (MQ=255)
                                                                                             gggggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                <  2:1917964/101‑1 (MQ=255)
                                                                                             gggggTAAA.GGGCGCA..GCCGGGGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                <  2:3553598/101‑1 (MQ=255)
                                                                                              ggggTGAA.TGGCGCA..GCCAGAGGGGGCCACC.CC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa               <  2:514208/101‑1 (MQ=255)
                                                                                               gggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt              <  2:982947/101‑1 (MQ=255)
                                                                                               gggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt              >  1:3355792/1‑101 (MQ=255)
                                                                                               gggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt              <  2:1706323/101‑1 (MQ=255)
                                                                                                ggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg             <  2:3299626/101‑1 (MQ=255)
                                                                                                 gTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt            >  2:148543/1‑101 (MQ=255)
                                                                                                 gTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt            <  1:3028008/101‑1 (MQ=255)
                                                                                                  tGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc           >  1:1124005/1‑101 (MQ=255)
                                                                                                   ttccgatctcA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc           <  2:4117340/92‑1 (MQ=255)
                                                                                                   gAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGACGCGCCGCGCGCCGCAACAGGCAGTGCCCCTTTGTCATTCTGTCGTCCCCCTGTCa          >  2:1900639/1‑101 (MQ=255)
                                                                                                     a.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa        <  1:3246538/101‑1 (MQ=255)
                                                                                                       tGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa       <  2:3818028/101‑1 (MQ=255)
                                                                                                       tGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa       <  1:675004/101‑1 (MQ=255)
                                                                                                       tGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa       <  2:289922/101‑1 (MQ=255)
                                                                                                          cgcA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCt    <  1:2475575/101‑1 (MQ=255)
                                                                                                           gcA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg   <  1:1251423/101‑1 (MQ=255)
                                                                                                            |                                                                                                        
GGCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGG  >  NC_014034/2249274‑2249472

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