Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	20,590	T→C	intergenic (‑82/+225)	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,590	0	T	C	90.0%	21.8 / ‑3.6	10	intergenic (‑82/+225)	insA1/rpsT	IS1 family protein InsA/30S ribosomal subunit protein S20
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (6/3);  total (7/3)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

TCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGAT  >  NC_000913/20530‑20815                                                             |                                                                                                                                                                                                                                                tcatACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGCGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCAGGAAGAGTTCAGATAGTACAACGGCATGTCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAACAACCCGCTTgcg                                     >  2:454868/1‑251 (MQ=255)           aaaatagagtggtcttaagtttaatacctcaatgtgtatcacaatatccacatTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgctt                                        >  1:305898‑M1/25‑163 (MQ=255)                       tcttaagtttaatacctcaatgtgtatcaCAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttt                            >  1:232909‑M1/13‑151 (MQ=255)                           aagtttaatacctcaatgtgtATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttaac                        >  1:224553‑M1/9‑147 (MQ=255)                                  atacctcAATGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagctt                 >  1:169582‑M1/2‑140 (MQ=255)                                   taccgCAATGGGTATCACAATATACACATTCTTTGTCGGGG‑GTCTGGCGTTTGAGTCGATATTCTTGTTCAGATTGTATCAGCCGATGGTTCTACGATTCTTAAGCCACTAAGAGTTCAGATACTACCACGGCATGTCTGCAAGCACCTGCAACTTTTGACTATCTGGCCACCGGCAGTGTGTTCTCTACCGCATCACAAAAGCATCCGGCCCAAAAAATACCGCTTTTGCAGGCTTTTTCACAAAACTTc                >  2:31573/1‑251 (MQ=255)                                        cAATGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTCCAGAATGTATCAGCCGATGGTTCTACTATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagca            <  1:156224‑M1/250‑117 (MQ=255)                                         aaTGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaa           >  1:348856‑M1/1‑133 (MQ=255)                                         aaTGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgt                                                                                    >  1:147465‑M1/1‑133 (MQ=255)                                         aaTGTGTATCACAATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgt                                                                                    <  2:147465‑M1/177‑45 (MQ=255)                                                  cacaATATCCACATTCTTTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGGTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgat  <  2:232909‑M1/250‑127 (MQ=255)                                                             |                                                                                                                                                                                                                                                TCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGAT  >  NC_000913/20530‑20815

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