Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 460,816 G→A 2.4% intergenic (‑246/+150) creA ← / ← mdtB creatinase/multidrug resistance protein MdtB

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034460,8160G→A 2.4% 2.8 ‑1.2 ‑0.7168intergenic (‑246/+150)creA/mdtBcreatinase/multidrug resistance protein MdtB
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (112/51):  total (116/52)
Fisher's exact test strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.07e-01

Only 100 of 224 total aligned reads displayed.

CAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGTTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCTTTGGGCCCCGGCGGGTCTTGCTGCAGAT  >  NC_014034/460716‑460916
                                                                                                    |                                                                                                       
cAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTg                                                                                                         >  2:5486171/1‑101 (MQ=255)
 aaGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGa                                                                                                        >  1:2454679/1‑101 (MQ=255)
 aaGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGa                                                                                                        >  1:1716783/1‑101 (MQ=255)
  aGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGTGACCGTGTCCTCAAGCGTGAc                                                                                                       <  2:4798306/101‑1 (MQ=255)
   ggggCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACa                                                                                                      <  2:2887888/101‑1 (MQ=255)
     ggcggcGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACaga                                                                                                    <  1:3415830/101‑1 (MQ=255)
      gcggcgATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACagag                                                                                                   >  1:1045147/1‑101 (MQ=255)
        ggcgATGTGCCGGGTGCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGc                                                                                                 >  1:4882278/1‑101 (MQ=255)
         gcgATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGcc                                                                                                <  1:5720742/101‑1 (MQ=255)
          cgATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCg                                                                                               >  2:4870025/1‑101 (MQ=255)
           gATGTGCCGGGTTGGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGcccg                                                                                              >  1:3850864/1‑99 (MQ=255)
             tgtgCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGgcg                                                                                            <  1:4231636/101‑1 (MQ=255)
             tgtgCCGGGTTCGGGCAGATCAGGAAAAGCCCCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGgcg                                                                                            <  2:2775140/101‑1 (MQ=255)
              gtgCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCTTGACAGAGGCCGgcgc                                                                                           >  1:1138107/1‑101 (MQ=255)
              gtgCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCgagat                                                                                           >  2:5011188/1‑97 (MQ=255)
               tgCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcc                                                                                          <  2:3120715/101‑1 (MQ=255)
                gCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGGgccc                                                                                         >  2:5700444/1‑101 (MQ=255)
                gCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgccc                                                                                         <  2:3956665/101‑1 (MQ=255)
                 ccGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcccg                                                                                        <  1:3020206/101‑1 (MQ=255)
                   gggTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcccggc                                                                                      >  2:511561/1‑101 (MQ=255)
                    ggTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcccggcc                                                                                     >  2:1225309/1‑101 (MQ=255)
                     gTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCTTGTCCTCAAGCGTGACAGAGGCCGGCgcccggtcc                                                                                    <  1:3304456/101‑1 (MQ=255)
                     gTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcccggccc                                                                                    <  1:4712922/101‑1 (MQ=255)
                      ttCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcccggcccg                                                                                   <  1:1311060/101‑1 (MQ=255)
                      ttCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcccggcccg                                                                                   <  2:5310417/101‑1 (MQ=255)
                            cGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCggcgggc                                                                             <  2:1345527/101‑1 (MQ=255)
                             ggATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCggcgggcg                                                                            >  1:1425612/1‑101 (MQ=255)
                             ggATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAACCGTGACAGAGGCCGGGGCCCGGCCCggcgggcg                                                                            >  1:5017305/1‑101 (MQ=255)
                             ggATCAGGAAAAGCACCGCTTTCCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCggcgggcg                                                                            <  2:1881536/101‑1 (MQ=255)
                              gATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGGGCCCGGCCCagcgggcgg                                                                           >  2:3420269/1‑100 (MQ=255)
                              gATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCggcgggggg                                                                           >  1:5497697/1‑98 (MQ=255)
                                tCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGGGCCCGGCCCCGCGGGCgaga                                                                         >  2:388091/1‑101 (MQ=255)
                                 cAGGAAAAGCCCCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGaa                                                                        <  1:4555875/101‑1 (MQ=255)
                                 cAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGaa                                                                        <  2:35636/101‑1 (MQ=255)
                                   ggAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAgc                                                                      <  2:407962/101‑1 (MQ=255)
                                     aaaaGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAgc.gc                                                                   >  2:5321535/1‑101 (MQ=255)
                                      aaaGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGGGCCCGCCCCGGCGGGCGAGAAAc.gcc                                                                  >  2:322296/1‑101 (MQ=25)
                                      aaaGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGc.gcc                                                                  <  1:5080628/101‑1 (MQ=255)
                                       aaGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGGGGCCCGCCCGGGGGGGGGGGAGc.gccc                                                                 >  2:814446/1‑101 (MQ=37)
                                        aGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGc.gcccg                                                                <  1:4499235/101‑1 (MQ=255)
                                         gCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGGGCCCCGCCCGGCGGGCGGGAAGc.gccccc                                                               >  1:3745002/1‑99 (MQ=37)
                                         gCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGGGAGAAGc.gcccgc                                                               >  2:3317239/1‑101 (MQ=255)
                                            ccGCTTTCCCCCCCCCGTCGGCGAAACGGCTTTCCCGCGACCGTGTCCTAAAGCGGGACAGAGGCGGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCAt                                                            <  2:4378084/101‑1 (MQ=21)
                                             cgcttttCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCCGCCCCGCGGGGGGGAAGc.gcccggcccg                                                           >  1:1987000/4‑96 (MQ=25)
                                              gctttcCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATgg                                                          <  2:1091323/99‑1 (MQ=255)
                                              gctttcCCCCCCCCGTCGGCGCACCGGCTGTTCCGCGCCCGTGCCCTCAAGCGTGACAGAGGCCGGCCCCCGCCCCGGCGGGCGAGAAGC.GCCCGCCATgg                                                          <  2:676586/99‑1 (MQ=21)
                                                 tttCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGACCCGCCCGGCGGGCGAAAAAC.GCCCGCCAAcggcg                                                       >  2:3615884/1‑101 (MQ=37)
                                                 ttcccccGCCCGCCGGCGCAACGGCTGTTCCGCGCCCGTCTCCTCAAGCGTGACAGAGGCCGGCCCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATgggcg                                                       <  2:347254/98‑1 (MQ=21)
                                                    ccccGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCCGCGGGCGAGAACC.GCCCCCCCTgggtgggg                                                    >  2:4245202/1‑101 (MQ=37)
                                                     cccgcccgTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGCCGGGGGGAAAGC.GCCCGCCACgggcggggg                                                   >  2:920494/1‑100 (MQ=37)
                                                     cccccccgGGGGCGAAACGGCTGTTCCGCGACCGTGTCCTCAGGCGTGACAGAGGCCGGCGCCCGGCCGGGCGGGCGAGAAGC.GCCCGCCATgggcggggc                                                   <  1:282867/101‑1 (MQ=18)
                                                      ccgccggCCGGCGCACCGGCTGTTCCGCGCCCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCGGGCGGGCGAGAAGC.GCCCGCCATgggcggggcg                                                  <  2:3843313/101‑1 (MQ=21)
                                                      ccgcccgTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATgggcggggcg                                                  <  1:1789691/101‑1 (MQ=255)
                                                        gcccgTCGGCGCAACGGCTGTTCCGCGCCCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCGGGCGGGGGAGAAGC.GCCCGCCATGGgcggggcggg                                                <  2:2928565/101‑1 (MQ=255)
                                                        gcccgTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGGGAAGC.GCCCGCCCCGGgcgggggggg                                                >  1:2747459/1‑101 (MQ=25)
                                                         cccgTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCCTGGgcggggggggg                                               >  1:739955/1‑100 (MQ=255)
                                                         cccgGCGGCGCAAGGGCTGTCCCGCGCCCGTGCCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGgggcgggc                                               <  1:1592706/101‑1 (MQ=18)
                                                            gTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGGGAAAC.GCCCGCCATGGGCGGGGCGGgcgcg                                            >  2:4477927/1‑100 (MQ=255)
                                                            gTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCCGGCCGGGGGGCGAGAAGC.CCCCGCCAGCGGCGGGGGGGGCGC.t                                           >  1:4303982/1‑101 (MQ=25)
                                                             tagtaccAACGGCTGCTCCGCCTCCGTGTCGTTAAGGTTGAAAAAGGCGGTCCCCCGGCCCGGCGGGCAAGGCTC.TCCCGCCATGGGCGGGGCGGGCGC.Tg                                          <  1:61489/95‑1 (MQ=9)
                                                               ggCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCCTGGGCGGGGCGGGCGC.GGG.c                                       >  2:5460484/1‑101 (MQ=255)
                                                               ggCGCAACGGCTGTCCCGCGACGGTGTCCTCAAGCGTGACAGGGGCCGGCGCCCGGCCGGGGGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.c                                       <  2:179860/101‑1 (MQ=21)
                                                                gcgcAACGGCTGTTCCTCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGGGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.cc                                      <  2:5508033/101‑1 (MQ=255)
                                                                 cgcAAGGGCTTTCCCGCGCCCGTGTCCTCAGGCGTGACAGAGGCCGGCCCCGGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCg                                     <  2:1547797/101‑1 (MQ=17)
                                                                 cgcAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGGGAAGC.GCCCGCCCTGGGCGGGGGGGGCGG.CGG.CCg                                     >  1:4413697/1‑101 (MQ=21)
                                                                    aaCGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGGGGGAAGC.GCCCGCCACGGGGGGGGCGGGGGC.GGG.CCgggg                                  >  1:1170041/1‑101 (MQ=21)
                                                                    aaCGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAAC.GCCCCCCCTGGGCGGGGGGGGGGC.CGG.GCcggg                                  >  1:4981598/1‑101 (MQ=21)
                                                                     aCGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGGGAAGG.GCCCGCCATGGGCGGGGGGGGGGG.CGG.CCggggg                                 >  1:2921982/1‑101 (MQ=21)
                                                                        gCTGTTCCGCCCCCGTGTCCTCAAGGGGGACAGGGGCCGGCCCCCGGCCCGGGGGGCGAGAAGC.CCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCt                              <  1:5379021/101‑1 (MQ=17)
                                                                         cTGTTCCGCGACGGTGCCCTAAGGCGTGACGGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCtt                             <  2:1499846/101‑1 (MQ=17)
                                                                         cTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.CGGCCCGGGGGCCt                              >  1:451443/1‑101 (MQ=255)
                                                                         cTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCCCCATGGGCGGGGCGGGGGC.CGGCCCGGGGGCCt                              >  1:4367314/1‑101 (MQ=18)
                                                                         cTGTCCCGGGCCCGTGTCCTCAGGCGTGCCAGGGGCCGGCGCCCGGCCCGGGGGGGGAGAGGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCtt                             <  1:5671027/101‑1 (MQ=17)
                                                                          tgctCCGCGGCCGTGTCCACAGGCGCGACGGAGCCCGGCGCCCGGCCCGGCGGGGGAGAGCC.CCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCttt                            <  2:111973/98‑1 (MQ=14)
                                                                           gTTCCGCGCCCGTGCCCTCAGGCGGGACAGGGGCCGGCGCCCGGCCGGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCTTTg                           <  2:2906626/101‑1 (MQ=17)
                                                                           gTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.GCGGGGGCCttta                           >  1:5319103/1‑100 (MQ=255)
                                                                           gTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.CGG.GCGGGGGGCCTTg                           >  2:2164378/1‑101 (MQ=18)
                                                                             tCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.CGG.CCGGGGGGCCTTggg                         >  1:2079657/1‑101 (MQ=255)
                                                                             gccaagACCGTTTCCCCAGGCGCAGCAGAGGCCGGCGCCCGCCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGGGGGCGCGGGG.CCCGGGGccctggg                          >  1:2144551/6‑96 (MQ=17)
                                                                             cccGCGGCCGTGCCCCCAAGCGTGACAGAGGCCGGCCCCCGGCCGGGCGGGGGAGAAGC.CCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCTTTggg                         <  1:4850569/100‑1 (MQ=14)
                                                                                gcgACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCCTGGGCGGGGCGGGGGC.GGG.GCGGGGGGCCTTGGGccc                      >  2:4816146/1‑101 (MQ=17)
                                                                                 cgACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGGGC.GGG.CCGGGGGCCTTTGGGgccc                     >  2:1671378/1‑101 (MQ=255)
                                                                                  gACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGGGGGCGC.CGG.CCGGGGGCCTTTTGGCGCcg                    >  2:3336579/1‑101 (MQ=17)
                                                                                  gACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCAAGGGCGGGGCGGGCGC.CGG.GCGGGGGCCCTTTGGCCCcg                    >  2:513611/1‑101 (MQ=17)
                                                                                  ccgatcttCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCTTTGGGCCCcg                    <  2:2550721/94‑1 (MQ=255)
                                                                                   aCCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGGCTTTGGGgccccg                   >  2:3790960/1‑99 (MQ=255)
                                                                                    cggtgtCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCTTTGGGCCCcggc                  <  2:4612151/99‑1 (MQ=255)
                                                                                     cGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCCGCGGGCGAGAAGCGGCCCCCCATGGGGGGGCGGGGGGC.CGG.CCGGGGGGCCTTGGGcccccgg                  >  1:3946708/1‑98 (MQ=17)
                                                                                         tCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGGCTTTGGGGCCcgagatcgg             >  1:2550863/1‑94 (MQ=17)
                                                                                          ccccacgCGTGACGGGGGCCGGCCCCGGGCCCGGCGGGAGAGAGGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCTTTGGGCCCCGGCGGGTCt            <  1:4318545/95‑1 (MQ=17)
                                                                                           cTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCCTTGGGCCCCGGCGGGGCtt           >  2:5410086/1‑101 (MQ=255)
                                                                                            tCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGGGGGGGCGGGG.CCGGGGGCCTTTGGGCCCCGGGGGGGCtt           >  1:4034686/1‑101 (MQ=17)
                                                                                            ccaggCGGGACGGAGGCGGGCCCCGGGCCCGGCGGGGGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCTTTGGGCCCCGGCGGGTCTtg          <  2:3820266/97‑1 (MQ=17)
                                                                                             cAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGGGC.CTG.CCGGGGGCCGTTGGGACCCGGCGGGGCGtgc         >  1:3153878/1‑101 (MQ=17)
                                                                                             cAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGGCCTTGGGCCCCGGCGGGGCTtgc         >  2:1968374/1‑101 (MQ=255)
                                                                                              aaGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCTTTGGGGCCCGGCGGGGCTtggt        >  2:4463290/1‑99 (MQ=255)
                                                                                                ggtggACAGGGCCCGGCGCCGGGCCCGGCGGGCGAGAGCC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCTTTGGGCCCCGGCGGGTCTtgctgc      <  2:302074/97‑1 (MQ=17)
                                                                                                 cGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCCTTGGGCCCCGGCGGGGCTtgctgcc     >  1:5345256/1‑100 (MQ=255)
                                                                                                   ggACAGAGGCCGGCGCCCGGCCGGGGGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCTTTGGGCCCCGGCGGGTCTTGCTGCAGa   <  1:920570/100‑1 (MQ=255)
                                                                                                    gACAGAGGCCGGCCCCCGGCCCGGGGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCTTTGGGCCCCGGCGGGTCTTGCTGCAGAt  <  1:3238114/101‑1 (MQ=255)
                                                                                                    |                                                                                                       
CAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGTTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGC.GCCCGCCATGGGCGGGGCGGGCGC.TGG.CCGGGGGCCTTTGGGCCCCGGCGGGTCTTGCTGCAGAT  >  NC_014034/460716‑460916

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