Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,659,831 Δ1 bp intergenic (+71/+28) alpR → / ← PA14_RS21395 transcriptional regulator AlpR/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,659,8280G.93.9% 1434.1 / 135.6 357intergenic (+68/+31)alpR/PA14_RS21395transcriptional regulator AlpR/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (141/180);  minor base A (9/12);  total (165/192)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 389 total aligned reads displayed.

GGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAAAAAAATCCCGGGCAAAAAAAAAGGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACGCTTG  >  NC_008463/4659679‑4659974
                                                                                                                                                     |                                                                                                                                                  
ggACGAGCGCTACGAGCCCCAGGCAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGGTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAg                                                                                                                                                    <  2:4557938/150‑1 (MQ=255)
  aCGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAggg                                                                                                                                                  >  2:9741238/1‑150 (MQ=255)
   cGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggc                                                                                                                                                 >  2:9135654/1‑149 (MQ=255)
    gAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCACCGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcc                                                                                                                                                <  2:2076412/150‑3 (MQ=255)
    gAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcc                                                                                                                                                >  2:6131043/1‑148 (MQ=255)
    gAGCGCGACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAgggcc                                                                                                                                                <  2:4786028/150‑3 (MQ=255)
         cTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCAt                                                                                                                                          >  2:2240787/1‑150 (MQ=255)
          tACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCAtt                                                                                                                                         <  2:8729217/150‑1 (MQ=255)
           aCGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtg                                                                                                                                        >  2:5007400/1‑150 (MQ=255)
            cGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtgc                                                                                                                                       >  1:9908141/1‑150 (MQ=255)
            cGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATtgc                                                                                                                                       >  1:3685151/1‑150 (MQ=255)
                ccccAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCa                                                                                                                                   <  1:7194948/150‑1 (MQ=255)
                        aGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCttt                                                                                                                           <  1:6421228/150‑1 (MQ=255)
                        aGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCttt                                                                                                                           >  1:9319795/1‑150 (MQ=255)
                        aGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGGTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCttt                                                                                                                           <  1:10025663/150‑1 (MQ=255)
                         gccgccGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAACAGGGCCTGCATTGCTGCAGGCCCTTTc                                                                                                                          <  1:10583963/150‑1 (MQ=255)
                         gccgccGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTc                                                                                                                          <  1:1767457/150‑1 (MQ=255)
                           cgccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCgg                                                                                                                        <  2:3782633/150‑1 (MQ=255)
                             ccgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTa                                                                                                                      <  1:3629836/150‑1 (MQ=255)
                              cgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTAt                                                                                                                     >  2:7115024/1‑150 (MQ=255)
                              cgAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTAt                                                                                                                     >  2:5326919/1‑150 (MQ=255)
                               gAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCACGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAATGAAAAAGGGCCTGCATTGCTGCAGGCCCGGTCGGTATg                                                                                                                    >  2:8427797/1‑150 (MQ=255)
                                     aTCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCAc                                                                                                              <  2:9270590/150‑1 (MQ=255)
                                      tCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACt                                                                                                             >  2:8060901/1‑150 (MQ=255)
                                         aCGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCag                                                                                                          <  1:9714563/150‑1 (MQ=255)
                                          cGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCagg                                                                                                         <  1:9450568/150‑1 (MQ=255)
                                          cGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCagg                                                                                                         >  2:7258628/1‑150 (MQ=255)
                                              aTCGGCCGGGTGTTCTGGTACTCGGAGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAggagga                                                                                                     <  1:2790026/150‑1 (MQ=255)
                                                 ggCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAATTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTc                                                                                                  >  1:10407/1‑150 (MQ=255)
                                                 ggCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTc                                                                                                  <  1:8664640/150‑1 (MQ=255)
                                                    cGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGGTCCGGCGGAAGAAAAATCCCGGGCAAAGATAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGaa                                                                                               <  2:3219160/150‑1 (MQ=255)
                                                      ggTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAAcc                                                                                             <  1:10081865/150‑1 (MQ=255)
                                                             tGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGAcc                                                                                      <  1:1243891/150‑1 (MQ=255)
                                                               gTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGc                                                                                    >  1:3896935/1‑150 (MQ=255)
                                                                  cTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGAGTCGAACCTCCGACCGCTCg                                                                                 >  2:8800384/1‑150 (MQ=255)
                                                                     ggTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCACAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGtt                                                                              <  2:10517419/150‑1 (MQ=255)
                                                                     ggTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGtt                                                                              <  1:4392471/150‑1 (MQ=255)
                                                                     ggTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGtt                                                                              >  1:8038646/1‑150 (MQ=255)
                                                                      gTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTc                                                                             <  1:6062677/150‑1 (MQ=255)
                                                                      gTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTc                                                                             <  2:2584133/150‑1 (MQ=255)
                                                                        gCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGt                                                                           >  1:9870367/1‑150 (MQ=255)
                                                                               tGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAg                                                                    >  1:6998099/1‑150 (MQ=255)
                                                                               tGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAg                                                                    <  1:4416504/150‑1 (MQ=255)
                                                                               tGACCCCGTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAg                                                                    <  2:7645820/150‑1 (MQ=255)
                                                                                gACCCCGTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGt                                                                   <  2:6700130/150‑1 (MQ=255)
                                                                                 aCCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTa                                                                  <  2:8895711/150‑1 (MQ=255)
                                                                                  ccccTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTAc                                                                 <  2:3394060/150‑1 (MQ=255)
                                                                                    ccTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTTGCCGAGTActc                                                               >  1:8235083/1‑150 (MQ=255)
                                                                                       gcccgccccgcCTGACGACGGGGCCGCGGTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCAGTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGATGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTActctct                                                            <  2:7471476/149‑3 (MQ=37)
                                                                                        cccgccccgcCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                           <  1:4254965/150‑1 (MQ=255)
                                                                                        cccgccccgcCTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCAGAGTACTCTATc                                                           >  2:7949610/1‑150 (MQ=255)
                                                                                         ccgccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAACAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATcc                                                          <  2:124891/150‑1 (MQ=255)
                                                                                         ccgccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATcc                                                          <  1:7598885/150‑1 (MQ=255)
                                                                                          cgccccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCa                                                         >  1:6287159/1‑150 (MQ=255)
                                                                                             cccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                      >  1:1505636/1‑150 (MQ=255)
                                                                                             cccgccTGACGACGGGGCCGCGGTCCGGCGGAAGACAAATCCCGGGCAACGAAACAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                      <  2:10627523/150‑1 (MQ=25)
                                                                                              ccgccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTg                                                     >  2:512622/1‑150 (MQ=255)
                                                                                                gccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAACAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAg                                                   <  1:4973085/150‑1 (MQ=255)
                                                                                                gccTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAg                                                   <  1:5190615/150‑1 (MQ=255)
                                                                                                  cTGACGACGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAATGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                 >  1:4621489/1‑150 (MQ=255)
                                                                                                    gacgacGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                               >  2:1759830/1‑150 (MQ=255)
                                                                                                    gacgacGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCAGAGCTaa                                               >  1:2650220/1‑149 (MQ=25)
                                                                                                    gacgacGGGGCCGCGTTCCGGCGGAAGAAAAATCCCGCGCAAAGAAACAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                               <  2:6487817/150‑1 (MQ=25)
                                                                                                      cgacgGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                             >  1:5080695/1‑150 (MQ=255)
                                                                                                        acgGGGCCGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGt                                           <  2:4225549/150‑1 (MQ=255)
                                                                                                              ccGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAACAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCgtga                                     <  1:4878368/150‑2 (MQ=21)
                                                                                                              ccGCGTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCgtgt                                     <  2:769850/150‑1 (MQ=255)
                                                                                                                 cgTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCgtcgtgg                                  <  1:2328309/150‑6 (MQ=21)
                                                                                                                 cgTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCgg                                  <  1:6185523/150‑1 (MQ=255)
                                                                                                                  gTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGt                                 <  1:2157656/150‑1 (MQ=255)
                                                                                                                  gTTCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGt                                 <  1:1644786/150‑1 (MQ=255)
                                                                                                                   ttCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTg                                >  2:128274/1‑150 (MQ=255)
                                                                                                                   ttCCGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTg                                >  1:3228604/1‑150 (MQ=255)
                                                                                                                     ccGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCt                              <  1:7655315/150‑1 (MQ=255)
                                                                                                                     ccGGCGGAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCt                              <  1:6595254/150‑1 (MQ=255)
                                                                                                                       ggcggAAGAAACATCCCGGGCAAAGAAAAAGCGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCttt                            <  1:10566712/150‑1 (MQ=21)
                                                                                                                       ggcggAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCttt                            <  1:2462671/150‑1 (MQ=255)
                                                                                                                        gcggAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGGCGGTGCtttt                           >  2:10016787/1‑150 (MQ=255)
                                                                                                                         cggAAGAAAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCttatt                          >  1:5523538/1‑150 (MQ=255)
                                                                                                                           gaagaaAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTcc                        <  1:8399021/150‑1 (MQ=255)
                                                                                                                           gaagaaAAATCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGGCGGTGCTGTTTcc                        >  1:10138115/1‑150 (MQ=21)
                                                                                                                           gaagaaAAATCCCGGGCAAAGAAAAAGGGCCGGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCTGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTcc                        <  1:963586/150‑1 (MQ=21)
                                                                                                                              gaaaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTACGGc                     >  1:227453/2‑150 (MQ=255)
                                                                                                                                aaaaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCTGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCa                   <  2:10114826/150‑1 (MQ=255)
                                                                                                                                  aaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGTGTGCGTGTCGGTGCTTTTTCCGGCCAcc                 <  1:1970837/150‑1 (MQ=255)
                                                                                                                                  aaTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCAcc                 <  1:9261844/150‑1 (MQ=255)
                                                                                                                                   aTCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGATTTTTCCGGCCACCt                >  1:7411912/1‑150 (MQ=255)
                                                                                                                                   aTCCCGGGCAAAGAAAATGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCt                >  1:5687172/1‑150 (MQ=255)
                                                                                                                                    tCCCGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTg               >  2:2400026/1‑150 (MQ=255)
                                                                                                                                      ccGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGcc             >  1:4804549/1‑150 (MQ=255)
                                                                                                                                      ccGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCATGATCCGGCCACCTggc             >  2:1814174/1‑148 (MQ=21)
                                                                                                                                       cGGGCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCg            >  1:9047025/1‑150 (MQ=255)
                                                                                                                                        gggCAAAGAAAAAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCgg           <  1:4725478/150‑1 (MQ=255)
                                                                                                                                            caaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTaa        <  2:6947652/149‑1 (MQ=255)
                                                                                                                                            caaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTaa        <  2:2879734/149‑1 (MQ=255)
                                                                                                                                             aaagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAAc       >  2:5647145/1‑150 (MQ=255)
                                                                                                                                              aagaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACg      <  1:10205308/147‑1 (MQ=255)
                                                                                                                                               agaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCGTGTTCCGGCCAACTGCCGGGAACGc     >  2:302299/3‑150 (MQ=21)
                                                                                                                                                gaaaaAGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACGCt    <  2:151656/149‑1 (MQ=255)
                                                                                                                                                  aaaaGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCAATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACGCTTg  <  1:5265337/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGACGAGCGCTACGAGCCCCAGGAAGCCGCCGAGCACATCAACGTGATCGGCCGGGTGTTCTGGTACTCGGTGCTGGTCTGACCCCTTCCCGCCCCGCCTGACGACGGGGCCGCGTTCCGGCGGAAAAAAAATCCCGGGCAAAAAAAAAGGGGCCTGCATTGCTGCAGGCCCTTTCGGTATGGCGCACTCAGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTGCGTGTCGGTGCTTTTTCCGGCCACCTGCCGGTAACGCTTG  >  NC_008463/4659679‑4659974

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: