Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,622,797	Δ1 bp	intergenic (+56/+20)	mgtS → / ← mgrR	small protein MgtS/small regulatory RNA MgrR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,622,794	0	G	.	100.0%	41.6 / NA	11	intergenic (+53/+23)	mgtS/mgrR	small protein MgtS/small regulatory RNA MgrR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (7/4);  total (7/4)

TTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCT  >  NC_000913/1622717‑1622850                                                                              |                                                         tttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAgggc                                                       >  2:268547/1‑80 (MQ=255) tttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAggg                                                        <  1:268522/80‑1 (MQ=255)  ttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAgggcc                                                      >  1:268549/1‑79 (MQ=255)        cacaAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAA‑GGGCCTGAttt                                               >  1:268542/1‑81 (MQ=255)         acaAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAA‑GGGCCTGAtttt                                              >  1:268544/1‑81 (MQ=255)            aaTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAA‑GGGCCTGATTTTTAt                                           >  1:268543/1‑81 (MQ=255)                     aCTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAA‑GGGCCTGATTTTTATGACGGCGa                                   >  1:268545/1‑80 (MQ=255)                             aCGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAA‑GGGCCTGATTTTTATGACGGCGAAAAAAAAcc                          <  1:268523/81‑1 (MQ=255)                                    tAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAA‑GGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTa                   <  1:268524/81‑1 (MQ=255)                                            cACCTAACCGGCCCCTTGTTACAGTTGTGTACAA‑GGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcg           >  2:268552/1‑81 (MQ=255)                                                 aaCCGGCCCCTTGTTACAGTTGTGTACAA‑GGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGaa      >  2:268551/1‑81 (MQ=255)                                                  aCCGGCCCCTTGTTACAGTTGTGTACAA‑GGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGa       <  2:268525/79‑1 (MQ=255)                                                   ccGGCCCCTTGTTACAGTTGTGTACAA‑GGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATg    <  2:268526/81‑1 (MQ=255)                                                     ggCCCCTTGTTACAGTTGTGTACAA‑GGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCt  >  1:268550/1‑81 (MQ=255)                                                                              |                                                         TTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCT  >  NC_000913/1622717‑1622850

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