Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 572,454 C→T C29C (TGC→TGT ylcH → uncharacterized protein, DLP12 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913572,4540CT100.0% 71.4 / NA 23C29C (TGC→TGTylcHuncharacterized protein, DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/13);  total (10/13)

CCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACC  >  NC_000913/572320‑572578
                                                                                                                                      |                                                                                                                            
ccGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAg                                                                                                                          <  2:209992/139‑1 (MQ=255)
                        cTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGa                                                                                                  <  2:714895/139‑1 (MQ=255)
                             gcgcTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                             <  2:377684/139‑1 (MQ=255)
                             gcgcTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                             <  2:314028/139‑1 (MQ=255)
                                   aaCGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTCCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAgggg                                                                                                                       >  1:345973/1‑107 (MQ=255)
                                   aaCGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTCCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAgggg                                                                                                                       <  2:345973/107‑1 (MQ=255)
                                                        aTACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCg                                                                  >  2:446429/1‑139 (MQ=255)
                                                                  aaCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                              <  1:594009/85‑1 (MQ=255)
                                                                  aaCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCCTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGaa                                                                                                              >  2:594009/1‑85 (MQ=255)
                                                                   aCTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGAACTTCACCGCTATCGGTCGGc                                                      <  1:392863/139‑1 (MQ=255)
                                                                        cttctCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGAt                                                  >  2:4045/1‑139 (MQ=255)
                                                                              ccTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGcc                                            <  1:4045/139‑1 (MQ=255)
                                                                                         gCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTggagg                                 >  1:81214/1‑139 (MQ=255)
                                                                                                     cTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACg                                                                                                <  1:200858/64‑1 (MQ=255)
                                                                                                     cTGCAATTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACg                                                                                                >  2:200858/1‑64 (MQ=255)
                                                                                                          aTTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                             <  2:743006/62‑1 (MQ=255)
                                                                                                          aTTCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa                                                                                             >  1:743006/1‑62 (MQ=255)
                                                                                                            tCTTTGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCAt              <  1:520913/139‑1 (MQ=255)
                                                                                                                tGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATg             <  1:225412/136‑1 (MQ=255)
                                                                                                                tGCATCTAGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATg             >  2:225412/1‑136 (MQ=255)
                                                                                                                       aGCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCACATGGTGCTTAAAc   >  1:8361/1‑139 (MQ=255)
                                                                                                                        gCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCACATGGTGCTTAAAcc  <  2:8361/139‑1 (MQ=255)
                                                                                                                        gCGTTTGACCTCTGTGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATATGGAGGAGGGAAAATGCTTTCGCATGGTGCTTAAAcc  >  2:533720/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                            
CCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACC  >  NC_000913/572320‑572578

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