Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,372 G→C 3.5% 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 3.5% 6.2 0.0 ‑0.3184intergenic (‑168/‑137)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (95/81):  total (99/85)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.57e-01

Only 100 of 205 total aligned reads displayed.

GGGCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAG  >  NC_014034/2249273‑2249468
                                                                                                             |                                                                                                    
gggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGAGAA.TGgcg                                                                                                      >  1:638457/1‑101 (MQ=255)
 ggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGgcgc                                                                                                     <  1:2191454/101‑1 (MQ=255)
 ggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGgcgc                                                                                                     >  1:2461494/1‑101 (MQ=255)
  gCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCa                                                                                                    <  2:262919/101‑1 (MQ=255)
  gCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCa                                                                                                    <  2:2587633/101‑1 (MQ=255)
     ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..Gcc                                                                                               <  2:3403963/97‑1 (MQ=255)
        cccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCaga                                                                                            >  1:649210/2‑101 (MQ=255)
         ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCagag                                                                                           >  1:3318288/1‑101 (MQ=255)
           ccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggg                                                                                         >  1:1183387/1‑101 (MQ=255)
           ccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAggg                                                                                         >  1:4304189/1‑101 (MQ=255)
            cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                                        >  1:679210/1‑101 (MQ=255)
            cgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggg                                                                                        >  1:493782/1‑101 (MQ=255)
             gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCGGAggggg                                                                                       <  2:1448065/101‑1 (MQ=255)
             gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                       >  1:3628503/1‑101 (MQ=255)
              tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                      >  1:2773954/1‑101 (MQ=255)
               cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                     <  2:3544450/101‑1 (MQ=255)
               cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAggggggg                                                                                     >  1:351778/1‑100 (MQ=255)
                ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                    <  2:3733945/101‑1 (MQ=255)
                 cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg                                                                                   <  2:454360/101‑1 (MQ=255)
                 cgtCCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...TGGG..C..GCTTTCATCTTCCCGAAAATATCCCGGGGGGTGAA.TGGCCCA..GCCAGAGGGGGGCAg                                                                                   <  2:3852007/101‑1 (MQ=25)
                  gcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAgggggg                                                                                      >  2:401988/7‑101 (MQ=18)
                   tGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg                                                                                <  2:4216260/101‑1 (MQ=255)
                   caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGc                                                                                     <  2:4058991/96‑1 (MQ=18)
                    gCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAg.cg.c                                                                              >  1:3413164/1‑101 (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)
                      cGGGACG.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccc                                                                            <  1:3475932/101‑1 (MQ=255)
                       gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccc                                                                           <  2:4247913/101‑1 (MQ=255)
                       gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccc                                                                           >  1:4100865/1‑101 (MQ=255)
                       gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.cccc                                                                           >  1:3856632/1‑101 (MQ=255)
                       cggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.cg                                                                                <  1:308289/100‑1 (MQ=18)
                       cggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..G..CCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCAg.cg                                                                                <  1:2545834/100‑1 (MQ=18)
                            t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccccttc                                                                      <  2:741661/101‑5 (MQ=255)
                            t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccccttc                                                                      <  1:1780179/101‑5 (MQ=255)
                            t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccccttc                                                                      <  1:3700401/101‑5 (MQ=255)
                              cccccGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.ccccccttcc                                                                     >  2:1358030/1‑96 (MQ=255)
                              cccccGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAA.TGGGGCA..GCCAGAGGGGGGCAG.CG.ccccccttcg                                                                     >  2:3336467/1‑96 (MQ=255)
                               ccccGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTcgc                                                                     <  2:1977061/101‑1 (MQ=255)
                                  cGGGATATTTGCGGCAAGATGAA...GGGGCAT..GATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCCtgccc                                                                     >  2:3971647/1‑97 (MQ=0)
                                  cGGGATATTTGCGGCAAGATGAA...GGGG..CATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCCtgccc                                                                     >  1:3361751/1‑97 (MQ=0)
                                  cGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCgcggg                                                                  >  1:1804492/1‑101 (MQ=255)
                                   gggATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggc                                                                 >  1:4656833/1‑101 (MQ=255)
                                    ggATATTTGCGGCAAGATGAA...GGGGCAT..GATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCAG.CG.CCCCCtgccccg                                                                   <  2:3054626/101‑7 (MQ=0)
                                    ggATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggcg                                                                >  1:2624461/1‑101 (MQ=255)
                                      atatTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGggcggc                                                              <  2:1624353/101‑1 (MQ=255)
                                       tatTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGcc                                                             >  2:4603816/1‑101 (MQ=255)
                                        atTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCcg                                                            >  1:3942320/1‑101 (MQ=255)
                                        atTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCcg                                                            >  1:4445668/1‑101 (MQ=255)
                                              ggCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgc                                                      <  1:1910438/101‑1 (MQ=255)
                                                cAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.AGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCTTGCGGCCGCGAcgcgc                                                    >  1:779436/1‑101 (MQ=25)
                                                  aGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGCCgcgcgc                                                  >  1:4456455/1‑101 (MQ=255)
                                                   gATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgcc                                                 <  1:2975977/101‑1 (MQ=255)
                                                   gATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgcc                                                 >  2:2469247/1‑101 (MQ=255)
                                                     tGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgccgc                                               >  2:1996427/1‑101 (MQ=255)
                                                       aa...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGAGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCaa                                             >  1:3556136/1‑101 (MQ=255)
                                                            gggg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAG.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                           >  1:4271465/1‑101 (MQ=255)
                                                             ggg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                          >  2:2431284/1‑101 (MQ=255)
                                                             ggg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                          <  2:2208392/101‑1 (MQ=255)
                                                             ggg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                          >  1:124686/1‑101 (MQ=255)
                                                              gg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAgg                                         <  2:4111817/101‑1 (MQ=255)
                                                              gg..g..gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAgg                                         <  1:2115580/101‑1 (MQ=255)
                                                               g..g..gCTTTCATCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGt                                        <  1:3767092/101‑1 (MQ=255)
                                                                      cTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGt                                     >  1:4052644/1‑101 (MQ=255)
                                                                        ttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                                   >  1:1385125/1‑101 (MQ=255)
                                                                          cATCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                 <  1:4206859/101‑1 (MQ=255)
                                                                          cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                 >  1:2900076/1‑101 (MQ=255)
                                                                          cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                 >  2:2637993/1‑101 (MQ=255)
                                                                          cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                                 >  1:4154009/1‑101 (MQ=255)
                                                                          cATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCCCAACAGGTCGTGccc                                 >  2:1122348/1‑101 (MQ=255)
                                                                           aTCTTGCCGAAAATATCCGGGGGGGTGAA.TGGCTCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACATGTCGTGcccc                                <  1:4032399/101‑1 (MQ=37)
                                                                           aTCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGCCGTGcccc                                >  1:2445144/1‑101 (MQ=255)
                                                                             cTTGGCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCATAGGGTGGCAG.CGCCCCCCTGCGCGGGCGGCCGCGACTCTCGCCGCACCAGGTCGTGCCCCtt                              <  2:2329171/101‑1 (MQ=21)
                                                                             cTTGCGGAAAATATCCGGGGGGGTGAA.GGGCGCA..GCCGGGGGGGGGCAGCCC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                              <  1:997103/101‑1 (MQ=255)
                                                                             cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                              >  2:2570841/1‑101 (MQ=255)
                                                                             cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                              >  1:379997/1‑101 (MQ=255)
                                                                                  cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                         <  2:3694959/101‑1 (MQ=255)
                                                                                  cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                         >  1:1066059/1‑101 (MQ=255)
                                                                                   gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                        >  2:2994249/1‑101 (MQ=255)
                                                                                     aattcgccaggggggGGAG.TGCGGCA..GCGAGCGGGGGGCAC.CGCCCCCCTCCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                      <  1:203184/92‑1 (MQ=255)
                                                                                      aaTGTCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                     <  1:4579511/101‑1 (MQ=255)
                                                                                       atatCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGCCACGCGCCGCAACAGGTCGTGCCCCTTGGTCAGTCTg                    >  2:2276934/1‑101 (MQ=37)
                                                                                        tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                   >  2:3475903/1‑101 (MQ=255)
                                                                                        tatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGt                   <  1:523221/101‑1 (MQ=255)
                                                                                            ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa               >  1:927373/1‑101 (MQ=255)
                                                                                             cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc              <  2:3898244/101‑1 (MQ=255)
                                                                                             cGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGAGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc              >  2:1807724/1‑101 (MQ=255)
                                                                                               gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc            >  2:3150915/1‑101 (MQ=255)
                                                                                               gggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc            <  2:1395583/101‑1 (MQ=255)
                                                                                                ggggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa           >  2:3183693/1‑101 (MQ=255)
                                                                                                 gggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt          >  2:862886/1‑101 (MQ=255)
                                                                                                  ggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg         <  2:1147735/101‑1 (MQ=255)
                                                                                                  ggTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg         >  1:1978177/1‑101 (MQ=255)
                                                                                                   gTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt        <  1:2279293/101‑1 (MQ=255)
                                                                                                    tGAA.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc       >  2:1055919/1‑101 (MQ=255)
                                                                                                      aa.TGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt     >  2:2807285/1‑101 (MQ=255)
                                                                                                         tGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa   <  2:649149/101‑1 (MQ=255)
                                                                                                         tGGCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa   <  2:3062261/101‑1 (MQ=255)
                                                                                                         gggCGCA..GCGAGGGGGGGGCAGCCC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa   <  1:157848/100‑1 (MQ=255)
                                                                                                          ggCGCA..GCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg  <  1:326527/101‑1 (MQ=255)
                                                                                                             |                                                                                                    
GGGCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..G..GCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCAG.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAG  >  NC_014034/2249273‑2249468

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