Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,372 G→C 2.2% intergenic (‑168/‑137) 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,3720G→C 2.2% 3.9 ‑0.9 ‑1.3229intergenic (‑168/‑137)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (4/1):  ref base (125/99):  total (129/100)
Fisher's exact test strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.75e-02

Only 100 of 265 total aligned reads displayed.

GGCTCTG.CCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTG  >  NC_014034/2249274‑2249471
                                                                                                           |                                                                                                        
ggCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgcgc                                                                                                         <  2:1902201/101‑1 (MQ=255)
ggCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgcgc                                                                                                         <  1:3463851/101‑1 (MQ=255)
ggCTCTGCCCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgcgc                                                                                                         >  2:5143000/1‑101 (MQ=255)
 gctctccccc.cGTCCCGTGCCGGGACC.CCCCCGGGATTTTTGAGGAAAGAGGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCa                                                                                                        <  1:169553/94‑1 (MQ=37)
    ctG.CCCGCGTCCCGGGCCGGGACT.CCCGCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcc                                                                                                   <  1:3555305/101‑1 (MQ=37)
     tgcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCa                                                                                                  >  2:2893922/4‑101 (MQ=255)
      gcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCGAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcccg                                                                                                 >  2:3091455/3‑99 (MQ=255)
      gcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                                 >  2:3855116/3‑101 (MQ=255)
      gcccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                                 >  2:4361792/3‑101 (MQ=255)
       cccc.cGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCaga                                                                                                <  2:2317432/100‑1 (MQ=255)
         cc.cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgg                                                                                              <  2:4689248/101‑1 (MQ=255)
           ccatcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggg                                                                                             <  1:2896283/98‑1 (MQ=255)
            cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                                            >  2:251324/1‑101 (MQ=255)
              tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                          >  1:4342072/1‑101 (MQ=255)
                ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.Ca                                                                                       <  1:4941439/101‑1 (MQ=255)
                ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.Ca                                                                                       >  1:3778898/1‑101 (MQ=255)
                 cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg                                                                                      <  2:4629769/101‑1 (MQ=255)
                 cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg                                                                                      >  1:3071048/1‑101 (MQ=255)
                  gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg.c                                                                                    >  1:3575561/1‑101 (MQ=255)
                  gcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggggg                                                                                          >  2:4767926/7‑101 (MQ=18)
                   tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg.cg                                                                                   >  2:5475188/1‑101 (MQ=255)
                   tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg.cg                                                                                   >  2:5691173/1‑101 (MQ=255)
                    gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg.cg.c                                                                                 >  1:2606070/1‑101 (MQ=255)
                      cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccc                                                                               >  1:5282299/1‑101 (MQ=255)
                       gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.cccc                                                                              >  1:3995097/1‑101 (MQ=255)
                       gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.cccc                                                                              >  1:1861457/1‑101 (MQ=255)
                        ggccT.CCCCCGGGATATTGGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccccc                                                                             <  1:1024759/98‑1 (MQ=255)
                        ggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccccc                                                                             >  2:3739006/1‑101 (MQ=255)
                        ggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccccc                                                                             >  2:4051320/1‑101 (MQ=255)
                         gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.cccccc                                                                            >  1:1672024/1‑100 (MQ=255)
                           cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GTCTTTCATCTGGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.cccccctt                                                                          <  2:5105679/101‑4 (MQ=255)
                           cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.cccccctt                                                                          <  1:3172465/101‑4 (MQ=255)
                            t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccccccttc                                                                         >  1:3511685/1‑97 (MQ=255)
                            t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccccccttc                                                                         >  2:580828/1‑97 (MQ=255)
                            t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.ccccccttc                                                                         >  2:3809318/1‑97 (MQ=255)
                               ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTcgc                                                                        >  2:3213230/1‑101 (MQ=255)
                                cccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGG.CAG.CG.CCCCCtgc                                                                          >  1:1041015/1‑99 (MQ=2)
                                cccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTcgcg                                                                       >  2:3663376/1‑101 (MQ=255)
                                 ccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCgcgg                                                                      <  2:1759336/101‑1 (MQ=255)
                                   gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCgcgggc                                                                    <  1:3268213/101‑1 (MQ=255)
                                   gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCgcgggc                                                                    >  1:4592779/1‑101 (MQ=255)
                                   gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCgcgggc                                                                    >  1:2198455/1‑101 (MQ=255)
                                    ggATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGG.CAG.CG.CCCCCtgcaccg                                                                      <  2:4900346/101‑7 (MQ=0)
                                    ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCgcgggcg                                                                   >  1:1103886/1‑101 (MQ=255)
                                     gATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCgcgggcgg                                                                  >  1:4712468/1‑101 (MQ=255)
                                       ttttggAGGAAAGAAAAG...GGGG..GGTTTTAATCTTGCCAAAAATATCCGGGGGGGTGAA.TGGCGCA..GCAAGGGGGGGG.CAGCCC.CCCCCTTCGCGGGCGGcc                                                                <  2:1175788/96‑1 (MQ=255)
                                       tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAAGATCCGGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGcc                                                                <  1:1374223/101‑1 (MQ=37)
                                            gAGGCAAGATGAA...GGGG..TGCGTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGAc                                                           <  2:4309113/101‑1 (MQ=25)
                                                cAAGATGAA...GGGG..GCCTTTCATCTGGCCGAAAATATGCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgc                                                       <  2:2175851/101‑1 (MQ=37)
                                                 aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg                                                      >  1:85512/1‑101 (MQ=255)
                                                  aGATGAG...GGGG..TGTTTTAATCTGGCAAAAAATACCCGGGGGGGTGAA.GGGCGCA..GCCAGGGGGGGG.CAGCCC.CCCCCTTCGCGGGCGGCCGCGACgcgcgc                                                     <  2:3958692/101‑1 (MQ=255)
                                                  aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACgcgcgc                                                     >  1:5437517/1‑101 (MQ=255)
                                                   gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgcc                                                    <  1:428638/101‑1 (MQ=255)
                                                   gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgcc                                                    >  1:4795376/1‑101 (MQ=255)
                                                      gga...GGGG..GGTTTTTCTTTTTGCGAAAATTGCCCGGGGGGGGAAGGGGCCCC..GCCGGGGGGGGGCCCC.CC.CCCCCTTCGCGGGCGGCCGCGACGCGcgccgc                                                  <  1:1965046/99‑1 (MQ=255)
                                                         gatggtg..ggCTTTCAACTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                 <  2:744295/95‑1 (MQ=25)
                                                             ggg..ggTTTTCATCTTGCCGAAAAAACCCGGGGGGGTGAA.TGGCGCA..GCCAGGGGGGGC.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                             <  1:4228043/101‑1 (MQ=37)
                                                             ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                             <  1:4987746/101‑1 (MQ=255)
                                                             ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                             >  2:1409554/1‑101 (MQ=255)
                                                               g..ggCTTTCATCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGt                                           <  1:2292518/101‑1 (MQ=255)
                                                               g..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGGGCA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCTCGGGCGGCAGGGACGCGGGACGCAACAgga                                           >  2:4375029/1‑100 (MQ=21)
                                                                   gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCACCAGGTCg                                         >  1:1329735/1‑101 (MQ=255)
                                                                    cTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGt                                        <  1:4143724/101‑1 (MQ=255)
                                                                         aTCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGGCGTGcccc                                   >  2:4949254/1‑101 (MQ=37)
                                                                          tCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                  <  2:2237122/101‑1 (MQ=255)
                                                                          tCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                  >  2:563883/1‑101 (MQ=255)
                                                                           cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                 <  2:4269389/101‑1 (MQ=255)
                                                                            ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                >  1:5744331/1‑101 (MQ=255)
                                                                            ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                <  2:3581247/101‑1 (MQ=255)
                                                                             tGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGAGGCGCGCCGCAACAGGTCGTGCCCCgttc                               >  1:1270002/1‑100 (MQ=37)
                                                                             tGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                               >  2:5059159/1‑101 (MQ=255)
                                                                              gCCGAAAATATCCCAGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGt                              >  2:1020012/1‑101 (MQ=255)
                                                                               ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                             >  2:2255458/1‑101 (MQ=255)
                                                                                 gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                           >  2:422137/1‑101 (MQ=255)
                                                                                 gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                           <  1:3739804/101‑1 (MQ=255)
                                                                                 gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                           <  2:1671925/101‑1 (MQ=255)
                                                                                 gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGCGCCCCTTTGTCAt                           >  1:2707275/1‑101 (MQ=255)
                                                                                  aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGTG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                          >  1:5656514/1‑101 (MQ=255)
                                                                                    aaTATCCGGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                        <  1:1160701/101‑1 (MQ=38)
                                                                                    aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                        >  2:4098105/1‑101 (MQ=255)
                                                                                     atatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTg                       <  2:5282009/101‑1 (MQ=255)
                                                                                     atatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGACGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTg                       >  2:341770/1‑101 (MQ=255)
                                                                                        tCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCg                    <  1:2677914/101‑1 (MQ=255)
                                                                                          ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa                  <  2:2575031/101‑1 (MQ=255)
                                                                                           cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                 >  2:2127671/1‑101 (MQ=255)
                                                                                           cGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGG.CAGCCC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                 <  1:1499431/101‑1 (MQ=255)
                                                                                            ggggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                >  2:1095929/1‑101 (MQ=255)
                                                                                             gggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc               >  2:2077836/1‑101 (MQ=255)
                                                                                             gggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc               >  2:5557038/1‑101 (MQ=255)
                                                                                              ggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa              >  1:62254/1‑101 (MQ=255)
                                                                                              ggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa              >  1:559764/1‑101 (MQ=255)
                                                                                               gggTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt             <  2:4478240/101‑1 (MQ=255)
                                                                                               gggTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt             >  2:1107115/1‑101 (MQ=255)
                                                                                                ggTGAA.TGGCCCA..CCCGGGGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg            <  1:5608619/101‑1 (MQ=255)
                                                                                                    aa.TGGCGCA..GCCAGAGGGGGG.CGG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt        >  2:2709452/1‑101 (MQ=255)
                                                                                                        ggCGCA..GCGAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg     <  1:2474238/101‑1 (MQ=255)
                                                                                                        ggCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg     >  2:2675298/1‑101 (MQ=255)
                                                                                                        ggCGCA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg     >  2:1006920/1‑101 (MQ=255)
                                                                                                           gcA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg  >  1:2018738/1‑101 (MQ=255)
                                                                                                           gcA..GCCAGAGGGGGG.CAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg  >  1:5488282/1‑101 (MQ=255)
                                                                                                           |                                                                                                        
GGCTCTG.CCC.CGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAG.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTG  >  NC_014034/2249274‑2249471

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