Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,389 +C 3.1% intergenic (‑185/‑120) 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,3891.→C 3.1% 7.1 ‑4.2 0.0161intergenic (‑185/‑120)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (90/66):  total (90/71)
Fisher's exact test strand distribution p-value = 1.54e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 193 total aligned reads displayed.

CGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTCC  >  NC_014034/2249289‑2249486
                                                                                                                |                                                                                                      
cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                         <  1:2241603/101‑1 (MQ=255)
cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                                         <  2:454360/101‑1 (MQ=255)
 gtGCCGGGACT.CCCCTGCGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.c                                                                                                       <  1:3869997/101‑1 (MQ=255)
  tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                                                                      <  2:4216260/101‑1 (MQ=255)
   gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg..c                                                                                                   >  2:594752/1‑101 (MQ=255)
    ccGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc                                                                                                  <  2:4312041/101‑1 (MQ=255)
     ccgggctcccccaGGGATATTTGAGGCAAGATGAG...GGGG..G..GCTTTCATCTTGACGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc                                                                                                  <  1:739917/93‑1 (MQ=25)
     cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc.c                                                                                                >  2:3060633/1‑101 (MQ=255)
     cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc.c                                                                                                <  2:4022810/101‑1 (MQ=255)
      gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc.cc                                                                                               >  1:4100865/1‑101 (MQ=255)
      gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc.cc                                                                                               >  2:2321797/1‑101 (MQ=255)
      cggccT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc.cc                                                                                               <  2:3154963/97‑1 (MQ=255)
      cggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.cg                                                                                                      <  1:2545834/100‑1 (MQ=18)
       ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.cg..c                                                                                                   <  1:2451327/101‑1 (MQ=17)
        gactccccccGGGATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAG.CC..cc                                                                                                  <  1:3450227/96‑1 (MQ=17)
         aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc.cccct                                                                                            >  1:3062505/1‑99 (MQ=255)
         aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGACGGGGGCGG.CG..cc.cccct                                                                                            >  2:3064665/1‑99 (MQ=255)
           t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc.ccccttc                                                                                          <  1:3700401/101‑5 (MQ=255)
           t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc.ccccttc                                                                                          >  2:157834/1‑97 (MQ=255)
           t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc.ccccttc                                                                                          <  2:741661/101‑5 (MQ=255)
           t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc.ccccttc                                                                                          <  1:1780179/101‑5 (MQ=255)
             cccccGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..cc.ccccttcg                                                                                         <  2:1288833/101‑6 (MQ=255)
             cccccGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGGGCA..GCCAGAGGGGGGCAG.CG..cc.ccccttcg                                                                                         >  2:3336467/1‑96 (MQ=255)
               cccGGGATATTTGCGGCAAGATGAA...GGGGCAT..GATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG..CC.CCC..tgc                                                                                         >  2:1905000/1‑99 (MQ=0)
                 cGGGATATTTGCGGCAAGATGAA...GGGGCAT..GATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG..CC.CCC..tgccc                                                                                       >  2:3971647/1‑97 (MQ=0)
                 cGGGATATTTGCGGCAAGATGAA...GGGG..CATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG..CC.CCC..tgccc                                                                                       >  1:3361751/1‑97 (MQ=0)
                 cGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCgcggg                                                                                    >  1:1804492/1‑101 (MQ=255)
                  gggATATTTGCGGCAAGATGAA...GGGGCAT..GATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG..CC.CCC..tgcccc                                                                                      >  2:3008472/1‑96 (MQ=2)
                  gggATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCgcgggc                                                                                   >  1:3520456/1‑101 (MQ=255)
                   ggATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCgcgggcg                                                                                  >  1:2624461/1‑101 (MQ=255)
                     atatTTGCGGCAAAATGAA...GGGGCAT..GATTTCATCTTGCCCAAAATATCCC.GGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG..CC.CCC..tgccccgcc                                                                                   <  2:1006226/101‑9 (MQ=9)
                       atTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCcg                                                                              <  2:227603/101‑1 (MQ=255)
                         ttGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCcgcg                                                                            >  2:3916653/1‑101 (MQ=255)
                          ttaGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGa                                                                           <  2:2773756/99‑1 (MQ=255)
                          tGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGa                                                                           <  2:3821064/101‑1 (MQ=255)
                               cAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.AGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCTTGCGGCCGCGAcgcgc                                                                      >  1:779436/1‑101 (MQ=25)
                                aaGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGAcgcgcc                                                                     >  1:516387/1‑100 (MQ=255)
                                 aGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGCCgcgcgc                                                                    >  1:4456455/1‑101 (MQ=255)
                                 aGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACgcgcgc                                                                    >  1:3308320/1‑101 (MQ=255)
                                  gATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGcgcc                                                                   >  2:1985429/1‑101 (MQ=255)
                                  gATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGcgcc                                                                   >  2:2469247/1‑101 (MQ=255)
                                  gATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGcgcc                                                                   <  1:2975977/101‑1 (MQ=255)
                                      aa...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGAGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                               >  1:3556136/1‑101 (MQ=255)
                                           gggg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAG.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                             >  1:4271465/1‑101 (MQ=255)
                                           gggg..g..gCTTTCATCTTGCCGAAAATAGCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                                             <  1:4622264/101‑1 (MQ=255)
                                            ggg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                            >  1:124686/1‑101 (MQ=255)
                                            ggg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                            >  1:3794724/1‑101 (MQ=255)
                                            ggg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                            >  2:2431284/1‑101 (MQ=255)
                                             gg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAgg                                                           <  2:4111817/101‑1 (MQ=255)
                                             gg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAgg                                                           <  1:2115580/101‑1 (MQ=255)
                                                     cTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGt                                                       >  1:4052644/1‑101 (MQ=255)
                                                      tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                                      <  1:321237/101‑1 (MQ=255)
                                                      tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                                      <  1:4387094/101‑1 (MQ=255)
                                                      tttCATCTTGCAGAAAATACCCGGGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCCGCCC..CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                                      <  1:943922/101‑1 (MQ=255)
                                                       ttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                                                     >  1:1385125/1‑101 (MQ=255)
                                                         cATCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                                   <  1:4206859/101‑1 (MQ=255)
                                                         cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                                   >  1:2855236/1‑101 (MQ=255)
                                                         cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                                   >  1:4154009/1‑101 (MQ=255)
                                                         cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCCCAACAGGTCGTGccc                                                   >  2:1122348/1‑101 (MQ=255)
                                                          aTCTTGCCGAAAATATCCCGGGGGGGGAA.TGGGGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGGCGGGcccc                                                  >  2:1368781/1‑101 (MQ=25)
                                                            cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                                >  2:2570841/1‑101 (MQ=255)
                                                             ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                               <  1:3628035/101‑1 (MQ=255)
                                                             ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                               >  2:2635338/1‑101 (MQ=255)
                                                                 cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                                           <  2:3694959/101‑1 (MQ=255)
                                                                 cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                                           >  1:1066059/1‑101 (MQ=255)
                                                                  gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                          <  1:2321959/101‑1 (MQ=255)
                                                                  gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                          >  1:682579/1‑101 (MQ=255)
                                                                   aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                         >  1:3955649/1‑101 (MQ=255)
                                                                    aattcgccaggggggGGAG.TGCGGCA..GCGAGCGGGGGGCAC.CGC.CC.CCC..TCCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                        <  1:203184/92‑1 (MQ=255)
                                                                     aagcccgtggggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                                       <  1:594813/93‑1 (MQ=38)
                                                                     aaTGTCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                                       <  1:4579511/101‑1 (MQ=255)
                                                                       tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                                     >  2:3475903/1‑101 (MQ=255)
                                                                            cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                                >  2:388471/1‑101 (MQ=255)
                                                                            cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGAGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                                >  2:1807724/1‑101 (MQ=255)
                                                                              gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                              <  1:644304/101‑1 (MQ=255)
                                                                                 ggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                           >  1:1978177/1‑101 (MQ=255)
                                                                                 ggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg                           <  2:1147735/101‑1 (MQ=255)
                                                                                     aa.TGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt                       >  2:2807285/1‑101 (MQ=255)
                                                                                        tGGCGCA..GCCAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa                     <  2:649149/101‑1 (MQ=255)
                                                                                        gggCGCA..GCGAGGGGGGGGCAGCCC..CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa                     <  1:157848/100‑1 (MQ=255)
                                                                                                 gCCAGGGGGGGGCCC.CGCCCC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGc               <  1:1658215/101‑1 (MQ=255)
                                                                                                  ccggagGGGGGCAGCCC..CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCg             <  1:1049738/98‑1 (MQ=255)
                                                                                                  ccAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCg             >  2:612521/1‑101 (MQ=255)
                                                                                                   cAGAGGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGc            >  2:2580717/1‑101 (MQ=255)
                                                                                                    agagGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa           >  2:897599/1‑101 (MQ=255)
                                                                                                    agagGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa           >  2:4287808/1‑101 (MQ=255)
                                                                                                    agagGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa           >  2:3520522/1‑101 (MQ=255)
                                                                                                     gagGGGGGCAGCCC..CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt          <  1:3375958/101‑1 (MQ=255)
                                                                                                     gagGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt          <  1:349240/101‑1 (MQ=255)
                                                                                                     gagGGGGGCAG.CG..CC.CCCCTTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa           <  2:1649322/101‑1 (MQ=255)
                                                                                                      gggggggCGG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATc         <  2:1586010/100‑1 (MQ=255)
                                                                                                      agGGGGGCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATc         <  1:3819061/101‑1 (MQ=255)
                                                                                                        gggggCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGa       >  1:526913/1‑101 (MQ=255)
                                                                                                         ggggCAG.CG..CCACCC..TTCGCGGGCGGCCGCGACGTGCGCGGCAACAGGTCGTGCCCCTTTGACATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGAt      >  1:3851508/1‑101 (MQ=255)
                                                                                                         gcggCAG.CGC.CC.CCC..TGCGCGGGCGGCCCCGCCGCGCGCCGCACCAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGAt      <  2:516339/99‑1 (MQ=255)
                                                                                                          gggCAG.CCC.CC.CCC..TCCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATg     <  1:1736739/101‑1 (MQ=255)
                                                                                                           ggCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGt    <  2:4234862/101‑1 (MQ=255)
                                                                                                           ggCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGt    <  2:2095963/101‑1 (MQ=255)
                                                                                                           ggCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGt    <  2:2965618/101‑1 (MQ=255)
                                                                                                            gCAG.CGC.CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTc   <  2:3412910/101‑1 (MQ=255)
                                                                                                                |                                                                                                      
CGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG..CC.CCC..TTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCATCGATGTCC  >  NC_014034/2249289‑2249486

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