Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,029,884	C→T	R29H (CGC→CAC)	rlmI ←	23S rRNA m(5)C1962 methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,029,884	0	C	T	100.0%	107.1 / NA	28	R29H (CGC→CAC)	rlmI	23S rRNA m(5)C1962 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/14);  total (14/14)

AAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGCAAAGAAATGGGCGTCATTGTCCGGAGTTTTCCATGGATTTG  >  NC_000913/1029737‑1030031                                                                                                                                                    |                                                                                                                                                    aaaGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTggg                                                                                                                                                  >  2:291885/1‑151 (MQ=255)    gCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCaa                                                                                                                                               <  2:610817/151‑1 (MQ=255)         gTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGcc                                                                                                                                          <  1:410403/151‑1 (MQ=255)           cGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGcccc                                                                                                                                        <  1:339247/151‑1 (MQ=255)                 aCTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGaaa                                                                                                                                   >  2:56393/1‑150 (MQ=255)                     aGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGAcc                                                                                                                              >  2:448633/1‑151 (MQ=255)                                        aCGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGt                                                                                                            >  1:590652/1‑150 (MQ=255)                                           cgcgcCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATg                                                                                                        <  2:21566/151‑1 (MQ=255)                                             cgcCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAt                                                                                                      >  1:383347/1‑151 (MQ=255)                                               cccGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCCTTTTCCCTGATGCTCAACAATATCGATGGTTTCACCGCGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGAt                                                                                                      <  2:105591/149‑1 (MQ=255)                                                                  cGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAac                                                                                 >  1:440370/1‑151 (MQ=255)                                                                  cGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAac                                                                                 >  1:534015/1‑151 (MQ=255)                                                                                  gCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCa                                                                 <  2:442344/151‑1 (MQ=255)                                                                                           tttCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTAc                                                         <  1:392403/150‑1 (MQ=255)                                                                                                 tgatgaTCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTg                                                  >  2:584145/1‑151 (MQ=255)                                                                                                  gatgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCt                                                   <  2:25837/149‑1 (MQ=255)                                                                                                  gatgatCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTg                                                  <  2:84014/150‑1 (MQ=255)                                                                                                        caacaaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGc                                           >  2:442512/1‑151 (MQ=255)                                                                                                           caaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGCaaa                                        <  1:148606/151‑1 (MQ=255)                                                                                                           caaTATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGCaa                                         >  2:446116/1‑150 (MQ=255)                                                                                                                    tGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGCAAAGAAATGGGc                               >  2:417842/1‑151 (MQ=255)                                                                                                                    tGGTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGCAAAGAAATGGGc                               >  1:497837/1‑151 (MQ=255)                                                                                                                     ggTTTCACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGCAAAGAAATGGGCg                              <  2:47956/151‑1 (MQ=255)                                                                                                                          cACCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGCAAAGAAATGGGCGTCAtt                         >  2:486789/1‑151 (MQ=255)                                                                                                                           aCCGAGGCTGGCTTTACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGCAAAGAAATGGGCGTCATTg                        >  1:211402/1‑151 (MQ=255)                                                                                                                                         tACCTTCCATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGCAAAGAAATGGGCGTCATTGTCCGGAGTTTTcc           <  1:388820/150‑1 (MQ=255)                                                                                                                                               ccATGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGCAAAGAAATGGGCGTCATTGTCCGGAGTTTTCCATGGAtt    <  2:497827/151‑1 (MQ=255)                                                                                                                                                 aTGTGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGCAAAGAAATGGGCGTCATTGTCCGGAGTTTTCCATGGATTTg  <  2:38857/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCTGGCTTTACCTTCCATGCGGGCAACGGCCCCGGAAAAGACCCACGGATGGCGACGAAGTAATGATTTTTCGCGCCCTTTGGCTAACACTAAACGTACACTCATAATTTACTTTTCTGTCGATGCAAAGAAATGGGCGTCATTGTCCGGAGTTTTCCATGGATTTG  >  NC_000913/1029737‑1030031

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