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.4% 2056.4 / NA 527P346P (CCC→CCG) PA14_RS17955MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base G (237/286);  total (239/288)
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 533 total aligned reads displayed.

GATCGTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCCGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCGGCTTCGT  >  NC_008463/3933205‑3933498
                                                                                                                                                   |                                                                                                                                                  
gATCGTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCggc                                                                                                                                                  <  1:3909302/150‑1 (MQ=255)
 aTCGTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCggcg                                                                                                                                                 >  1:5090725/1‑150 (MQ=255)
  tCGTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGa                                                                                                                                                >  2:3430372/1‑150 (MQ=255)
    gTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACCAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATg                                                                                                                                              >  1:3305392/1‑150 (MQ=255)
     ttCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCCGCGATGa                                                                                                                                             <  2:9321445/150‑1 (MQ=255)
       cccACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGacc                                                                                                                                           >  1:6260523/1‑150 (MQ=255)
       cccACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGacc                                                                                                                                           >  1:2241440/1‑150 (MQ=255)
        ccACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGacca                                                                                                                                          >  2:7475085/1‑150 (MQ=255)
          aCCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGAACCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGaccacc                                                                                                                                        <  2:421970/150‑1 (MQ=255)
            cGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCgg                                                                                                                                      <  1:7805772/150‑1 (MQ=255)
              ccTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGgca                                                                                                                                    <  1:9052878/150‑1 (MQ=255)
               cTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGgcag                                                                                                                                   >  1:2408516/1‑150 (MQ=255)
                 tCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGgcagca                                                                                                                                 >  1:8343652/1‑150 (MQ=255)
                  cGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCAt                                                                                                                                <  1:1736343/150‑1 (MQ=255)
                  cGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCGTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCAt                                                                                                                                <  2:3707515/150‑1 (MQ=255)
                    aGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTc                                                                                                                              >  1:6849669/1‑150 (MQ=255)
                     ggTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTcc                                                                                                                             <  2:5855312/150‑1 (MQ=255)
                          ccGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAg                                                                                                                        <  1:4725357/150‑1 (MQ=255)
                          ccGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAg                                                                                                                        <  1:10575937/150‑1 (MQ=255)
                           cGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGt                                                                                                                       >  1:2821009/1‑150 (MQ=255)
                           cGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGt                                                                                                                       <  1:9701401/150‑1 (MQ=255)
                           cGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGt                                                                                                                       <  1:1613490/150‑1 (MQ=255)
                              ccTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGt                                                                                                                    <  1:1489308/150‑1 (MQ=255)
                                    cGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCct                                                                                                              <  1:993497/150‑1 (MQ=255)
                                    cGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCct                                                                                                              >  2:3296936/1‑150 (MQ=255)
                                     ggCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCctg                                                                                                             >  1:8115341/1‑150 (MQ=255)
                                       cccGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCctgct                                                                                                           <  2:1389138/150‑1 (MQ=255)
                                       cccGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCAGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCctgct                                                                                                           >  2:6120147/1‑150 (MQ=255)
                                          gATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGtt                                                                                                        <  1:5699967/150‑1 (MQ=255)
                                          gATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGtt                                                                                                        >  1:10021591/1‑150 (MQ=255)
                                           aTCCTGATCGGCTTCGCCTCGGACACGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGACTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTc                                                                                                       >  2:11024310/1‑150 (MQ=255)
                                            tCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCa                                                                                                      <  2:6369800/150‑1 (MQ=255)
                                               tGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCatga                                                                                                   <  1:10155295/150‑1 (MQ=255)
                                                 aTCGGCTTCGCCTCGGACAAGCTGTTCGGCGCCCGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCatgatg                                                                                                 >  2:3149499/1‑150 (MQ=255)
                                                 aTCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCatgatg                                                                                                 >  2:7981171/1‑150 (MQ=255)
                                                  tCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATCGTGGTCGGCCTGCTGTTCATGATgg                                                                                                <  2:301949/150‑1 (MQ=255)
                                                    ggCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGa                                                                                              <  1:1322829/150‑1 (MQ=255)
                                                     gCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGAc                                                                                             <  2:4964276/150‑1 (MQ=255)
                                                       ttCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTg                                                                                           >  1:9849377/1‑150 (MQ=255)
                                                         cGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGAc                                                                                         <  1:10743202/150‑1 (MQ=255)
                                                            cTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCt                                                                                      <  2:1833325/150‑1 (MQ=255)
                                                            cTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCt                                                                                      >  2:2376368/1‑150 (MQ=255)
                                                               ggACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTa                                                                                   >  1:1254526/1‑150 (MQ=255)
                                                               ggACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTa                                                                                   >  1:7683019/1‑150 (MQ=255)
                                                                  cAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACgg                                                                                >  2:4363375/1‑150 (MQ=255)
                                                                     gCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGccc                                                                             <  2:9886595/150‑1 (MQ=255)
                                                                       tGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCgg                                                                           >  2:10300083/1‑150 (MQ=255)
                                                                         ttCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGAc                                                                         <  2:5098590/150‑1 (MQ=255)
                                                                         ttCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGAc                                                                         <  1:8546233/150‑1 (MQ=255)
                                                                           cggcgcacggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTc                                                                       >  1:6101548/1‑150 (MQ=255)
                                                                             gcgcacggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCga                                                                     >  1:10848411/1‑150 (MQ=255)
                                                                              cgcacggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCgat                                                                    <  1:8427596/150‑1 (MQ=255)
                                                                              cgcacggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCgat                                                                    <  1:8692315/150‑1 (MQ=255)
                                                                              cgcacggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCgat                                                                    <  1:4524967/150‑1 (MQ=255)
                                                                                cacggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGGTCAGGATGGGACTGACCCTGTACGGCCCGGACTCgatga                                                                  >  2:8552414/1‑150 (MQ=255)
                                                                                   ggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCGCGGCGATGACCACCGGCAGCATTCCGATTGTGGTCGGCCTTCTGTTCATGATGGGAATGACCCTGTTCGGCCCGGACTCGAGGTTCt                                                               <  2:7110358/150‑1 (MQ=255)
                                                                                   ggcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCt                                                               <  2:9019658/150‑1 (MQ=255)
                                                                                    gcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTg                                                              <  1:5945155/150‑1 (MQ=255)
                                                                                    gcgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTg                                                              >  1:6523613/1‑150 (MQ=255)
                                                                                     cgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGGTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGc                                                             >  1:10983397/1‑150 (MQ=255)
                                                                                     cgcaTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGc                                                             >  2:4891890/1‑150 (MQ=255)
                                                                                       caTGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCgg                                                           <  1:2353784/150‑1 (MQ=255)
                                                                                          gCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTc                                                        <  2:1607208/150‑1 (MQ=255)
                                                                                          gCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTc                                                        <  1:3668912/150‑1 (MQ=255)
                                                                                           ccGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTcc                                                       <  2:4633582/150‑1 (MQ=255)
                                                                                            cGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTccg                                                      <  2:9499291/150‑1 (MQ=255)
                                                                                            cGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACGGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTccg                                                      <  2:10592625/150‑1 (MQ=255)
                                                                                             ggCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTccgc                                                     <  1:213869/150‑1 (MQ=255)
                                                                                             ggCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTccgc                                                     <  1:9469460/150‑1 (MQ=255)
                                                                                               ccTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTccgccg                                                   >  2:6465444/1‑150 (MQ=255)
                                                                                               ccTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTccgccg                                                   >  1:7396840/1‑150 (MQ=255)
                                                                                                  gCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCa                                                <  1:7245404/150‑1 (MQ=255)
                                                                                                       aTCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGAc                                           <  2:2773750/150‑1 (MQ=255)
                                                                                                        tCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACt                                          >  2:1974357/1‑150 (MQ=255)
                                                                                                         cAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACtt                                         <  1:4885274/150‑1 (MQ=255)
                                                                                                           gCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCg                                       >  1:10300090/1‑150 (MQ=255)
                                                                                                            ccTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCgg                                      <  2:6871094/150‑1 (MQ=255)
                                                                                                               gctgGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCAc                                   >  1:5357350/1‑150 (MQ=255)
                                                                                                                ctgGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCAcc                                  <  2:8283537/150‑1 (MQ=255)
                                                                                                                 tgGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCg                                 <  1:1971034/150‑1 (MQ=255)
                                                                                                                 tgGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCg                                 <  1:4137463/150‑1 (MQ=255)
                                                                                                                  ggCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCgg                                >  2:7753662/1‑150 (MQ=255)
                                                                                                                  ggCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCgg                                >  2:2618070/1‑150 (MQ=255)
                                                                                                                  ggCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACAACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCgg                                >  1:6349099/1‑150 (MQ=255)
                                                                                                                   gcgcTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGc                               <  1:5475250/150‑1 (MQ=255)
                                                                                                                     gcTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGa                             <  2:2125386/150‑1 (MQ=255)
                                                                                                                      cTGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAg                            >  2:8749031/1‑150 (MQ=255)
                                                                                                                       tGACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTC‑GCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGgc                          <  1:3068004/150‑1 (MQ=255)
                                                                                                                        gACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGgc                          <  2:8877448/150‑1 (MQ=255)
                                                                                                                        gACCGTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGgc                          <  1:2530848/150‑1 (MQ=255)
                                                                                                                            gTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGgccgcc                      <  1:6726639/150‑1 (MQ=255)
                                                                                                                            gTCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGgccgcc                      >  1:4362464/1‑150 (MQ=255)
                                                                                                                             tCTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGccgccg                     >  1:3116608/1‑150 (MQ=255)
                                                                                                                              cTGCCTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCGGCCgg                    >  2:5605142/1‑150 (MQ=255)
                                                                                                                                  cTGGCGCTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCAcc                <  2:4498902/150‑1 (MQ=255)
                                                                                                                                       gcTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTc           <  2:7220976/150‑1 (MQ=255)
                                                                                                                                        cTGTTCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTcc          >  2:4063990/1‑150 (MQ=255)
                                                                                                                                           ttCGTCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCGGc       >  2:10501131/1‑150 (MQ=255)
                                                                                                                                               tCCCGGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCGGCTTCg   >  2:7394380/1‑150 (MQ=255)
                                                                                                                                                cccgGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCGGCTTCGt  <  1:3748026/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GATCGTTCCCACCGCCTTCGAGGTAGCCGGCCTGCTCGGCCCGATCCTGATCGGCTTCGCCTCGGACAAGCTGTTCGGCGCACGGCGCATGCCGGCCTGCGTGATCAGCCTGCTGGCGCTGACCGTCTGCCTGGCGCTGTTCGTCCCCGCGATGACCACCGGCAGCATTCCGATAGTGGTCGGCCTGCTGTTCATGATGGGACTGACCCTGTACGGCCCGGACTCGATGATCTGCGGCTCCGCCGCCATCGACTTCGGCACCGGCGAGGCCGCCGGCACCGCCTCCGGCTTCGT  >  NC_008463/3933205‑3933498

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