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% 133.2 / NA 42intergenic (+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 (22/20);  total (22/20)

GGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920068‑1920330
                                                                                                                                   |                                                                                                                                    
ggAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGCACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGtttttttt                                                                                                                              <  1:1083072/139‑2 (MQ=255)
ggAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGtttttttt                                                                                                                              <  2:303604/139‑2 (MQ=255)
             ttCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          >  1:1476067/1‑131 (MQ=255)
             ttCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          <  2:1476067/131‑1 (MQ=255)
                cTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCt                                                                                                               <  1:500348/139‑1 (MQ=255)
                     gCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTgg                                                                                                          >  2:1110003/1‑139 (MQ=255)
                             gTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacac                                                                                                  >  2:1361063/1‑139 (MQ=255)
                                aaCCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCTGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCACCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACAcc                                                                                               >  2:287633/1‑139 (MQ=255)
                                     aaaGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTa                                                                                          <  1:1240791/139‑1 (MQ=255)
                                     aaaGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgcca                                                                                                                          <  1:250913/106‑6 (MQ=255)
                                     aaaGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgcca                                                                                                                          >  2:250913/1‑101 (MQ=255)
                                               gTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGAc                                                                                <  1:1361063/139‑1 (MQ=255)
                                                            tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          >  2:779233/1‑84 (MQ=255)
                                                            tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          <  1:779233/84‑1 (MQ=255)
                                                            tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGacac                                                                   >  2:589090/1‑139 (MQ=255)
                                                            tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGacac                                                                   >  2:1420019/1‑139 (MQ=255)
                                                                 ttCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCg                                                                                                             <  2:705657/92‑1 (MQ=255)
                                                                 ttCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCg                                                                                                             >  1:705657/1‑92 (MQ=255)
                                                                             gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          >  2:709288/1‑67 (MQ=255)
                                                                             gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          <  1:709288/67‑1 (MQ=255)
                                                                             gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACa                                                                                                                   >  2:1815413/1‑74 (MQ=255)
                                                                             gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACa                                                                                                                   <  1:1815413/74‑1 (MQ=255)
                                                                              cccACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCg                                                                                                             >  2:942833/1‑79 (MQ=255)
                                                                              cccACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCg                                                                                                             <  1:942833/79‑1 (MQ=255)
                                                                                cACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgc                                                                                                                            <  1:903336/61‑4 (MQ=255)
                                                                                cACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgc                                                                                                                            >  2:903336/1‑58 (MQ=255)
                                                                                   tGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgc                                                                                                                            <  1:666019/58‑4 (MQ=255)
                                                                                   tGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgc                                                                                                                            >  2:666019/1‑55 (MQ=255)
                                                                                    gCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCg                                                     >  1:550691/1‑129 (MQ=255)
                                                                                    gCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCg                                                     <  2:550691/129‑1 (MQ=255)
                                                                                      tCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCACCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACg                                         <  1:287633/139‑1 (MQ=255)
                                                                                        cGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAaca                                                                                                   >  1:579722/1‑79 (MQ=255)
                                                                                        cGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAaca                                                                                                   <  2:579722/79‑1 (MQ=255)
                                                                                             tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACg                                                                               >  2:326014/1‑94 (MQ=255)
                                                                                             tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACg                                                                               <  1:326014/94‑1 (MQ=255)
                                                                                             tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGcc                                                                             <  2:1829260/96‑1 (MQ=255)
                                                                                             tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGcc                                                                             >  1:1829260/1‑96 (MQ=255)
                                                                                               ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgc                                                                                                                            >  2:362858/1‑43 (MQ=255)
                                                                                               ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgc                                                                                                                            <  1:362858/46‑4 (MQ=255)
                                                                                                  gcgcCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          >  2:398125/1‑46 (MQ=255)
                                                                                                  gcgcCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          <  1:398125/46‑1 (MQ=39)
                                                                                                         aGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTc                      >  2:388288/1‑139 (MQ=255)
                                                                                                               tATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTgccg                <  1:388288/139‑1 (MQ=255)
                                                                                                                     aaaCTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGa                                           <  2:427922/106‑1 (MQ=255)
                                                                                                                     aaaCTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGa                                           >  1:427922/1‑106 (MQ=255)
                                                                                                                      aaCTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCa         >  2:1298133/1‑139 (MQ=255)
                                                                                                                       aCTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTgccgc               >  2:608090/1‑132 (MQ=255)
                                                                                                                       aCTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTgccgc               <  1:608090/132‑1 (MQ=255)
                                                                                                                           cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacaca                                                                                                 <  2:510301/46‑1 (MQ=255)
                                                                                                                           cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacaca                                                                                                 >  1:510301/1‑46 (MQ=39)
                                                                                                                            gcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCTACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGccc   >  1:1099493/1‑139 (MQ=255)
                                                                                                                             cgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCg  <  2:1099493/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
GGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920068‑1920330

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: