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.7 / 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 (13/11);  total (13/11)

AGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCC  >  NC_000913/1920091‑1920328
                                                                                                            |                                                                                                                                  
aGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGcc                                                                                                      >  1:888554/1‑139 (MQ=255)
        gAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                              <  2:888554/139‑1 (MQ=255)
                    tGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                        <  2:1083662/101‑1 (MQ=255)
                    tGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                        >  1:1083662/1‑101 (MQ=255)
                           aaTGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGcc                                                                           >  2:531576/1‑139 (MQ=255)
                                                    atGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                              <  2:698708/95‑1 (MQ=255)
                                                    atGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                              >  1:698708/1‑95 (MQ=255)
                                                    atGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGa                                                  <  1:531576/139‑1 (MQ=255)
                                                            tGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACg                                                                             <  1:508640/104‑1 (MQ=255)
                                                            tGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACg                                                                             >  2:508640/1‑104 (MQ=255)
                                                            tGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgc                                                                                                                          >  2:1329626/1‑55 (MQ=255)
                                                            tGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgc                                                                                                                          <  1:1329626/58‑4 (MQ=255)
                                                            tGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGtttttttt                                                                                                                            >  1:500561/1‑55 (MQ=255)
                                                            tGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGtttttttt                                                                                                                            <  2:500561/56‑2 (MQ=255)
                                                                      tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                              <  1:146054/77‑1 (MQ=255)
                                                                      tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                              >  2:146054/1‑77 (MQ=255)
                                                                        ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGAc                                                                              <  1:839826/91‑1 (MQ=255)
                                                                        ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGAc                                                                              >  2:839826/1‑91 (MQ=255)
                                                                        ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATg                                                >  2:119252/1‑121 (MQ=255)
                                                                        ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATg                                                <  1:119252/121‑1 (MQ=255)
                                                                                 gAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACg                                                                             <  1:694470/83‑1 (MQ=255)
                                                                                 gAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACg                                                                             >  2:694470/1‑83 (MQ=255)
                                                                                              aaaCTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCg                                                                                                           <  1:73382/40‑1 (MQ=255)
                                                                                              aaaCTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCg                                                                                                           >  2:73382/1‑40 (MQ=38)
                                                                                              aaaCTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTc        >  1:610/1‑139 (MQ=255)
                                                                                                    cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTccc                                                                                  >  1:720643/1‑59 (MQ=255)
                                                                                                    cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTccc                                                                                  <  2:720643/59‑1 (MQ=255)
                                                                                                    cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGcc  >  2:82654/1‑139 (MQ=255)
                                                                                                            |                                                                                                                                  
AGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCC  >  NC_000913/1920091‑1920328

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