Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,700,807 Δ1 bp intergenic (+165/‑150) PA14_RS21585 → / → pbpG aspartate/tyrosine/aromatic aminotransferase/D‑alanyl‑D‑alanine endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,700,8070G.99.8% 2269.2 / ‑3.6 459intergenic (+165/‑150)PA14_RS21585/pbpGaspartate/tyrosine/aromatic aminotransferase/D‑alanyl‑D‑alanine endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (201/257);  total (202/257)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 459 total aligned reads displayed.

AGCCGGCTGTGATGGCCGGCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCGCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCTTAGTGAGCAAC  >  NC_008463/4700661‑4700952
                                                                                                                                                  |                                                                                                                                                 
aGCCGGCTGTGATGGCCGGCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAAc                                                                                                                                              <  1:4830163/151‑1 (MQ=255)
   cGGCTGTGATGGCCGGCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACctc                                                                                                                                           <  2:50852/151‑1 (MQ=255)
    ggCTGTGATGGCCGGCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACctct                                                                                                                                          <  1:6519210/151‑1 (MQ=255)
       tgtgATGGCCGGCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCg                                                                                                                                       >  2:4146293/1‑151 (MQ=255)
         tgATGGCCGGCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTAGTCGCCCTACCGCACCTCGACCGTGTGGGCGTTGGGCGTATTCCGCCCCAGGCGTTATCCGCCGCACCTCGCCGac                                                                                                                                     >  2:7704844/1‑150 (MQ=37)
           aTGGCCGGCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAAcc                                                                                                                                   >  2:9660926/1‑151 (MQ=255)
             ggCCGGCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGt                                                                                                                                 <  2:6029435/151‑1 (MQ=255)
                cGGCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGcc                                                                                                                              >  2:1064261/1‑151 (MQ=255)
                 ggCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCa                                                                                                                             >  1:4652567/1‑151 (MQ=255)
                 ggCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCa                                                                                                                             >  2:10062268/1‑151 (MQ=255)
                  gCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCaa                                                                                                                            >  1:6878265/1‑151 (MQ=255)
                  gCTTGTTTTCGCCCCAGCGACTGCCGCGCACTCCCGCCGAGGCGGCGGTTAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACTGGCGTTATCCGCCGTACCTCTCCGAACCGTGCCaa                                                                                                                            <  1:8259027/151‑1 (MQ=255)
                    tttttttCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAAcc                                                                                                                          <  1:515875/151‑1 (MQ=255)
                    tttttttCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAAcc                                                                                                                          >  2:1828780/1‑151 (MQ=255)
                       ttttCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTGGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATa                                                                                                                       >  2:4603759/1‑151 (MQ=255)
                        tttCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATAGCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATaa                                                                                                                      >  1:1367971/1‑151 (MQ=255)
                        tttCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATaa                                                                                                                      >  2:8405483/1‑151 (MQ=255)
                          tCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGCACCGGGCCAACCATAAGt                                                                                                                    >  1:174998/1‑151 (MQ=255)
                          tCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGt                                                                                                                    <  1:3969459/151‑1 (MQ=255)
                              cccAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGc                                                                                                                <  2:11137256/151‑1 (MQ=255)
                              cccAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGc                                                                                                                >  2:8132388/1‑151 (MQ=255)
                                cAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCAc                                                                                                              <  1:10422715/151‑1 (MQ=255)
                                cAGCGAATGCCGCGCACTCCCGCCGACGCGGCGGATAACCGCGTCGCGGTCATTCGCCCTCCCGCCCCACGCCCTAGAGGGCCCAGGGCGAATACCGCCACAGGCCTTATCCGCCGAACCACTCCGAACCGTGCCAACCATAAGAcaccac                                                                                                              <  2:5554271/151‑1 (MQ=255)
                                 aGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACt                                                                                                             <  2:9131670/151‑1 (MQ=255)
                                  ccGAATGCCGCGCACTCCCGCCGACGCGGCGGAGAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACtt                                                                                                            <  1:10872321/150‑1 (MQ=255)
                                   cGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACttt                                                                                                           <  2:2313160/151‑1 (MQ=255)
                                   cGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACttt                                                                                                           <  1:10570701/151‑1 (MQ=255)
                                   cGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGGCAGCACttt                                                                                                           >  2:5225831/1‑151 (MQ=255)
                                    gAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTc                                                                                                          <  1:8589705/151‑1 (MQ=255)
                                    gAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCGTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGATCCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTc                                                                                                          >  2:10542181/1‑151 (MQ=255)
                                       tGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCaaa                                                                                                       <  2:9870387/151‑1 (MQ=255)
                                        gCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAg                                                                                                      >  1:587514/1‑151 (MQ=255)
                                           gcgcACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGccc                                                                                                   <  2:7568190/151‑1 (MQ=255)
                                               aCTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACg                                                                                               <  1:7513709/151‑1 (MQ=255)
                                               aCTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACg                                                                                               <  2:9175438/151‑1 (MQ=255)
                                                cTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGt                                                                                              <  2:7028959/151‑1 (MQ=255)
                                                cTCCCGCCGAGGCGGCGGATAACCGCGGTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTGATCCGCCGAAGCTCTCCGTACCGTGCCAACCATAAGGCTGCTCTGTCTACGCCCTACGt                                                                                              >  2:7827222/1‑151 (MQ=37)
                                                 tCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTg                                                                                             <  2:7984115/151‑1 (MQ=255)
                                                  cccGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGc                                                                                            >  1:5968934/1‑151 (MQ=255)
                                                     gccgAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCttt                                                                                         <  2:613599/151‑1 (MQ=255)
                                                      ccgAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCtttt                                                                                        <  2:4886882/151‑1 (MQ=255)
                                                          ggcggcGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTa                                                                                    <  1:7309781/151‑1 (MQ=255)
                                                          ggcggcGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTa                                                                                    <  2:362789/151‑1 (MQ=255)
                                                              gcggATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACaaa                                                                                <  1:6334864/151‑1 (MQ=255)
                                                               cggATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACaaaa                                                                               <  1:709050/151‑1 (MQ=255)
                                                                 gATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAaga                                                                             >  1:167252/1‑151 (MQ=255)
                                                                   tAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGc                                                                           <  2:1460475/151‑1 (MQ=255)
                                                                   tAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGc                                                                           >  2:3330874/1‑151 (MQ=255)
                                                                     aCCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCAt                                                                         >  1:2145834/1‑151 (MQ=255)
                                                                      ccGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATa                                                                        >  2:935917/1‑151 (MQ=255)
                                                                        gcgTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATaaa                                                                      <  1:10203808/151‑1 (MQ=255)
                                                                            tGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGATCCTCTCCGAACCGGGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAACATGc                                                                  <  1:6217151/151‑1 (MQ=255)
                                                                            tGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCAAAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGc                                                                  <  1:9484158/151‑1 (MQ=255)
                                                                             gCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCt                                                                 <  2:2958937/151‑1 (MQ=255)
                                                                              cGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTc                                                                <  2:2471923/151‑1 (MQ=255)
                                                                               ggTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCa                                                               <  1:227073/151‑1 (MQ=255)
                                                                                 ttATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAAt                                                             >  2:5707248/1‑151 (MQ=255)
                                                                                  tATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAAtt                                                            <  1:1953679/151‑1 (MQ=255)
                                                                                   aTTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTa                                                           <  1:5561812/151‑1 (MQ=255)
                                                                                   aTTCGCCCGACCGCACCACGACCTTGTGGGCGTAGGGCGACTACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTa                                                           <  2:10296212/151‑1 (MQ=255)
                                                                                    ttCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTAt                                                          <  1:1128252/151‑1 (MQ=255)
                                                                                    atCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTAt                                                          <  1:1520469/150‑1 (MQ=255)
                                                                                       gCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGCACCTCTCCGACCCGTGCCACCCCTCAGTCAGCACTTTCAAAGCCCTACGTGCATTTCCTACAAAAGAGCATAAAATGCTCAATTATggc                                                       <  1:1231585/151‑2 (MQ=255)
                                                                                       gCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATggg                                                       <  1:10061597/151‑1 (MQ=255)
                                                                                       gCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCACTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATggg                                                       >  2:2899523/1‑151 (MQ=255)
                                                                                        cccTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTCCGTGCTTTTCCTCCAAAAGAGCATAAAATGCTCAATGATGGGc                                                      >  2:5986125/1‑151 (MQ=255)
                                                                                        cccTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGc                                                      <  1:8361361/151‑1 (MQ=255)
                                                                                         ccTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCa                                                     >  1:9856891/1‑151 (MQ=255)
                                                                                         ccTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCa                                                     <  1:2549275/151‑1 (MQ=255)
                                                                                         ccTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCa                                                     >  1:4234106/1‑151 (MQ=255)
                                                                                          cTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCAc                                                    >  2:3926169/1‑151 (MQ=255)
                                                                                           tACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCAcc                                                   >  1:7673789/1‑151 (MQ=255)
                                                                                            aCCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCa                                                  <  1:9275144/151‑1 (MQ=255)
                                                                                            aCCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCa                                                  <  2:6823794/151‑1 (MQ=255)
                                                                                             ccGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAg                                                 >  2:1342118/1‑151 (MQ=255)
                                                                                              cGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGt                                                <  1:5321176/151‑1 (MQ=255)
                                                                                                caccacGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTc                                              >  2:8957858/1‑151 (MQ=255)
                                                                                                  ccacGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGAGCTTTTCCTACAAAAGAGCATAAACTGCTCAATTATGGGCACCAGTTcac                                            <  2:6671486/151‑1 (MQ=255)
                                                                                                   cacGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTcaca                                           >  2:3967730/1‑151 (MQ=255)
                                                                                                     cGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACt                                         <  1:4448524/151‑1 (MQ=255)
                                                                                                        ccTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGATCCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTgcg                                      <  2:2775403/151‑1 (MQ=255)
                                                                                                           tgtgGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGccc                                   >  1:10731097/1‑151 (MQ=255)
                                                                                                              gggCGTAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCAc                                <  1:9791981/151‑1 (MQ=255)
                                                                                                                   tAGGGCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGccc                           >  2:1741483/1‑151 (MQ=255)
                                                                                                                     gggCGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCACCCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGccccc                         <  2:6356569/151‑1 (MQ=255)
                                                                                                                        cGAATACCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACt                      <  1:2513705/151‑1 (MQ=255)
                                                                                                                             aCCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTtctc                 >  1:8698115/1‑151 (MQ=255)
                                                                                                                             aCCGCCACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTtctc                 >  1:4630849/1‑151 (MQ=255)
                                                                                                                                 ccACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAACAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCt             <  2:2661496/151‑1 (MQ=255)
                                                                                                                                 ccACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCt             <  2:4321517/151‑1 (MQ=255)
                                                                                                                                 ccACAGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCt             <  2:857440/151‑1 (MQ=255)
                                                                                                                                  cacaGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTGTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAATTGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTGTCTCGTCtt            >  1:3681864/1‑151 (MQ=255)
                                                                                                                                     aGGCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTCCGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCTTAGt         <  1:10799706/151‑1 (MQ=255)
                                                                                                                                       gCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAATGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCTTAGTGa       <  1:8859826/151‑1 (MQ=255)
                                                                                                                                       gCGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCTTAGTGa       <  1:3330650/151‑1 (MQ=255)
                                                                                                                                        cGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCTTAGTGAg      >  1:9146972/1‑151 (MQ=255)
                                                                                                                                        cGTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCTTAGTGAg      >  1:2986037/1‑151 (MQ=255)
                                                                                                                                         gTTATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCTTAGTGAGc     <  1:8044373/151‑1 (MQ=255)
                                                                                                                                          ttATCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCTTAGTGAGCa    >  2:8280998/1‑151 (MQ=255)
                                                                                                                                            aTCCGCCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCTTAGTGAGCAAc  >  2:1883013/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGCCGGCTGTGATGGCCGGCTTTTTTTCGCCCCAGCGAATGCCGCGCACTCCCGCCGAGGCGGCGGATAACCGCGTTGCGGTTATTCGCCCTACCGCACCACGACCTTGTGGGCGTAGGGCGAATACCGCCACAGGCGTTATCCGCGCGAACCTCTCCGAACCGTGCCAACCATAAGTCAGCACTTTCAAAGCCCTACGTGCTTTTCCTACAAAAGAGCATAAAATGCTCAATTATGGGCACCAGTTCACACTGCGCCCCACGGCCCCCACTTTCTCGTCTTAGTGAGCAAC  >  NC_008463/4700661‑4700952

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