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% 75.5 / NA 25*235* (TAG→TAAsfsAsugar fermentation stimulation protein A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/14);  total (11/14)

GAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTT  >  NC_000913/160648‑160888
                                                                                                                                      |                                                                                                          
gAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACa                                                                                                        <  2:1291278/139‑1 (MQ=255)
      tGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACa                                                                                                        >  2:974131/1‑133 (MQ=255)
      tGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACa                                                                                                        <  1:974131/133‑1 (MQ=255)
              cTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTTTGTGAAGAAAAGCGTATTTGCGCACGCGACCAGAAAGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCCGCAGt                                                                                          <  1:466854/139‑1 (MQ=255)
                       gCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTa                                                                                 <  2:497290/139‑1 (MQ=255)
                           aTGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAagag                                                                             >  1:1265999/1‑139 (MQ=255)
                           aTGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAagag                                                                             >  1:131394/1‑139 (MQ=255)
                             gCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGcc                                                                           >  1:559568/1‑139 (MQ=255)
                                                           acacTTGACAACTGTGTGAGGAAAAGTGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAACGTAACCGGCAGt                                                                                          <  1:1358615/94‑1 (MQ=255)
                                                           acacTTGACAACTGTGTGAGGAAAAGTGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAACGTAACCGGCAGt                                                                                          >  2:1358615/1‑94 (MQ=255)
                                                           acacTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTaa                                             <  2:590450/139‑1 (MQ=255)
                                                                         tgtgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACa                                                                                                        <  1:188079/66‑1 (MQ=255)
                                                                         tgtgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACa                                                                                                        >  2:188079/1‑66 (MQ=255)
                                                                           tgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGcc                                                                           <  1:1077193/93‑1 (MQ=255)
                                                                           tgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGcc                                                                           >  2:1077193/1‑93 (MQ=255)
                                                                           tgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCtt                                                                    <  2:1046243/100‑1 (MQ=255)
                                                                           tgAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCtt                                                                    >  1:1046243/1‑100 (MQ=255)
                                                                                   aGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagca                                                               <  1:628722/97‑1 (MQ=255)
                                                                                   aGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagca                                                               >  2:628722/1‑97 (MQ=255)
                                                                                   aGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGt                     >  2:1330101/1‑139 (MQ=255)
                                                                                            gcgcACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCt            <  2:559568/139‑1 (MQ=255)
                                                                                               cACGCGACCAGAATGTAAATTAACCAGTTACTTACTttatta                                                                                                          >  2:654984/1‑42 (MQ=255)
                                                                                                      ccAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGtt  <  2:131394/139‑1 (MQ=255)
                                                                                                      ccAGAATGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGtt  <  2:248496/139‑1 (MQ=255)
                                                                                                          aaTGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAccccc                           >  2:1090766/1‑110 (MQ=255)
                                                                                                          aaTGTAAATTAACCAGTTACTTACTTTATTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTAccccc                           <  1:1090766/110‑1 (MQ=255)
                                                                                                                                actttattACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGtt                    >  2:270253/1‑95 (MQ=255)
                                                                                                                                actttattACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGtt                    <  1:270253/95‑1 (MQ=255)
                                                                                                                                      |                                                                                                          
GAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAATTTCCGCTTTGTAAGCCAGAATTTCTACCCCCCTCTGTTGAGCTTCTGACAATAGTT  >  NC_000913/160648‑160888

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