Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,933,352 C→G P346P (CCC→CCG)  PA14_RS17955 → MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,933,3520CG99.8% 1877.4 / NA 493P346P (CCC→CCG) PA14_RS17955MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (231/261);  total (232/261)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 508 total aligned reads displayed.

CGATCGTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCCGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCGGCTTCGTC  >  NC_008463/3933204‑3933499
                                                                                                                                                    |                                                                                                                                                   
cgATCGTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCTGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCggc                                                                                                                                                   >  1:2100743/1‑151 (MQ=255)
cgATCGTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCggc                                                                                                                                                   >  2:5999794/1‑151 (MQ=255)
    cGTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATg                                                                                                                                               <  2:9415840/151‑1 (MQ=255)
     gTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGAtta                                                                                                                                              >  1:4091978/1‑149 (MQ=255)
     gTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGa                                                                                                                                              >  1:7835080/1‑151 (MQ=255)
     gTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGa                                                                                                                                              <  1:7507574/151‑1 (MQ=255)
     gTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGa                                                                                                                                              >  1:7063665/1‑151 (MQ=255)
     gTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGa                                                                                                                                              >  1:4283607/1‑151 (MQ=255)
      ttCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGac                                                                                                                                             <  2:5754766/151‑1 (MQ=255)
       tCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGacc                                                                                                                                            <  1:1853261/151‑1 (MQ=255)
        cccACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGacca                                                                                                                                           <  1:4582235/151‑1 (MQ=255)
          cACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGaccacc                                                                                                                                         >  1:8630059/1‑151 (MQ=255)
           cccGCCTTCGAGGTAGCCGGCCGGCTCGGCCCGATCCTGCTCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCg                                                                                                                                        <  2:6974503/150‑1 (MQ=255)
           aCCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCg                                                                                                                                        <  2:8223422/151‑1 (MQ=255)
              gCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGgca                                                                                                                                     >  2:6473446/1‑151 (MQ=255)
              gCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGgca                                                                                                                                     <  1:2911399/151‑1 (MQ=255)
               ccTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGgcag                                                                                                                                    >  1:6280189/1‑151 (MQ=255)
                   cGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCAtt                                                                                                                                <  1:1662263/151‑1 (MQ=255)
                       gTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGa                                                                                                                            >  1:8532647/1‑151 (MQ=255)
                        tAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGAt                                                                                                                           >  2:10585115/1‑151 (MQ=255)
                         aGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATa                                                                                                                          <  2:5329474/151‑1 (MQ=255)
                           ccGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGt                                                                                                                        <  1:2832132/151‑1 (MQ=255)
                                cTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCg                                                                                                                   <  1:311717/151‑1 (MQ=255)
                                 tGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCgg                                                                                                                  <  2:9700468/151‑1 (MQ=255)
                                    tCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCct                                                                                                               <  1:5623214/151‑1 (MQ=255)
                                    tCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCct                                                                                                               <  1:9513637/151‑1 (MQ=255)
                                    tCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTCCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCct                                                                                                               <  2:8918080/151‑1 (MQ=255)
                                        cccGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCctgctg                                                                                                           >  1:1391529/1‑151 (MQ=255)
                                         ccGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGt                                                                                                          >  2:8824262/1‑151 (MQ=255)
                                            aTCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCa                                                                                                       >  1:4991417/1‑151 (MQ=255)
                                               cTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATCGTGGTCGGCCTGCTGTTCatga                                                                                                    >  2:282905/1‑151 (MQ=255)
                                                tGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCatgat                                                                                                   <  2:7767262/151‑1 (MQ=255)
                                                tGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCatgat                                                                                                   >  1:9016880/1‑151 (MQ=255)
                                                 gATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGTCCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGTCCACCGGCTGCATTCCGATAGTGGTCGGCCTGCTGTTCttgttg                                                                                                  >  2:8612552/1‑151 (MQ=25)
                                                 gATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCatgatg                                                                                                  >  2:4129709/1‑151 (MQ=255)
                                                  aTCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGTTAGTGGTCGGCCTGCTGTTCATGATgg                                                                                                 <  2:9352605/151‑1 (MQ=255)
                                                  aTCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATgg                                                                                                 <  1:10809791/151‑1 (MQ=255)
                                                     ggCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGAc                                                                                              <  1:3269712/151‑1 (MQ=255)
                                                      gCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACt                                                                                             >  1:2584931/1‑151 (MQ=255)
                                                       cTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTg                                                                                            >  2:9527670/1‑151 (MQ=255)
                                                        ttCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGa                                                                                           <  2:10047614/151‑1 (MQ=255)
                                                         tCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGAc                                                                                          <  2:3638103/151‑1 (MQ=255)
                                                           gCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGAccc                                                                                        >  2:2999785/1‑151 (MQ=255)
                                                             cTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTg                                                                                      <  2:6414120/151‑1 (MQ=255)
                                                             cTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTg                                                                                      <  1:6682404/151‑1 (MQ=255)
                                                               cGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTa                                                                                    <  1:7360878/151‑1 (MQ=255)
                                                               cGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTa                                                                                    >  2:7265038/1‑151 (MQ=255)
                                                                ggACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTAc                                                                                   >  2:3856181/1‑151 (MQ=255)
                                                                 gACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACg                                                                                  >  1:3798421/1‑151 (MQ=255)
                                                                  aCAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACgg                                                                                 <  1:10385021/151‑1 (MQ=255)
                                                                   cAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGc                                                                                <  2:4044208/151‑1 (MQ=255)
                                                                   cAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGc                                                                                >  2:4319263/1‑151 (MQ=255)
                                                                    aaGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCATCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGGTCGTCCCGGCGATGACCACCGGCATCATTCCGAACGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGcc                                                                               >  2:8039013/1‑151 (MQ=25)
                                                                    aaGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGcc                                                                               >  2:6615253/1‑151 (MQ=255)
                                                                    aaGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGcc                                                                               >  2:2250626/1‑151 (MQ=255)
                                                                    aaGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGcc                                                                               >  2:8295723/1‑151 (MQ=255)
                                                                     aGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGccc                                                                              <  2:560319/151‑1 (MQ=255)
                                                                     aGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGccc                                                                              <  1:7668461/151‑1 (MQ=255)
                                                                      gCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCg                                                                             <  2:7368219/151‑1 (MQ=255)
                                                                       cTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCgg                                                                            >  1:4765201/1‑151 (MQ=255)
                                                                         gTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGAc                                                                          <  2:2076066/151‑1 (MQ=255)
                                                                              gcgcacggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCgat                                                                     <  2:1241535/151‑1 (MQ=255)
                                                                              gcgcacggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGACCGGACTCgat                                                                     >  2:7090397/1‑151 (MQ=255)
                                                                                 cacggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCgatgat                                                                  >  2:8962235/1‑151 (MQ=255)
                                                                                 cacggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCgatgat                                                                  >  1:7177737/1‑151 (MQ=255)
                                                                                  acggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATc                                                                 <  1:9239238/151‑1 (MQ=255)
                                                                                   cggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCt                                                                >  1:9947299/1‑151 (MQ=255)
                                                                                   cggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCt                                                                <  2:3321747/151‑1 (MQ=255)
                                                                                    ggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTg                                                               >  1:10007641/1‑151 (MQ=255)
                                                                                        caTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGc                                                           >  1:8072457/1‑151 (MQ=255)
                                                                                         aTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCt                                                          >  1:641248/1‑151 (MQ=255)
                                                                                         aTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCt                                                          <  2:6816533/151‑1 (MQ=255)
                                                                                           gCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTcc                                                        <  1:4369386/151‑1 (MQ=255)
                                                                                            ccGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTccg                                                       <  1:3573552/151‑1 (MQ=255)
                                                                                             cGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTccgc                                                      >  1:7293258/1‑151 (MQ=255)
                                                                                             cGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTccgc                                                      >  1:10841713/1‑151 (MQ=255)
                                                                                                 cTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACGCGATGAGCTGCGGCTCCgccgcc                                                  >  2:7568840/1‑151 (MQ=255)
                                                                                                     gTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCg                                              <  2:103930/151‑1 (MQ=255)
                                                                                                      tGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGa                                             <  1:5255363/151‑1 (MQ=255)
                                                                                                      tGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGa                                             >  2:10798770/1‑151 (MQ=255)
                                                                                                        aTCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACt                                           <  2:6996241/151‑1 (MQ=255)
                                                                                                           aGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCg                                        >  2:4483304/1‑151 (MQ=255)
                                                                                                                gctgGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCAcc                                   >  1:3494710/1‑151 (MQ=255)
                                                                                                                     cgcTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGa                              <  1:1950847/151‑1 (MQ=255)
                                                                                                                       cTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAgg                            >  2:5181482/1‑151 (MQ=255)
                                                                                                                       cTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAgg                            <  2:2509457/151‑1 (MQ=255)
                                                                                                                         gACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGgcc                          <  2:1357199/151‑1 (MQ=255)
                                                                                                                              tCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCgg                     <  2:10164186/151‑1 (MQ=255)
                                                                                                                               cTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGc                    >  2:6205070/1‑151 (MQ=255)
                                                                                                                                tGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCa                   >  1:5498510/1‑151 (MQ=255)
                                                                                                                                       cgcTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTc            >  2:8140710/1‑151 (MQ=255)
                                                                                                                                         ctgggcGTCCCGGCGATGACCACCGGCAGCATTCCGATCGTGGTCGGCCTGCGGTTCATCATGGGACTGACCCTGCACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACGTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCg          <  1:7877798/146‑1 (MQ=255)
                                                                                                                                         cTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCg          <  1:3196406/151‑1 (MQ=255)
                                                                                                                                         cTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCg          <  2:9844443/151‑1 (MQ=255)
                                                                                                                                            ttCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCGGCt       <  1:10559172/151‑1 (MQ=255)
                                                                                                                                              cGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCGGCTTc     <  2:4927989/151‑1 (MQ=255)
                                                                                                                                               gTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCCCCGGCGAGGCCGCCGGCACCGCCTCCGGCTTCg    <  2:1739698/151‑1 (MQ=255)
                                                                                                                                               gTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCGGCTTCg    <  1:6156831/151‑1 (MQ=255)
                                                                                                                                               gTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCGGCTTCg    >  2:268319/1‑151 (MQ=255)
                                                                                                                                               gTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCGGCTTCg    >  1:340125/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGATCGTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCCGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCGGCTTCGTC  >  NC_008463/3933204‑3933499

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: ‑