Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,327,342	G→A	A108V (GCA→GTA)	rlmE ←	23S rRNA 2'‑O‑ribose U2552 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,327,342	0	G	A	100.0%	109.2 / NA	34	A108V (GCA→GTA)	rlmE	23S rRNA 2'‑O‑ribose U2552 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/18);  total (16/18)

ATTTCTAGCGCCAGTTCCACCAGATACATGGCACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTGCTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACC  >  NC_000913/3327216‑3327478                                                                                                                               |                                                                                                                                         aTTTCTAGCGCCAGTTCCACCAGATACATGGCACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAg                                                                                                                             >  2:616018/1‑140 (MQ=255)   ttCTAGCGCCAGTTCCACCAGATACATGGCACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGtt                                                                                                                           >  1:654445/1‑140 (MQ=255)           ccAGTTCCACCAGATACATGGCACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGa                                                                                                                   >  2:521284/1‑140 (MQ=255)             aGTTCCACCAGATACATGGCACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGaaa                                                                                                                 >  1:126509/1‑140 (MQ=255)               ttCCACCAGATACATGGCACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAAt                                                                                                               >  2:505032/1‑140 (MQ=255)                    ccaGATACATGGCACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGccc                                                                                                          <  1:153438/140‑1 (MQ=255)                     caGATACATGGCACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCt                                                                                                         <  2:581769/140‑1 (MQ=255)                       gATACATGGCACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCAGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTga                                                                                                       >  1:620550/1‑140 (MQ=255)                         tACATGGCACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTgaag                                                                                                     >  2:473774/1‑140 (MQ=255)                          aCATGGCACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTgaaga                                                                                                    <  1:616018/140‑1 (MQ=255)                               agacaggggATATCCTCCGCCGATGATCCTCTCATGTATGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTc                                                                                               <  1:169987/135‑1 (MQ=255)                                cACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTcc                                                                                              <  2:77778/140‑1 (MQ=255)                                 aCGGGGGATATCCACTGCCGGTCTTGCGCCCATGTTTTGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCca                                                                                             <  1:196936/140‑1 (MQ=255)                                      ggATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCATTTCATCACGAAAATCGCCCTGCAGAAAGTCCACACCa                                                                                        >  1:257713/1‑140 (MQ=255)                                       gATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCaa                                                                                       <  2:55797/140‑1 (MQ=255)                                          atCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGa                                                                                    >  2:35897/1‑140 (MQ=255)                                             cACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAg                                                                                 <  2:639671/140‑1 (MQ=255)                                                cgccgGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCGGTTCATCAGGAAGATCGCCCTGAAGGAAGTCCACACCAACGactgtct                                                                              >  2:169987/1‑134 (MQ=255)                                                  ccgGTGTTCCGCTCACGTTTGGAGCGATAGCGGACATTACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATcc                                                                            <  1:316511/140‑1 (MQ=255)                                                    ggTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCAt                                                                          <  1:634026/140‑1 (MQ=255)                                                    ggTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCAt                                                                          <  1:391857/140‑1 (MQ=255)                                                                          ccATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACGAACGATACGATTCACAGGTAACAGATCGCAAGCgatg                                                    >  1:233755/1‑140 (MQ=255)                                                                             tatCGGACATGACATCCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGatgatg                                                 >  2:333438/1‑140 (MQ=255)                                                                                                    ttGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGt                          >  1:150020/1‑140 (MQ=255)                                                                                                    ttGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGt                          >  1:54082/1‑140 (MQ=255)                                                                                                     tGCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTg                         <  1:548288/140‑1 (MQ=255)                                                                                                      gCTGTCGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGa                        >  1:564362/1‑140 (MQ=255)                                                                                                       cTGTCGCCAACGCGCTCCAGCAATACTTTCATCACCAGTCCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGAc                       <  2:257713/140‑1 (MQ=255)                                                                                                           cGCCAACGCGCTCCAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACcaca                   <  2:576870/140‑1 (MQ=255)                                                                                                                       ccAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGaccaa       <  2:126509/140‑1 (MQ=255)                                                                                                                       ccAGCAGTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGaccaa       <  1:505032/140‑1 (MQ=255)                                                                                                                            agTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAaccacc  <  2:620550/140‑1 (MQ=255)                                                                                                                            agTACTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAaccacc  >  1:572988/1‑140 (MQ=255)                                                                                                                            agTACTTTCATCACCAGTTCATCACGAAAATAGCCTTGAGGAAAGTCCTCACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAaccacc  <  2:233755/140‑1 (MQ=255)                                                                                                                               |                                                                                                                                         ATTTCTAGCGCCAGTTCCACCAGATACATGGCACGGGGGATATCCACCGCCGGTGTTCCGCTCATGTTTGGTGCCATATCGGACATGACAACCTGGACTTTGCTGTCGCCAACGCGCTCCAGCAGTGCTTTCATCACCAGTTCATCACGAAAATCGCCCTGAAGAAAGTCCACACCAACGATAGGATCCATAGGTAAAAGATCGCAAGCGATGATGCGGCCTTTGCCGCCAATTTGGGTGACCACATATTGTGACCAACCACC  >  NC_000913/3327216‑3327478

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