Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,379 G→A 3.4% 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 3.4% 6.5 ‑0.4 ‑0.0175intergenic (‑175/‑130)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (4/2):  ref base (94/75):  total (98/77)
Fisher's exact test strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.45e-01

Only 100 of 203 total aligned reads displayed.

TGCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAT  >  NC_014034/2249279‑2249478
                                                                                                             |                                                                                                        
gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                                          >  2:3371068/3‑101 (MQ=255)
   cccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgg                                                                                                       <  2:2348383/101‑1 (MQ=255)
    ccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAggg                                                                                                      >  1:4304189/1‑101 (MQ=255)
     cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                                                     >  1:493782/1‑101 (MQ=255)
      gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCGGAggggg                                                                                                    <  2:1448065/101‑1 (MQ=255)
       tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                                   >  1:3661675/1‑101 (MQ=255)
        cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                                  <  1:301943/101‑1 (MQ=255)
        cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                                  <  2:3544450/101‑1 (MQ=255)
        cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                                  <  1:3756107/101‑1 (MQ=255)
        cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                                  <  1:1201492/101‑1 (MQ=255)
        cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                                  >  2:4579212/1‑101 (MQ=255)
         ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                                 <  2:3733945/101‑1 (MQ=255)
         ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                                 >  1:4291518/1‑101 (MQ=255)
         ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                                 <  2:3234555/101‑1 (MQ=255)
          cgtCCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...TGGG..CGCTTTCATCTTCCCGAAAATATCCCGGGGGGTGAA.TGGCCCA..GCCAGAGGGGGGCAg                                                                                                <  2:3852007/101‑1 (MQ=25)
          cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                                                    <  2:2140257/94‑1 (MQ=18)
           gtGCCGGGACT.CCCCTGCGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                                                                              <  1:3869997/101‑1 (MQ=255)
           gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                                                                              >  1:3345788/1‑101 (MQ=255)
           gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                                                                              >  1:4170793/1‑101 (MQ=255)
           gcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggggg                                                                                                   <  1:2974769/95‑1 (MQ=18)
           gcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggggg                                                                                                   >  2:401988/7‑101 (MQ=18)
            tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c.g                                                                                            <  2:4216260/101‑1 (MQ=255)
            caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                                                                  <  2:4058991/96‑1 (MQ=18)
             gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c.g.c                                                                                          >  1:3413164/1‑101 (MQ=255)
              ccGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cc                                                                                         >  2:3461441/1‑101 (MQ=255)
               cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.ccc                                                                                        >  2:3060633/1‑101 (MQ=255)
               cGGGACG.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.ccc                                                                                        <  1:3475932/101‑1 (MQ=255)
                gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccc                                                                                       >  2:2321797/1‑101 (MQ=255)
                gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccc                                                                                       <  2:4247913/101‑1 (MQ=255)
                cggccT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccc                                                                                       <  2:3154963/97‑1 (MQ=255)
                cggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.c.g                                                                                            <  1:2545834/100‑1 (MQ=18)
                   aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.cccccct                                                                                    >  1:3062505/1‑99 (MQ=255)
                   aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGACGGGGGCGG.C.G.cccccct                                                                                    >  2:3064665/1‑99 (MQ=255)
                     t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.ccccccttc                                                                                  >  2:157834/1‑97 (MQ=255)
                     t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.ccccccttc                                                                                  >  2:538975/1‑97 (MQ=255)
                     t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.ccccccttc                                                                                  <  2:741661/101‑5 (MQ=255)
                       cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.ccccccttcc                                                                                 >  2:1358030/1‑96 (MQ=255)
                         cccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.C.G.CCCCC..tgc                                                                                 >  2:1905000/1‑99 (MQ=0)
                         cccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTcgcg                                                                              <  1:1250189/101‑1 (MQ=255)
                           cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.C.G.CCCCC..tgccc                                                                               >  2:3971647/1‑97 (MQ=0)
                            gggATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.C.G.CCCCC..tgcccc                                                                              >  2:3008472/1‑96 (MQ=2)
                            gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCgcgggc                                                                           >  1:1799212/1‑101 (MQ=255)
                            gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCgcgggc                                                                           >  1:3520456/1‑101 (MQ=255)
                             ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCgcgggcg                                                                          >  1:2624461/1‑101 (MQ=255)
                                 atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCcg                                                                      >  1:4445668/1‑101 (MQ=255)
                                 atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCcg                                                                      <  2:227603/101‑1 (MQ=255)
                                  tttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCcgc                                                                     <  1:3508841/101‑1 (MQ=255)
                                       ggCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGAcgc                                                                <  1:1910438/101‑1 (MQ=255)
                                           aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACgcgcgc                                                            >  1:3308320/1‑101 (MQ=255)
                                           aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACgcgcgc                                                            <  1:6732/101‑1 (MQ=255)
                                            gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGcgcc                                                           >  2:2469247/1‑101 (MQ=255)
                                            gATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGcgcc                                                           >  2:1985429/1‑101 (MQ=255)
                                              tGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGTGAA.AGGCGCA..GCCAGAAGGGGGGAG.CGG.CCCCC..CTCGCCGGCGGCCGCGGCGCGCGCCGCa                                                        >  2:4505292/1‑101 (MQ=25)
                                               gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCa                                                        <  1:4265493/101‑1 (MQ=255)
                                                 a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGAAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                      <  1:4115703/101‑1 (MQ=255)
                                                     gggg..ggCTTTCATCTTGCCGAAAATATCCGGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                     <  1:3371318/101‑1 (MQ=37)
                                                     gggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                     <  1:1446647/101‑1 (MQ=255)
                                                     gggg..ggCTTTCATCTTGCCGAAAATAGCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                     <  1:4622264/101‑1 (MQ=255)
                                                      ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                    >  1:200484/1‑101 (MQ=255)
                                                      ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                    >  1:3794724/1‑101 (MQ=255)
                                                       gg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAgg                                                   <  2:4111817/101‑1 (MQ=255)
                                                              tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                              >  1:2049931/1‑101 (MQ=255)
                                                              tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                              <  1:321237/101‑1 (MQ=255)
                                                               ttCATCATGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                                             <  1:768392/101‑1 (MQ=255)
                                                                 cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                           >  2:2637993/1‑101 (MQ=255)
                                                                 cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                           >  1:2855236/1‑101 (MQ=255)
                                                                 cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCCCAACAGGTCGTGccc                                           >  2:1122348/1‑101 (MQ=255)
                                                                  aTCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCTCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACATGTCGTGcccc                                          <  1:4032399/101‑1 (MQ=37)
                                                                  aTCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGCCGTGcccc                                          >  1:2445144/1‑101 (MQ=255)
                                                                    cTTGCGGAAAATATCCGGGGGGGTGAA.GGGCGCA..GCCGGGGGGGGGCAGCC.C.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                        <  1:997103/101‑1 (MQ=255)
                                                                    cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                        >  2:2626868/1‑101 (MQ=255)
                                                                     ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                       <  1:3628035/101‑1 (MQ=255)
                                                                         cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                                   <  2:3694959/101‑1 (MQ=255)
                                                                          gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                  >  1:682579/1‑101 (MQ=255)
                                                                           aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                 >  1:2161595/1‑101 (MQ=255)
                                                                             aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                               >  1:113275/1‑101 (MQ=255)
                                                                             aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                               >  1:445804/1‑101 (MQ=255)
                                                                               tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                             >  2:2572905/1‑101 (MQ=255)
                                                                               tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                             <  1:523221/101‑1 (MQ=255)
                                                                               tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGAGCCCCTTTGTCATTCTgc                             >  1:3826008/1‑100 (MQ=255)
                                                                                   ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa                         >  1:927373/1‑101 (MQ=255)
                                                                                    cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                        <  2:3898244/101‑1 (MQ=255)
                                                                                    cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCCCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTTGCAc                        >  1:1608715/1‑101 (MQ=38)
                                                                                    cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGAGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                        >  2:1807724/1‑101 (MQ=255)
                                                                                      gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                      <  2:1395583/101‑1 (MQ=255)
                                                                                      gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                      >  2:349856/1‑101 (MQ=255)
                                                                                      gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                      <  1:644304/101‑1 (MQ=255)
                                                                                        gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                    <  2:4291230/101‑1 (MQ=255)
                                                                                         ggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                   <  2:1147735/101‑1 (MQ=255)
                                                                                          gTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                  <  1:2279293/101‑1 (MQ=255)
                                                                                           tGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc                 >  2:1055919/1‑101 (MQ=255)
                                                                                            gAA.TGGCGCA..GCCAGGGGGGGGCAGCC.C.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa                <  1:4066306/101‑1 (MQ=255)
                                                                                             aa.TGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt               >  2:2807285/1‑101 (MQ=255)
                                                                                                tGGCGCA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa             <  2:649149/101‑1 (MQ=255)
                                                                                                     cA..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg        >  1:1658780/1‑101 (MQ=255)
                                                                                                      a..GCCAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGc       >  2:2905234/1‑101 (MQ=255)
                                                                                                          ccggagGGGGGCAGCC.C.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCg     <  1:1049738/98‑1 (MQ=255)
                                                                                                          ccAGAGGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCg     >  2:612521/1‑101 (MQ=255)
                                                                                                            agagGGGGGCAG.C.GCCCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa   >  2:897599/1‑101 (MQ=255)
                                                                                                             gagGGGGGCAG.C.G.CCCCCCTTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa   <  2:1649322/101‑1 (MQ=255)
                                                                                                             |                                                                                                        
TGCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.C.G.CCCCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAT  >  NC_014034/2249279‑2249478

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