Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,392,890 G→A *300* (TAG→TAA pgrR → murein peptide degradation regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,392,8900GA100.0% 102.7 / NA 33*300* (TAG→TAApgrRmurein peptide degradation regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/15);  total (18/15)

AATTCACACCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAGCACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAACGT  >  NC_000913/1392781‑1393023
                                                                                                             |                                                                                                                                     
aaTTCACACCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttc                                                                                                          <  2:64215/139‑1 (MQ=255)
      cacCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAacac                                                                                                                                    <  2:1360839/107‑1 (MQ=255)
      cacCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAacac                                                                                                                                    >  1:1360839/1‑107 (MQ=255)
                                     tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttc                                                                                                          <  1:1697208/102‑1 (MQ=255)
                                     tCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAAttc                                                                                                          >  2:1697208/1‑102 (MQ=255)
                                       cacaCCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGt                                                                   >  1:337391/1‑139 (MQ=255)
                                          acCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCAt                                                                >  1:1346635/1‑139 (MQ=255)
                                            cGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAg                                                                                      >  2:781072/1‑115 (MQ=255)
                                            cGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAg                                                                                      <  1:781072/115‑1 (MQ=255)
                                               cGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCt                                                                                                                            >  2:1664660/1‑74 (MQ=255)
                                               cGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCt                                                                                                                            >  2:1682570/1‑74 (MQ=255)
                                               cGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCt                                                                                                                            <  1:1664660/74‑1 (MQ=255)
                                               cGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCt                                                                                                                            <  1:1682570/74‑1 (MQ=255)
                                               cGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAt                                                                                                                     >  1:188623/1‑81 (MQ=255)
                                               cGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACAt                                                                                                                     <  2:188623/81‑1 (MQ=255)
                                                    tGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTAttct                                                      <  2:337391/139‑1 (MQ=255)
                                                         gCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCtt                                                                                                   <  2:1106385/89‑1 (MQ=255)
                                                         gCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCtt                                                                                                   >  1:1106385/1‑89 (MQ=255)
                                                              gCATACTCATTGTTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATg                                            >  2:1005662/1‑139 (MQ=255)
                                                                        tatttatAGATAGGCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATg                                  <  1:1005662/139‑1 (MQ=255)
                                                                                    ggCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATa                      <  2:1728000/139‑1 (MQ=255)
                                                                                     gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAg                                                                                      <  2:1707498/74‑1 (MQ=255)
                                                                                     gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAg                                                                                      <  1:612005/74‑1 (MQ=255)
                                                                                     gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAg                                                                                      >  2:612005/1‑74 (MQ=255)
                                                                                     gCTGAAGTATAAAGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAg                                                                                      >  1:1707498/1‑74 (MQ=255)
                                                                                               aaaGGTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACggg           >  1:143409/1‑139 (MQ=255)
                                                                                                   gTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTGTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTg       >  2:1225572/1‑139 (MQ=255)
                                                                                                   gTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGTTTTg       >  2:920390/1‑139 (MQ=255)
                                                                                                   gTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGtt         <  2:1106087/137‑1 (MQ=255)
                                                                                                   gTGCTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGtt         >  1:1106087/1‑137 (MQ=255)
                                                                                                      cTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGt                                                                                              >  2:1165371/1‑49 (MQ=255)
                                                                                                      cTGTTTAACACTACTTGCTGATACATTAATTTAATTCTTATCTTAACGt                                                                                              <  1:1165371/49‑1 (MQ=255)
                                                                                                       tGTTTAACACTACTTGCTGATACATTAATTTAACTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCAtttttt                                                           >  2:574229/1‑83 (MQ=255)
                                                                                                       tGTTTAACACTACTTGCTGATACATTAATTTAACTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCAtttttt                                                           <  1:574229/83‑1 (MQ=255)
                                                                                                        gTTTAACACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAACGt  >  2:398626/1‑139 (MQ=255)
                                                                                                             |                                                                                                                                     
AATTCACACCAGATTTACCCGGTTATCACCTGTACTATCCACACCGTCGACATGCTGGCTCGGCATTCTCATTATTTATAGATAGGCTGAAGTATAAAGGTGCTGTTTAGCACTACTTGCTGATACATTAATTTAATTCTTCTCTTAACGTATTCTCAGTTCCTTTCAACGTTTTGGTCATTTTTTATTCTTCGTACAATGGCGACAGATGCTGATTATGATACCGAAAACGGGTTTGAACGT  >  NC_000913/1392781‑1393023

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