Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,409,911	T→C	L213L (TTG→CTG)	rlmB →	23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,409,911	0	T	C	100.0%	203.0 / NA	54	L213L (TTG→CTG)	rlmB	23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/21);  total (33/21)

CGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGA  >  NC_000913/4409790‑4410041                                                                                                                          |                                                                                                                                   cgGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATg                                                                                                                   >  2:77061/1‑139 (MQ=255)               aggcgGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTttc                                                                                                     <  1:168586/139‑1 (MQ=255)                 gcgGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTttctt                                                                                                   >  2:1072774/1‑139 (MQ=255)                     aTCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCt                                                                                               <  1:77061/139‑1 (MQ=255)                               ctATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCg                                                                                     >  2:367358/1‑139 (MQ=255)                               ctATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCg                                                                                     >  2:630183/1‑139 (MQ=255)                                tATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCgg                                                                                    >  1:159889/1‑139 (MQ=255)                                tATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCgg                                                                                    >  1:782964/1‑139 (MQ=255)                                 aTCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGt                                                                                   >  1:1194837/1‑139 (MQ=255)                                      agCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGttt                                                                                                      <  1:1063355/115‑1 (MQ=255)                                      agCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGttt                                                                                                      >  2:1063355/1‑115 (MQ=255)                                      agCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGa                                                                              >  2:634804/1‑139 (MQ=255)                                       gCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGAc                                                                             >  2:589249/1‑139 (MQ=255)                                        cAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGAcc                                                                            >  2:615367/1‑139 (MQ=255)                                             tGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAAt                                                                       <  1:428839/139‑1 (MQ=255)                                                ccGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTg                                                                    <  1:634804/139‑1 (MQ=255)                                                                        gTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAgcg                                            <  2:415406/139‑1 (MQ=255)                                                                           cGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTg                                                                                 >  2:640814/1‑99 (MQ=255)                                                                           cGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTg                                                                                 <  1:640814/99‑1 (MQ=255)                                                                           cGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCg                                                       <  2:1069086/125‑1 (MQ=255)                                                                           cGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCg                                                       >  1:1069086/1‑125 (MQ=255)                                                                           cGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAg                                         >  1:741473/1‑139 (MQ=255)                                                                           cGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAg                                         >  2:933451/1‑139 (MQ=255)                                                                           cGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATCGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTg                                                                                              >  1:503557/1‑86 (MQ=255)                                                                           cGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATCGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTg                                                                                              <  2:503557/86‑1 (MQ=255)                                                                               aggtgaaggtATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTaa                                     <  2:215670/139‑1 (MQ=255)                                                                                 gtgaaggtATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGAt                                                                                                                                <  2:1207420/46‑1 (MQ=255)                                                                                 gtgaaggtATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGAt                                                                                                                                >  1:1207420/1‑46 (MQ=255)                                                                                   gaaggtATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTtc                                 >  1:198466/1‑139 (MQ=255)                                                                                   gaaggtATGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTtc                                 >  2:423096/1‑139 (MQ=255)                                                                                         aTGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGaa                           <  2:435084/139‑1 (MQ=255)                                                                                         aTGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGaa                           >  1:435084/1‑139 (MQ=255)                                                                                         aTGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGaa                           <  1:615367/139‑1 (MQ=255)                                                                                          tGCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAAt                          >  2:507519/1‑139 (MQ=255)                                                                                           gCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAAt                                                                       <  1:807936/93‑1 (MQ=255)                                                                                           gCGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAAt                                                                       >  2:807936/1‑93 (MQ=255)                                                                                            cGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGTGCCAGCGCAGCTAATTtct                                >  2:1030856/1‑131 (MQ=255)                                                                                            cGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGTGCCAGCGCAGCTAATTtct                                <  1:1030856/131‑1 (MQ=255)                                                                                            cGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTtct                                <  2:340966/131‑1 (MQ=255)                                                                                            cGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTtct                                <  1:850658/131‑1 (MQ=255)                                                                                            cGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTtct                                <  2:501748/131‑1 (MQ=255)                                                                                            cGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTtct                                >  2:850658/1‑131 (MQ=255)                                                                                            cGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTtct                                >  1:501748/1‑131 (MQ=255)                                                                                            cGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTtct                                >  1:340966/1‑131 (MQ=255)                                                                                            cGTCGCCTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTa                        >  2:1195488/1‑139 (MQ=255)                                                                                                 ccTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTtct                                >  2:999302/1‑126 (MQ=255)                                                                                                 ccTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTtct                                <  1:999302/126‑1 (MQ=255)                                                                                                 ccTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGaaa                   >  1:667842/1‑139 (MQ=255)                                                                                                 ccTGACTCGTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTCTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGaaa                   >  1:1031786/1‑139 (MQ=255)                                                                                                         gTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTaa                                     >  2:585490/1‑113 (MQ=255)                                                                                                         gTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTaa                                     <  1:585490/113‑1 (MQ=255)                                                                                                         gTGAACATTGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATcc           >  1:509724/1‑139 (MQ=255)                                                                                                                 tGCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTgg   >  2:238956/1‑139 (MQ=255)                                                                                                                  gCGATGAGCTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGa  <  2:509724/139‑1 (MQ=255)                                                                                                                          |                                                                                                                                   CGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGA  >  NC_000913/4409790‑4410041

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: ‑