Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	226,641	(G)5→6	noncoding (883/2904 nt)	rrlH →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	226,636	1	.	G	100.0%	55.8 / NA	27	noncoding (878/2904 nt)	rrlH	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (15/12);  total (15/12)

CTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAA‑GGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGAAGAGGG  >  NC_000913/226498‑226736                                                                                                                                            |                                                                                                     cTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAA‑g                                                                                                     <  1:149361/140‑1 (MQ=32)          aaaaaTTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAA‑ggggggtcat                                                                                            >  2:372873/1‑135 (MQ=17)                        gACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAAc                                                                               <  2:10995/140‑1 (MQ=17)                             gtgGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTTGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGAt                                                                          <  1:441050/140‑1 (MQ=12)                             gtgGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTTGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGAt                                                                          >  2:441050/1‑140 (MQ=12)                                    ggggTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACt                                                                   >  2:54440/1‑140 (MQ=17)                                     gggTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTg                                                                  >  1:275769/1‑140 (MQ=17)                                      ggTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGc                                                                 <  1:185334/140‑1 (MQ=17)                                           aaGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATa                                                            <  1:496681/140‑1 (MQ=17)                                                      aaCCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGt                                                 <  1:323809/140‑1 (MQ=17)                                                            gagaTAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCAc                                           >  1:341223/1‑140 (MQ=17)                                                            gagaTAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCAc                                           <  2:495989/140‑1 (MQ=17)                                                                         ctcCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACAcgg                              >  1:244744/1‑140 (MQ=18)                                                                          tcCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACAcggc                             <  1:54440/140‑1 (MQ=17)                                                                          tcCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACAcggc                             <  2:341223/140‑1 (MQ=17)                                                                             ccGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCggg                          <  2:461317/140‑1 (MQ=17)                                                                                aaaGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGc                       >  2:16496/1‑140 (MQ=18)                                                                                      aTTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAAcg                  >  2:263174/1‑139 (MQ=18)                                                                                      aTTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAAcg                  <  1:263174/139‑1 (MQ=17)                                                                                          aGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAAcgtccgt             <  2:244744/140‑1 (MQ=17)                                                                                              aGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTg         >  2:277428/1‑140 (MQ=18)                                                                                              aGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTg         <  1:16496/140‑1 (MQ=17)                                                                                               gcgcCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGa        >  2:510690/1‑140 (MQ=17)                                                                                               gcgcCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGGGCTATCGTCCGGCGTGa        >  2:406752/1‑140 (MQ=11)                                                                                                cgcCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGTTAACGTCCGTCGTGaa       >  2:407347/1‑140 (MQ=16)                                                                                                 gcCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGAag      >  2:395878/1‑140 (MQ=18)                                                                                                 gcCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGAag      >  2:421640/1‑140 (MQ=17)                                                                                                 gcCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGAag      >  1:462210/1‑140 (MQ=18)                                                                                                     cGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGAAGAggg  >  1:253053/1‑140 (MQ=18)                                                                                                                                            |                                                                                                     CTAATGTTGAAAAATTAGCGGATGACTTGTGGCTGGGGGTGAAAGGCCAATCAAACCGGGAGATAGCTGGTTCTCCCCGAAAGCTATTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAA‑GGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGAAGAGGG  >  NC_000913/226498‑226736

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