Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,379 G→A 2.2% intergenic (‑175/‑130) 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,3790G→A 2.2% 3.9 ‑0.4 ‑0.1229intergenic (‑175/‑130)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (4/1):  ref base (132/92):  total (136/93)
Fisher's exact test strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.61e-01

Only 100 of 270 total aligned reads displayed.

TGCCCCGTCCCGTGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAT  >  NC_014034/2249279‑2249478
                                                                                                            |                                                                                                     
gcccccGTCCCGTGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                                       >  2:3855116/3‑101 (MQ=255)
 cccccGTCCCGTGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCaga                                                                                                      <  2:2317432/100‑1 (MQ=255)
      gtcccgtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                                 >  1:2258931/1‑101 (MQ=255)
       tcccgtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                                <  2:4370485/101‑1 (MQ=255)
        cccgtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                               <  1:3253733/101‑1 (MQ=255)
         ccgtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                              >  1:3778898/1‑101 (MQ=255)
         ccgtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGcg                                                                                              >  1:1104936/1‑100 (MQ=255)
          cgtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                             <  2:5103290/101‑1 (MQ=255)
           gtGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                                                                           >  2:60989/1‑101 (MQ=255)
            tGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                                                          >  2:5691173/1‑101 (MQ=255)
            tGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                                                          >  2:5475188/1‑101 (MQ=255)
              ccGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cc                                                                                       >  2:2292392/1‑101 (MQ=255)
               cGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc                                                                                      >  1:1279654/1‑101 (MQ=255)
                gggACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccc                                                                                     >  1:3995097/1‑101 (MQ=255)
                 ggccTCCCCCGGGATATTGGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccc                                                                                    <  1:1024759/98‑1 (MQ=255)
                 ggACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccc                                                                                    >  2:4023660/1‑101 (MQ=255)
                 ggACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccc                                                                                    <  2:5644335/101‑1 (MQ=255)
                 gactccccccGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAG.CG.cc                                                                                       <  2:1610315/96‑1 (MQ=17)
                  gACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccccc                                                                                   >  1:1672024/1‑100 (MQ=255)
                    cTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccccctt                                                                                 <  1:3172465/101‑4 (MQ=255)
                       ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTcgc                                                                               >  2:5608327/1‑101 (MQ=255)
                          cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCCtgccc                                                                               >  2:4065635/1‑97 (MQ=2)
                           gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCGC.GG.CCCCCCTCgggggcg                                                                          >  1:2812371/1‑101 (MQ=37)
                           gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggc                                                                           <  2:4246316/101‑1 (MQ=255)
                           gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggc                                                                           >  1:4592779/1‑101 (MQ=255)
                             gATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCCtgccccgc                                                                            >  2:4463974/1‑94 (MQ=2)
                             gATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggcgg                                                                         <  1:5228203/101‑1 (MQ=255)
                               ttttggAGGAAAGAAAAG...GGGG..GGTTTTAATCTTGCCAAAAATATCCGGGGGGGTGAA.TGGCGCA..GCAAGGGGGGGGCAGCCC.CCCCCTTCGCGGGCGGcc                                                                       <  2:1175788/96‑1 (MQ=255)
                                 tttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCcgc                                                                     >  2:2914001/1‑101 (MQ=255)
                                  ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGccaga                                                                    >  2:5220671/1‑98 (MQ=255)
                                    gAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAc                                                                  >  1:5427229/1‑101 (MQ=255)
                                     aGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcg                                                                 >  2:4950824/1‑101 (MQ=255)
                                        cAAGATGAA...GGGG..GCCTTTCATCTGGCCGAAAATATGCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgc                                                              <  2:2175851/101‑1 (MQ=37)
                                         aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAccggcg                                                             >  1:3044340/1‑96 (MQ=25)
                                           gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgcc                                                           >  1:4795376/1‑101 (MQ=255)
                                                a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                      <  1:2572678/101‑1 (MQ=255)
                                                a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                      <  2:1951910/101‑1 (MQ=255)
                                                 gatggtg..ggCTTTCAACTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                        <  2:744295/95‑1 (MQ=25)
                                                     ggg..ggTTTTCATCTTGCCGAAAAAACCCGGGGGGGTGAA.TGGCGCA..GCCAGGGGGGGCCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                    <  1:4228043/101‑1 (MQ=37)
                                                     ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                    <  1:4987746/101‑1 (MQ=255)
                                                          ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc                                                 <  2:5398942/101‑1 (MQ=255)
                                                            cTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGt                                               <  1:1916433/101‑1 (MQ=255)
                                                            cTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGt                                               <  1:4143724/101‑1 (MQ=255)
                                                             tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                              >  1:2357232/1‑101 (MQ=255)
                                                             tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                              <  1:1366222/101‑1 (MQ=255)
                                                             tttCATCTTGCCAAAAATACCGGGGGGGGTAAA.GGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                              <  1:61000/101‑1 (MQ=37)
                                                              ttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                                             >  2:4698769/1‑101 (MQ=255)
                                                                 aTCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcccc                                          <  1:1243856/101‑1 (MQ=255)
                                                                 aTCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcccc                                          >  1:34726/1‑101 (MQ=255)
                                                                  tCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                         >  1:4329083/1‑101 (MQ=255)
                                                                  tCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                         <  2:2237122/101‑1 (MQ=255)
                                                                   cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                        >  2:1313696/1‑101 (MQ=255)
                                                                    ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCGG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCGGCAACAGGTCGTGCCCCttt                                       >  2:2555882/1‑101 (MQ=37)
                                                                    ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                       <  2:3581247/101‑1 (MQ=255)
                                                                    ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                       >  1:5744331/1‑101 (MQ=255)
                                                                     tGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                      <  2:4277762/101‑1 (MQ=255)
                                                                      caaaaaaaTCTGCGGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCGGCCC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGt                                     <  2:1037156/97‑1 (MQ=255)
                                                                       ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                                    <  2:1861347/101‑1 (MQ=255)
                                                                         gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                  <  1:3739804/101‑1 (MQ=255)
                                                                         gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                  >  1:269632/1‑101 (MQ=255)
                                                                         gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                  >  2:422137/1‑101 (MQ=255)
                                                                         gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                  <  2:374572/101‑1 (MQ=255)
                                                                          aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGTGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                 >  1:5656514/1‑101 (MQ=255)
                                                                           aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                >  1:4728132/1‑101 (MQ=255)
                                                                           aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                >  1:3514480/1‑101 (MQ=255)
                                                                            aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                               >  2:4098105/1‑101 (MQ=255)
                                                                             atatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGACGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTg                              >  2:341770/1‑101 (MQ=255)
                                                                              tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                             <  1:1828752/101‑1 (MQ=255)
                                                                               atCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                            <  1:885779/101‑1 (MQ=255)
                                                                                 cccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGc                          >  1:2011492/1‑101 (MQ=255)
                                                                                  ccGGGGGGTGAA.TGGCGCA..GCCAGGGGAGGTCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa                         <  2:3687640/101‑1 (MQ=255)
                                                                                   cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                        >  2:2127671/1‑101 (MQ=255)
                                                                                   cGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAGCCC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                        <  1:1499431/101‑1 (MQ=255)
                                                                                     gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                      >  2:2077836/1‑101 (MQ=255)
                                                                                     gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGCTGCAccc                      >  2:4867085/1‑101 (MQ=38)
                                                                                      ggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                     >  1:559764/1‑101 (MQ=255)
                                                                                      ggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                     >  1:4904079/1‑101 (MQ=255)
                                                                                      ggggGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                     <  2:3228511/101‑1 (MQ=255)
                                                                                       gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                    >  1:1050514/1‑101 (MQ=255)
                                                                                        ggTGAA.TGGCGCA..GCCGGGGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                   <  2:1104861/101‑1 (MQ=255)
                                                                                        ggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                   <  2:2605921/101‑1 (MQ=255)
                                                                                        ggTGAA.TGGCCCA..CCCGGGGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                   <  1:5608619/101‑1 (MQ=255)
                                                                                         gTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                  >  2:3496758/1‑101 (MQ=255)
                                                                                          ggAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc                 <  1:2174160/100‑1 (MQ=255)
                                                                                            aa.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt               >  1:1424690/1‑101 (MQ=255)
                                                                                            aa.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCACTCTGTCGCACCCATGTCAt               <  1:3553284/101‑1 (MQ=255)
                                                                                             a.TGGCGCA..GCCAGAGGGGGGCGG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa              <  2:4647387/101‑1 (MQ=255)
                                                                                                ggCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg            >  2:1006920/1‑101 (MQ=255)
                                                                                                ggCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg            >  2:2675298/1‑101 (MQ=255)
                                                                                                ggCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg            <  1:2675447/101‑1 (MQ=255)
                                                                                                   gcA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg         <  2:4537448/101‑1 (MQ=255)
                                                                                                   gcA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg         >  1:5488282/1‑101 (MQ=255)
                                                                                                    cA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCGCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg        >  2:5005509/1‑101 (MQ=255)
                                                                                                        cccAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTCGcc      >  2:1714662/2‑101 (MQ=255)
                                                                                                        atctgagGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGcc      <  2:4355326/97‑1 (MQ=255)
                                                                                                          cAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc    >  2:3110152/1‑101 (MQ=255)
                                                                                                          cAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc    >  2:276084/1‑101 (MQ=255)
                                                                                                           agagGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa   >  1:759998/1‑101 (MQ=255)
                                                                                                            gagGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGccagat  >  1:4355592/1‑97 (MQ=255)
                                                                                                            gagGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt  <  2:3778655/101‑1 (MQ=255)
                                                                                                            |                                                                                                     
TGCCCCGTCCCGTGCCGGGACTCCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAT  >  NC_014034/2249279‑2249478

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG