Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,670,537 Δ1 bp coding (537/1020 nt) PA14_RS16675 → ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,670,5340T.99.7% 1687.1 / ‑2.3 343coding (534/1020 nt)PA14_RS16675ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (156/186);  total (156/187)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 354 total aligned reads displayed.

AGCCCGCGACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTTCCTGGAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCC  >  NC_008463/3670384‑3670679
                                                                                                                                                      |                                                                                                                                                   
ggCCCGCGACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGt                                                                                                                                                     <  1:8243599/150‑1 (MQ=255)
       gACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTgg                                                                                                                                             >  2:3238566/1‑151 (MQ=255)
       gACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGAGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGtttccggg                                                                                                                                              >  2:1404674/1‑146 (MQ=255)
             cTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGaagtct                                                                                                                                      <  2:6305261/151‑6 (MQ=255)
              tCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTg                                                                                                                                     <  2:3866618/151‑1 (MQ=255)
              tCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTg                                                                                                                                     <  2:2767075/151‑1 (MQ=255)
              tCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTg                                                                                                                                     >  1:5618281/1‑151 (MQ=255)
                gCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGt                                                                                                                                   <  2:84400/151‑1 (MQ=255)
                gCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGt                                                                                                                                   >  2:1696406/1‑151 (MQ=255)
                 cccGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTc                                                                                                                                  <  2:4575621/151‑1 (MQ=255)
                  ccGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCg                                                                                                                                 >  1:9283732/1‑151 (MQ=255)
                  ccGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTATCCTGGAAGTCTGGTCg                                                                                                                                 >  1:10548358/1‑151 (MQ=255)
                   cGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCgg                                                                                                                                >  2:820183/1‑151 (MQ=255)
                    gggTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCggg                                                                                                                               >  2:232825/1‑151 (MQ=255)
                     ggTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCgggg                                                                                                                              <  1:6186590/151‑1 (MQ=255)
                     ggTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCgggg                                                                                                                              >  1:340870/1‑151 (MQ=255)
                       tGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCt                                                                                                                            <  1:32977/151‑1 (MQ=255)
                            tACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCgg                                                                                                                       <  2:4252893/151‑1 (MQ=255)
                             aCGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGt                                                                                                                      <  2:4745130/151‑1 (MQ=255)
                              cGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGCCCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTg                                                                                                                     <  2:6410586/151‑1 (MQ=255)
                              cGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTg                                                                                                                     >  2:10133538/1‑151 (MQ=255)
                               ggCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGc                                                                                                                    <  1:870813/151‑1 (MQ=255)
                               ggCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGc                                                                                                                    >  2:6019987/1‑151 (MQ=255)
                                gCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGct                                                                                                                   <  1:9445794/151‑1 (MQ=255)
                                  ttCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGctct                                                                                                                 >  2:5156570/1‑151 (MQ=255)
                                      gcgTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACct                                                                                                             <  2:10084531/151‑1 (MQ=255)
                                      gcgTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACct                                                                                                             <  1:5980756/151‑1 (MQ=255)
                                      gcgTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACct                                                                                                             <  2:1809467/151‑1 (MQ=255)
                                        gtgtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgc                                                                                                           <  1:2589517/151‑1 (MQ=255)
                                         tgtCGATCCTCTTCGCCCTGACCCTGTCCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCAGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgct                                                                                                          >  2:5899320/1‑151 (MQ=255)
                                         tgtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgct                                                                                                          <  1:4989895/151‑1 (MQ=255)
                                         tgtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgct                                                                                                          <  2:2686929/151‑1 (MQ=255)
                                          gtCGATCCTCATCGCCCTGACCCGGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgctg                                                                                                         <  2:10199146/151‑1 (MQ=255)
                                            cGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGat                                                                                                       >  1:3358673/1‑151 (MQ=255)
                                            cGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGat                                                                                                       >  1:9726954/1‑151 (MQ=255)
                                            cGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGat                                                                                                       <  1:3910368/151‑1 (MQ=255)
                                            cGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTGTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTCCCTGCTGat                                                                                                       >  1:2253246/1‑151 (MQ=255)
                                              aTCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatca                                                                                                     >  1:4602459/1‑151 (MQ=255)
                                               tCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatcat                                                                                                    >  2:540525/1‑151 (MQ=255)
                                               tCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatcat                                                                                                    <  1:7394139/151‑1 (MQ=255)
                                                ccTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatcatc                                                                                                   <  2:7130557/151‑1 (MQ=255)
                                                ccTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatcatc                                                                                                   <  1:4136733/151‑1 (MQ=255)
                                                 ctctTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATcc                                                                                                  >  1:2989181/1‑151 (MQ=255)
                                                   ctTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTg                                                                                                >  2:1328516/1‑151 (MQ=255)
                                                    ttCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTgg                                                                                               >  2:9473947/1‑151 (MQ=255)
                                                       gCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCa                                                                                            <  2:7861078/151‑1 (MQ=255)
                                                       gCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCa                                                                                            <  1:6725146/151‑1 (MQ=255)
                                                         cctgaccctgacCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGc                                                                                          >  2:3740840/1‑151 (MQ=255)
                                                          ctgaccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCt                                                                                         <  1:2017109/151‑1 (MQ=255)
                                                                ctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTcc                                                                                   <  1:8604007/151‑1 (MQ=255)
                                                                 tgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCa                                                                                  <  2:9115303/151‑1 (MQ=255)
                                                                 tgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCa                                                                                  <  1:267422/151‑1 (MQ=255)
                                                                 tgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCa                                                                                  <  2:8922694/151‑1 (MQ=255)
                                                                 tgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCa                                                                                  >  1:761825/1‑151 (MQ=255)
                                                                   accGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGc                                                                                >  1:845327/1‑151 (MQ=255)
                                                                    ccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGcc                                                                               >  1:633250/1‑151 (MQ=255)
                                                                     cGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCg                                                                              <  1:6577325/151‑1 (MQ=255)
                                                                          tCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACtt                                                                         >  2:6856489/1‑151 (MQ=255)
                                                                                  ctGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgc                                                                 >  1:6220492/1‑151 (MQ=255)
                                                                                   tGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgct                                                                >  1:364856/1‑151 (MQ=255)
                                                                                   cgATCGGCGTCACCGTCGGCGCCCTGCAAGTCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCGACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgct                                                                <  2:10509116/150‑1 (MQ=255)
                                                                                     aTCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgctgc                                                              >  1:2545251/1‑151 (MQ=255)
                                                                                     aTCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgctgc                                                              >  1:5831345/1‑151 (MQ=255)
                                                                                       cGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGctgctg                                                            >  1:8107435/1‑151 (MQ=255)
                                                                                        ggCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTgg                                                           <  2:10414737/151‑1 (MQ=255)
                                                                                        ggCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTgg                                                           >  2:416148/1‑151 (MQ=255)
                                                                                            tCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCGGCTGGGcag                                                       >  1:9845058/1‑150 (MQ=255)
                                                                                             cACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAGGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcatc                                                      >  2:2352110/1‑151 (MQ=255)
                                                                                              aCCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcatca                                                     >  1:5207664/1‑151 (MQ=255)
                                                                                               ccGTCGGCGCCCTGCACGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcatcat                                                    <  2:8560003/151‑1 (MQ=255)
                                                                                                 gTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgc                                                  <  1:1637791/151‑1 (MQ=255)
                                                                                                   cGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgctg                                                >  1:10276627/1‑151 (MQ=255)
                                                                                                   cGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgctg                                                <  2:5704512/151‑1 (MQ=255)
                                                                                                    ggCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgctgc                                               >  1:10107762/1‑151 (MQ=255)
                                                                                                     gcgcCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATgctgct                                              <  2:3531137/151‑1 (MQ=255)
                                                                                                       gcCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGt                                            >  2:4355674/1‑151 (MQ=255)
                                                                                                       gcCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGt                                            >  2:7629091/1‑151 (MQ=255)
                                                                                                        cccTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGtt                                           <  2:3076471/151‑1 (MQ=255)
                                                                                                        cccTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGtt                                           <  1:9108771/151‑1 (MQ=255)
                                                                                                        cccTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGtt                                           <  1:4943156/151‑1 (MQ=255)
                                                                                                           tGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTc                                        <  1:2903711/151‑1 (MQ=255)
                                                                                                             cAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCt                                      <  1:6885016/151‑1 (MQ=255)
                                                                                                              aaGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTg                                     <  1:3809894/151‑1 (MQ=255)
                                                                                                                 gCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGAt                                  <  1:8859085/151‑1 (MQ=255)
                                                                                                                   tactacGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGa                                >  1:10409673/1‑151 (MQ=255)
                                                                                                                    actacGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAg                               <  1:10013508/151‑1 (MQ=255)
                                                                                                                    actacGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAg                               >  1:3152517/1‑151 (MQ=255)
                                                                                                                     ctacGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGc                              <  1:8802920/151‑1 (MQ=255)
                                                                                                                        cggcggCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCtg                           >  2:10624640/1‑151 (MQ=255)
                                                                                                                          gcggcTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCtggt                         <  2:402791/151‑1 (MQ=255)
                                                                                                                              cTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGGTTTCCTGGATGAGCCTGGTGGAc                     >  1:4407029/1‑151 (MQ=255)
                                                                                                                                gggTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGt                   <  2:1126717/151‑1 (MQ=255)
                                                                                                                                  gTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTgg                 <  1:1304032/151‑1 (MQ=255)
                                                                                                                                    cGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTc               <  1:1797826/151‑1 (MQ=255)
                                                                                                                                     gACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTcc              <  1:3241064/151‑1 (MQ=255)
                                                                                                                                      aCCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCcg             <  1:6453074/151‑1 (MQ=255)
                                                                                                                                        cTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCcgcg           <  2:5482076/151‑1 (MQ=255)
                                                                                                                                        cTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCcgcg           >  1:9511620/1‑151 (MQ=255)
                                                                                                                                                ccAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtc   <  2:7653758/151‑1 (MQ=255)
                                                                                                                                                 cAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCAGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtcc  >  2:6636749/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
AGCCCGCGACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTTCCTGGAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCC  >  NC_008463/3670384‑3670679

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