Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,235,293	C→A	E195* (GAA→TAA)	rlmG ←	23S rRNA m(2)G1835 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,235,293	0	C	A	100.0%	79.0 / NA	24	E195* (GAA→TAA)	rlmG	23S rRNA m(2)G1835 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (14/10);  total (14/10)

CCATTACCGCAACCGAGATCGACAATCTCCCCTTCGAGATTCTCTGGCAGATGTTGCATAAAGAAGCGCGCGCCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGAACAGTTCCA  >  NC_000913/3235157‑3235421                                                                                                                                         |                                                                                                                                 ccATTACCGCAACCGAGATCGACAATCTCCCCTTCGAGATTCTCTGGCAGATGTTGCATAAAGAAGCGCGCGCCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACcttacag                                                                                                                               >  2:612773/1‑140 (MQ=255)   aTTACCGCAACCGAGATCGACAATCTCCCCTTCGAGATTCTCTGGCAGATGTTGCATAAAGAAGCGCGCGCCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACcttacagct                                                                                                                             >  1:547523/1‑140 (MQ=255)            aCCGAGATCGACAATCTCCCCTTCGAGATTCTCTGGCAGATGTTGCATAAAGAAGCGCGCGCCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTACAGCTTCCAGCTaa                                                                                                                    <  2:547523/140‑1 (MQ=255)                                cTTCGAGATTCTCTGGCAGATGTTGCATAAAGAAGCGCGCGCCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGcc                                                                                                <  1:255489/140‑1 (MQ=255)                                            ctGGCAGATGTTGCATAAAGAAGCGCGCGCCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCa                                                                                    >  1:356808/1‑140 (MQ=255)                                             tGGCAGATGTTGCATAAAGAAGCGCGCGCCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCAt                                                                                   <  1:612773/140‑1 (MQ=255)                                                   aTGTTGCATAAAGAAGCGCGCGCCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAg                                                                             <  1:376572/140‑1 (MQ=255)                                                              agaagCGCGCGCCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATc                                                                  >  2:219487/1‑140 (MQ=255)                                                                   cgcgcgcCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATc                                                                  <  1:364101/135‑1 (MQ=255)                                                                   cgcgcgcCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATc                                                                  >  2:364101/1‑135 (MQ=255)                                                                        gcCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCtt                                                        >  2:595149/1‑140 (MQ=255)                                                                                      cGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGc                                          <  2:530498/140‑1 (MQ=255)                                                                                        gTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCgt                                        >  1:110751/1‑140 (MQ=255)                                                                                                                 tgtgGATAGTCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACtttt               <  1:342512/140‑1 (MQ=255)                                                                                                                      aTAGTCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGaa          <  2:356808/140‑1 (MQ=255)                                                                                                                          tCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGaa          >  1:587067/1‑136 (MQ=255)                                                                                                                          tCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGaa          <  2:587067/136‑1 (MQ=255)                                                                                                                          tCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGAACAGt      >  2:645818/1‑140 (MQ=255)                                                                                                                          tCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGAACAGt      >  1:256845/1‑140 (MQ=255)                                                                                                                          tCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGAACAGt      >  2:527299/1‑140 (MQ=255)                                                                                                                          tCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGAACAGt      >  2:356615/1‑140 (MQ=255)                                                                                                                          tCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGAACAGt      >  1:345277/1‑140 (MQ=255)                                                                                                                          tCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGGGGTTGGGCCGAGCACTTTTTCGAACAGt      >  2:270851/1‑140 (MQ=255)                                                                                                                          tCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGCGGTGGTTGGGCCGAGCACTTTTTCGAACAGt      >  2:629594/1‑140 (MQ=255)                                                                                                                          tCCAGTCAGTACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGACCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGAACAGt      >  1:295905/1‑140 (MQ=255)                                                                                                                              gtcagtACCTTACAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGAACAGTTCCa  <  1:645818/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                 CCATTACCGCAACCGAGATCGACAATCTCCCCTTCGAGATTCTCTGGCAGATGTTGCATAAAGAAGCGCGCGCCGATATCAAGCCCGGTGCGGGAGAAGACATTCGCATGGTTGTGGATAGTCCAGTCAGTACCTTCCAGCTTCCAGCTAACGGTCTGCGGCGCATCGGCCAGCTGCGGTTCATTGAAAGTGCAATTAATCAGGCGCGCTTTTTTCCATGCCAGCGTGGTGGTGGTTGGGCCGAGCACTTTTTCGAACAGTTCCA  >  NC_000913/3235157‑3235421

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑