Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,731,034	G→A	noncoding (124/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,731,034	0	G	A	100.0%	63.2 / NA	25	noncoding (124/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/10);  total (15/10)

GGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAG  >  NC_000913/2730823‑2731264                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       ggACCGTGTCTCAATTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTCCTCACCCGTCCGCCACTCGTCAGCGAagcagc                                                                                                                                                                                                  >  2:334981/1‑250 (MQ=255) ggACCGTGTCCCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGACGCCTTGGTGAGCCGTTCCCCCACCAACCAGCTAATCCCGTATGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGAAGTTATCCCCCTCCATCAGGCCGCTTCCCCAACATTACTCCCCCGTCCGCCCCTCGTCACCGAcgcagcc                                                                                                                                                                                                 >  2:209661/1‑250 (MQ=11)                                             aGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTgcc                                                                                                                                                     >  2:317022/1‑251 (MQ=38)                                                          cgCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAAt                                                                                                                                        >  1:46550/1‑251 (MQ=38)                                                                gttacgacGTTACCCCACCAACAAGCTAATCTCATCTGGGCACATCCGCTGGCAAGCGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACCTTACTCACCTGTCCGCCACTCGACATCGACGCAGCCAGCTGAATCCTGCTACCTTTCATCATGTATTTGTTAGTCCAGCCGACAGCGCTCAATCtaatc                                                                                                                                  >  2:231069/8‑247 (MQ=2)                                                                    gccatttACCCCGACAACAAGCTAAGCCCATTTGGGCACAGCCGGTGGCAAGAGGCCAGACGGTCCCCCTCTTTGGTCTGGCGATATTATGCGGTATTAGCTCCCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACGTTTCTCACCCGTCCGCCACTCGGCAGCGAAGCAGCAAGCGGCTTCCGGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCAtt                                                                                                                              <  2:274844/246‑2 (MQ=11)                                                                            ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCCAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCACGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGCGCCATGATCAAAct                                                                                                                      >  1:371598/1‑251 (MQ=255)                                                                                 ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTc                                                                                                                  >  1:155103/1‑250 (MQ=39)                                                                                   aacaaGCTAATCCCATCTGGGCTCATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAAt                                                                                                               <  1:300704/251‑1 (MQ=255)                                                                                       agccatgcccacttggccacagccaaTGGAAAGAGGCCCGAAGGCCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCTTTCCCAGTAGTTCTCCCCCTCCATCAGGCAGCTTCCCAAACATTTCTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGGGCGACTTGAATGTGTTAGGCCTTCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTaa                                                                                                           <  2:107507/226‑1 (MQ=11)                                                                                       aGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTaa                                                                                                           <  1:98309/251‑1 (MQ=38)                                                                                             ccccATCTGGGCATATCCGATGACAAGAGGCCCGAAGGTCACCCTCTTTGGTCTTGCGACATTAGGCGGTATTAGCTACCGTTTCCATTAGTTATCCCCCTCCACCAGGCAGCTTCCCAAACATTGCGCGCCCGTCCGCCACTCGTCAGCGAAGCAGCATGCTGCTTCCTGTTACCGTTTGACTTGCATGTGTTAGTCCTGCCGGCAGCTTGCAATCTGTGCCATGATCAAACTCTTCAATTTAAAAGttt                                                                                                     <  1:57888/250‑1 (MQ=14)                                                                                                   atGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTGGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTACCCCCCTCCATCAGGCAGCTTCCAAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGc                                                                                                <  1:209661/249‑1 (MQ=255)                                                                                                              cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCg                                                                                                                                                                           <  2:136929/164‑1 (MQ=255)                                                                                                              cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCg                                                                                                                                                                           >  1:136929/1‑164 (MQ=255)                                                                                                                aTGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAAc                                                                                  >  1:185133/1‑251 (MQ=37)                                                                                                                                        ctctTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCAct                                                          >  2:90256/1‑251 (MQ=37)                                                                                                                                                              gcggcatttacTACCGTTTCCAGTCGTTATCCACCTCCATCCGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGCAGCAGCAAGCTCCTTCATGTTACCGTTCGACTTGCACGTGTTAGGCCTGCCGCCAGCGTCCAATCTGAGCCATGATCAAACTCTTCAATTTCAAAGTTTGATGCTCAAAGAATTACACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTAtacctttt                                    >  2:278062/11‑244 (MQ=11)                                                                                                                                                                             gTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGtt                     >  2:31137/1‑251 (MQ=21)                                                                                                                                                                                 ccattagtTATCCCCCTCCATCAGGCAGCTTCCCAAACATTGCTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAAGTACTTTTGCGCTCTTGAGATTTGGTATTCATTTTTCGTCTTGCGACGTTAAGa                 <  1:221146/247‑1 (MQ=11)                                                                                                                                                                                 ccAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACATTAAGa                 <  1:125402/251‑1 (MQ=21)                                                                                                                                                                                   agtagtTATCACCCTCCATCAGGCAGCTTCCCAAAAATTACTCACCCGTCCGCCACTCGTCAACGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAAt               >  2:263513/1‑251 (MQ=255)                                                                                                                                                                                     tagtTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATcc             >  1:104228/1‑251 (MQ=21)                                                                                                                                                                                              cccTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGCGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCg    >  2:5028/1‑251 (MQ=21)                                                                                                                                                                                                ccccATCCGGCAGCTTCCCAACCCTCCCTCCCCCGTCCGCCACGCGTCAGCGAAGCAGCAAGCGGCTTCCTGTTCCCGTTCGACTTGCATGTGGTAGGCCTGCCGCCAGCGTTCAAGCTGAGCCATTATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTTCGTGTTCACTCTTGCGACTTGGTATTCAGTTTGCGTCTTGCGCAGTTAAGAATCCGTATCTTCGAg  <  2:322197/249‑1 (MQ=11)                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       GGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAG  >  NC_000913/2730823‑2731264

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