Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 460,808 T→C 2.2% intergenic (‑238/+158) creA ← / ← mdtB creatinase/multidrug resistance protein MdtB

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034460,8080T→C 2.2% 2.1 ‑1.2 ‑1.3179intergenic (‑238/+158)creA/mdtBcreatinase/multidrug resistance protein MdtB
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (111/64):  total (115/64)
Fisher's exact test strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.02e-02

Only 100 of 230 total aligned reads displayed.

CCGCTGGCAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGTTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.GGGGGCCTTTGGGCCCCGGCGGGTCTTG  >  NC_014034/460709‑460908
                                                                                                   |                                                                                                      
ccGCTGGCAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCTCGACCGTGTCCTc                                                                                                       >  1:5508325/1‑101 (MQ=255)
ccGCTGGCAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTc                                                                                                       >  2:4712641/1‑101 (MQ=255)
 cGCTGGCAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCa                                                                                                      <  2:4825867/101‑1 (MQ=255)
 cGCTGGCAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCa                                                                                                      >  2:1789592/1‑101 (MQ=255)
 cGCTGGCAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCa                                                                                                      <  1:1868758/101‑1 (MQ=255)
 cGCTGGCAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCa                                                                                                      <  1:714072/101‑1 (MQ=255)
  gCTGGCAAGGGGCGGCTATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGTCGCAACGGCTGTTCCGCGACCGTGTCCTCaa                                                                                                     <  2:4219991/101‑1 (MQ=255)
  gCTGGCAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCaa                                                                                                     <  1:3717732/101‑1 (MQ=255)
    tGGCAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGc                                                                                                   <  1:1152571/101‑1 (MQ=255)
    tGGCAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGc                                                                                                   <  1:5740928/101‑1 (MQ=255)
         aGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGAc                                                                                              >  2:1108690/1‑101 (MQ=255)
          ggggCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACa                                                                                             <  2:2887888/101‑1 (MQ=255)
          ggggCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACa                                                                                             <  2:1209841/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                                                                                       <  2:4414451/101‑1 (MQ=255)
                gcgATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAAAGGcc                                                                                       <  1:1566838/101‑1 (MQ=255)
                 cgATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCg                                                                                      <  1:4343931/101‑1 (MQ=255)
                 cgATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCg                                                                                      >  2:4870025/1‑101 (MQ=255)
                 atctgtgCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCg                                                                                      <  1:5011469/98‑1 (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)
                      tgCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcc                                                                                 <  2:3120715/101‑1 (MQ=255)
                       gCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGGgccc                                                                                >  2:5700444/1‑101 (MQ=255)
                       gCCGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgccc                                                                                <  2:3956665/101‑1 (MQ=255)
                       ccgatcttCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgccc                                                                                <  1:5697995/95‑1 (MQ=255)
                        ccGGGTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcccg                                                                               <  2:5669971/101‑1 (MQ=255)
                          gggTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcccggc                                                                             <  1:4500418/101‑1 (MQ=255)
                          gggTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcccggc                                                                             >  2:511561/1‑101 (MQ=255)
                            gTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCTTGTCCTCAAGCGTGACAGAGGCCGGCgcccggtcc                                                                           <  1:3304456/101‑1 (MQ=255)
                            gTTCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcccggccc                                                                           <  1:4631117/101‑1 (MQ=255)
                             ttCGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCgcccggcccc                                                                          >  1:4479477/1‑100 (MQ=255)
                               cGGGCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCccggcccggc                                                                        >  1:122377/1‑101 (MQ=255)
                                gggCGGATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGGGCCCGGCCcggcg                                                                       >  1:2307782/1‑101 (MQ=255)
                                 gcgggATCAGGAAAAGCACCGCTTTTGCCCGCCCGTCGGCGCACCGGCTGTTCGGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCGGCCCcggcgg                                                                      <  2:4802009/98‑1 (MQ=25)
                                    ggATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCCGCCCggcgggcg                                                                   >  1:5037958/1‑101 (MQ=255)
                                     gATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGGGCCCGGCCCagcgggcgg                                                                  >  2:3420269/1‑100 (MQ=255)
                                     gATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCga                                                                  <  2:698871/101‑1 (MQ=255)
                                     gATCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGATGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCga                                                                  <  1:4052397/101‑1 (MQ=255)
                                       tCAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGGGCCCGGCCCCGCGGGCgaga                                                                >  2:388091/1‑101 (MQ=255)
                                        cAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGaa                                                               <  2:35636/101‑1 (MQ=255)
                                        cAGGAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCCGCCCGGCGGGCgagga                                                               >  1:3204373/1‑99 (MQ=255)
                                          ggAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAgc                                                             <  2:1019216/101‑1 (MQ=255)
                                          ggAAAAGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAgc                                                             <  2:407962/101‑1 (MQ=255)
                                           gAAAAGCCCCGGTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAgcg                                                            <  2:429709/101‑1 (MQ=255)
                                            aaaaGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCggcgggcggggagggc                                                           >  2:1880005/1‑93 (MQ=255)
                                             aaaGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGGGCCCGCCCCGGCGGGCGAGAAAcgcc                                                          >  2:322296/1‑101 (MQ=25)
                                             aaaGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGcgcc                                                          <  1:5080628/101‑1 (MQ=255)
                                              aaGCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGcgccc                                                         <  2:5341858/101‑1 (MQ=255)
                                               aGCCCCGCTTTTCCCCGCCGGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGcgcccg                                                        <  2:1346891/101‑1 (MQ=255)
                                                gCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGGGCCCCGCCCGGCGGGCGGGAAGcgccccc                                                       >  1:3745002/1‑99 (MQ=37)
                                                gCACCGCTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGGGAGAAGcgcccgc                                                       >  2:3317239/1‑101 (MQ=255)
                                                    cgcttttCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCCGCCCCGCGGGGGGGAAGcgcccggcccg                                                   >  1:1987000/4‑96 (MQ=25)
                                                      ctttcCCCCGCCCGTCGGCGAAACGGCTGTTCCGCGACCGTGTCCTCAAGCGGGACAGGGGCCGGCGCCGGGCCCGGCGGGCGAGAGGCGCCCGCCATggg                                                 <  2:1235879/100‑1 (MQ=25)
                                                        tttCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGACCCGCCCGGCGGGCGAAAAACGCCCGCCAAcggcg                                               >  2:3615884/1‑101 (MQ=37)
                                                          tCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCCGCGGGCGGGAAGCGCCCCCCCCggggggg                                             >  1:2078107/1‑101 (MQ=37)
                                                          tCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCCGCCCGGCGGGGGGGAAACCCCCCCCCAcggcggg                                             >  1:2121187/1‑101 (MQ=37)
                                                            cccgcccgTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGCCGGGGGGAAAGCGCCCGCCACgggcggggg                                           >  2:920494/1‑100 (MQ=37)
                                                            cccgcccgCGGGCGCAAGGGCTGTTCCGCGCCCGTGTCCTCAAGCGTGACAGGGGCCGGCGCCCGGCCGGGGGGGCGAGAAGCGCCCGCCATgggcggggc                                           <  1:1735484/101‑1 (MQ=21)
                                                            cccccccgGGGGCGAAACGGCTGTTCCGCGACCGTGTCCTCAGGCGTGACAGAGGCCGGCGCCCGGCCGGGCGGGCGAGAAGCGCCCGCCATgggcggggc                                           <  1:282867/101‑1 (MQ=18)
                                                             ccgccggCCGGCGCACCGGCTGTTCCGCGCCCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCGGGCGGGCGAGAAGCGCCCGCCATgggcggggcg                                          <  2:3843313/101‑1 (MQ=21)
                                                               gcccgTCGGCGCAACGGCTGTTCCGCGCCCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCGGGCGGGGGAGAAGCGCCCGCCATGGgcggggcggg                                        <  2:2928565/101‑1 (MQ=255)
                                                               gcccgTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGGGAAGCGCCCGCCCCGGgcgggggggg                                        >  1:2747459/1‑101 (MQ=25)
                                                                cccgGCGGCGCAAGGGCTGTCCCGCGCCCGTGCCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGgggcgggc                                       <  1:1592706/101‑1 (MQ=18)
                                                                   gTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGGGAAACGCCCGCCATGGGCGGGGCGGgcg.cg                                   >  2:4477927/1‑100 (MQ=255)
                                                                   gTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCCGGCCGGGGGGCGAGAAGCCCCCGCCAGCGGCGGGGGGGGCG.Ct                                   >  1:4303982/1‑101 (MQ=25)
                                                                    tagtaccAACGGCTGCTCCGCCTCCGTGTCGTTAAGGTTGAAAAAGGCGGTCCCCCGGCCCGGCGGGCAAGGCTCTCCCGCCATGGGCGGGGCGGGCG.CTg                                  <  1:61489/95‑1 (MQ=9)
                                                                      ggCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCCTGGGCGGGGCGGGCG.CGGGc                                >  2:5460484/1‑101 (MQ=255)
                                                                         gcaccgtgtGTTCCGCGCCCGTGTCCTCAAGGGGGACAGAGGCCGGCGCCCGCCCCGGGGGGGGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.gg                            <  1:1317827/93‑1 (MQ=17)
                                                                          caccGGCTTTTCCGCGACCGTGTCCCCAAGCGGGACAGAGGCCGGCCCCCGGCCGGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.ggg                           <  1:5321825/98‑1 (MQ=17)
                                                                           acgggCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.gggg                          <  1:5468200/98‑1 (MQ=255)
                                                                           aaCGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCCCCATGGGCGGGGGGGGGG.CTGGCC.cggg                          >  2:1885854/1‑101 (MQ=21)
                                                                            aCGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGGGAAGGGCCCGCCATGGGCGGGGGGGGGG.GCGGCC.ggggg                         >  1:2921982/1‑101 (MQ=21)
                                                                              ggCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCCTGGGCGGGGCGGGCG.CTGGCC.GGGGGcc                       >  1:4785026/1‑101 (MQ=255)
                                                                              ggCTGTCCCGCGACCGTCCCCTCAAGCGTGACAGGGGCCGGCCCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.GGGGGcc                       <  2:1718412/101‑1 (MQ=18)
                                                                              cggttcaccGCGACCGTTTCCCCAAGTGCAACAGAGGCCAGCGCCCGACCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGGGCCGGGCC.gggggt                        >  1:3986579/8‑100 (MQ=17)
                                                                               gCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGGG.CTGGCC.GGGGGccc                      >  1:3059650/1‑100 (MQ=255)
                                                                               gCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCCCCCTCGGCGGGGAGGGGG.CGGGCC.CGGGGGCt                      >  2:822185/1‑101 (MQ=18)
                                                                               gCTGTTCCGCCCCCGTGTCCTCAAGGGGGACAGGGGCCGGCCCCCGGCCCGGGGGGCGAGAAGCCCCCGCCATGGGCGGGGCGGGCG.CTGGCC.GGGGGCCt                      <  1:5379021/101‑1 (MQ=17)
                                                                                ccgTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCCTGGGCGGGGCGGGCG.CCGGCC.ggggggcct                     >  1:451444/3‑97 (MQ=18)
                                                                                cTGTTCCGCGACGGTGCCCTAAGGCGTGACGGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.GGGGGCCtt                     <  2:1499846/101‑1 (MQ=17)
                                                                                 tttcccGCGCCGGTGCCCTCAAGGGTGACGGGGCCGGGCCCCCGGCCGGGGGGGCGAAAGCCCCCCGCCATGGGCGGGTCGGGCG.CTGGCC.GGGGGCCttt                    <  2:5318813/97‑1 (MQ=12)
                                                                                 tGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGGGCAGGGCCCGCCAGGGGCGGGGCGGGGG.CTCGCC.CGGGGCCtct                    >  1:4386811/1‑99 (MQ=18)
                                                                                  gtcccGCGCCCGTGCCCTCAAGCGTGACAGGGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.GGGGGCCTTTg                   <  2:3478437/98‑1 (MQ=17)
                                                                                  gTTCCGCGCCCGTGCCCTCAGGCGGGACAGGGGCCGGCGCCCGGCCGGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.GGGGGCCTTTg                   <  2:2906626/101‑1 (MQ=17)
                                                                                  gTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CCGGGC.GGGGGGCCTTg                   >  2:2164378/1‑101 (MQ=18)
                                                                                  gTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCCGCCCGGCGGGCGAGAAGCGCCCCCCCCGGGGGGGGCGGGGG.GGGGGC.GGGGGGCCTTg                   >  1:24668/1‑101 (MQ=18)
                                                                                    cccGCGGCCGTGCCCCCAAGCGTGACAGAGGCCGGCCCCCGGCCGGGCGGGGGAGAAGCCCCCGCCATGGGCGGGGCGGGCG.CTGGCC.GGGGGCCTTTggg                 <  1:4850569/100‑1 (MQ=14)
                                                                                     ccGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCCTGGGGGGGGGGGGCG.CGCGCC.CGGGGCCCTTgggg                >  2:38095/1‑100 (MQ=17)
                                                                                         gACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGGGGGCG.CCGGCC.GGGGGCCTTTTGGCGCcg            >  2:3336579/1‑101 (MQ=17)
                                                                                         gACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCAAGGGCGGGGCGGGCG.CCGGGC.GGGGGCCCTTTGGCCCcg            >  2:513611/1‑101 (MQ=17)
                                                                                         ccgatcttCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.GGGGGCCTTTGGGCCCcg            <  2:2550721/94‑1 (MQ=255)
                                                                                           cggtgtCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.GGGGGCCTTTGGGCCCcggc          <  2:4612151/99‑1 (MQ=255)
                                                                                           ccGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CCGGCCGGGGGGCCTTGGGGCCCcgc           >  2:505878/1‑100 (MQ=17)
                                                                                              ggtCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.GGGGGCCTTTGGGCCCCGGCggg       <  2:1662092/100‑1 (MQ=255)
                                                                                                  cTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.GGGGGCCCTTGGGCCCCGGCGGGGCtt   >  2:5410086/1‑101 (MQ=255)
                                                                                                   tCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CCGGCC.CGGGGCCTTTGGGCCCCCGCcgggcctg  >  2:5500155/1‑96 (MQ=14)
                                                                                                   |                                                                                                      
CCGCTGGCAAGGGGCGGCGATGTGCCGGGTTCGGGCGGATCAGGAAAAGCACCGTTTTTCCCCGCCCGTCGGCGCAACGGCTGTTCCGCGACCGTGTCCTCAAGCGTGACAGAGGCCGGCGCCCGGCCCGGCGGGCGAGAAGCGCCCGCCATGGGCGGGGCGGGCG.CTGGCC.GGGGGCCTTTGGGCCCCGGCGGGTCTTG  >  NC_014034/460709‑460908

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG