Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 160,782 C→T *235* (TAG→TAA sfsA ← sugar fermentation stimulation protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913160,7820CT100.0% 81.9 / NA 26*235* (TAG→TAAsfsAsugar fermentation stimulation protein A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/10);  total (16/10)

GAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGGCGCG  >  NC_000913/160648‑160915
                                                                                                                                      |                                                                                                                                     
gAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACa                                                                                                                                   <  1:266626/139‑1 (MQ=255)
                       gCAGATGCTTTTCCGGATAGCAAGTATCTAAACGTAACACTTGACAACCGAGTGAGGAAAAGCGTATTTGCGCACGCGCCCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTa                                                                                                            >  2:604349/1‑139 (MQ=255)
                           aTGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAagag                                                                                                        >  1:826953/1‑139 (MQ=255)
                             gCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGcc                                                                                                      <  1:14501/139‑1 (MQ=255)
                                       aTAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcag                                                                                            <  1:534548/139‑1 (MQ=255)
                                                           acacTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTaa                                                                        >  2:687127/1‑139 (MQ=255)
                                                            cacTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAg                                                                       >  2:757830/1‑139 (MQ=255)
                                                                 gACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGa                                                                  <  1:757830/139‑1 (MQ=255)
                                                                    aaCTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTa                                                                                                            <  2:72043/94‑1 (MQ=255)
                                                                    aaCTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTa                                                                                                            >  1:72043/1‑94 (MQ=255)
                                                                    aaCTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGTAGAAATTTCCGCTTTGTAAGCCAGAAtt                                                               >  1:632742/1‑139 (MQ=255)
                                                                      cTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTc                                                             >  1:136490/1‑139 (MQ=255)
                                                                         tgtgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAc                                                          >  2:948191/1‑139 (MQ=255)
                                                                             aGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAccccc                                                      <  1:26253/139‑1 (MQ=255)
                                                                              ggAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAcccccc                                                     >  1:638302/1‑139 (MQ=255)
                                                                                aaaaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCctc                                                   >  1:23541/1‑139 (MQ=255)
                                                                                aaaaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCctc                                                   >  1:123339/1‑139 (MQ=255)
                                                                                    gCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGtt                                               <  1:948191/139‑1 (MQ=255)
                                                                                       tATTTGCGCACGCGACCAGAATGTAAATTAACCAGATACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAg                                            <  2:67237/139‑1 (MQ=255)
                                                                                                             gTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTAt                     <  2:1020993/139‑1 (MQ=255)
                                                                                                              tAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTAt                     >  2:523541/1‑139 (MQ=255)
                                                                                                                    aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCa               >  2:115928/1‑139 (MQ=255)
                                                                                                                    aaCCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCa               >  2:864444/1‑139 (MQ=255)
                                                                                                                      ccAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATc             >  1:1042873/1‑139 (MQ=255)
                                                                                                                          ttacttacTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGtt                                               <  1:676420/101‑1 (MQ=255)
                                                                                                                          ttacttacTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGtt                                               >  2:676420/1‑101 (MQ=255)
                                                                                                                                 ctttattACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGgcgcg  <  1:687127/139‑1 (MQ=255)
                                                                                                                                 ctttattACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGgcgcg  <  2:946360/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
GAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTTGCGCGTATTTCTCATCGATGTGGCGCG  >  NC_000913/160648‑160915

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: