Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	20,605	Δ1 bp	intergenic (‑97/+210)	insA1 ← / ← rpsT	IS1 family protein InsA/30S ribosomal subunit protein S20

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	20,605	0	A	.	100.0%	39.8 / NA	10	intergenic (‑97/+210)	insA1/rpsT	IS1 family protein InsA/30S ribosomal subunit protein S20
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (8/2);  total (8/2)

TCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGC  >  NC_000913/20542‑20838                                                                |                                                                                                                                                                                                                                          aatagagtggtcttaagtttaatacctcaatgtgtatcacaatatccacATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgc                                               >  1:61710‑M1/23‑161 (MQ=255)   tagagtggtcttaagtttaatacctcactgtgtatcacaatatccACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgc                                             >  2:274043‑M1/21‑159 (MQ=255)                 gtttaatacctcaatgtGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGCGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcaca                                >  1:43009‑M1/7‑145 (MQ=255)                  tttaatacctcaatgTGTATCTCAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGATTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgcctatctggca                                                                                                          <  2:263753‑M1/170‑32 (MQ=255)                  tttaatacctcaatgTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggca                                                                                                          >  1:263753‑M1/6‑144 (MQ=255)                                 tgtATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTTCGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTCAGCCACGAAGAATTCa                                                                                                                                                             >  2:297790/1‑109 (MQ=255)                                 tgtATCACAATATCAACAGTCTTGGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCa                                                                                                                                                             <  1:297790/109‑1 (MQ=255)                                        caATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgatta         >  2:439508‑M1/1‑122 (MQ=255)                                          atatCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagc      >  1:77048‑M1/1‑120 (MQ=255)                                              ccacattcttTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagt  >  1:26828‑M1/5‑116 (MQ=255)                                                                |                                                                                                                                                                                                                                          TCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGC  >  NC_000913/20542‑20838

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