Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,920,205 (T)7→8 intergenic (+62/‑18) yebT → / → rsmF MCE domain protein/16S rRNA m(5)C1407 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,920,1981.T100.0% 74.4 / NA 24intergenic (+55/‑25)yebT/rsmFMCE domain protein/16S rRNA m(5)C1407 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (15/9);  total (15/9)

GGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑TTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920068‑1920330
                                                                                                                                     |                                                                                                                                     
ggAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑tttttttt                                                                                                                              <  2:1035753/139‑2 (MQ=255)
                cTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCt                                                                                                               >  2:80952/1‑139 (MQ=255)
                  gTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCg                                                                                                             <  1:80952/139‑1 (MQ=255)
                       aGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGccgtccggg                                                                                                        >  2:750950/1‑131 (MQ=255)
                                    gAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTTGCCAGTGGTAc                                                                                                                    <  2:588557/115‑1 (MQ=255)
                                    gAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTTGCCAGTGGTAc                                                                                                                    >  1:588557/1‑115 (MQ=255)
                                                             gcgcTTCCCAAATAATGCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGc                                                                              >  1:512485/1‑127 (MQ=255)
                                                             gcgcTTCCCAAATAATGCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGc                                                                              <  2:512485/127‑1 (MQ=255)
                                                                      aaaTAATGCCCACTGCTCCGGCGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATgcgc                                                           >  1:634257/1‑139 (MQ=255)
                                                                           atGCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATgcgc                                                           <  1:126140/132‑1 (MQ=255)
                                                                           atGCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATgcgc                                                           >  2:126140/1‑132 (MQ=255)
                                                                             gCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑ttttttttg                                                                                                                             >  2:233129/1‑61 (MQ=255)
                                                                             gCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑ttttttttg                                                                                                                             <  1:233129/63‑3 (MQ=255)
                                                                                   tGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACAt                                                                                                                  >  2:1039031/1‑69 (MQ=255)
                                                                                   tGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACAt                                                                                                                  <  1:1039031/69‑1 (MQ=255)
                                                                                               tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCg                                                                                   >  2:382566/1‑90 (MQ=255)
                                                                                               tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCg                                                                                   <  1:382566/90‑1 (MQ=255)
                                                                                                                             cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGt                                                                                                            <  1:864576/35‑1 (MQ=37)
                                                                                                                             cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGt                                                                                                            >  2:864576/1‑35 (MQ=37)
                                                                                                                             cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCtt                                              <  1:827752/97‑1 (MQ=255)
                                                                                                                             cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCtt                                              >  2:827752/1‑97 (MQ=255)
                                                                                                                             cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGcc    >  1:910522/1‑139 (MQ=255)
                                                                                                                             cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGcc    >  2:920026/1‑139 (MQ=255)
                                                                                                                              gcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGa                                           <  1:612357/99‑1 (MQ=255)
                                                                                                                              gcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGa                                           >  2:612357/1‑99 (MQ=255)
                                                                                                                              gcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGccc   >  2:513012/1‑139 (MQ=255)
                                                                                                                               cgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCg  >  1:829145/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                     
GGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCG‑‑GCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑TTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920068‑1920330

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