Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	20,600	G→A	intergenic (‑92/+215)	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,600	0	G	A	90.9%	25.8 / ‑3.3	12	intergenic (‑92/+215)	insA1/rpsT	IS1 family protein InsA/30S ribosomal subunit protein S20
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (6/4);  total (8/4)
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.12e-01

TCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGT  >  NC_000913/20530‑20824                                                                       |                                                                                                                                                                                                                                               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)                                                          ccacattcttTGTAGGGG‑GTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTgcaagaacctgcaacttttgactatctggcaaccggcagtgtgttctctcacgcatcacaaaagcagcaggcataaaaaaacccgcttgcgcgggctttttcacaaagcttcagcaaattggcgattaagccagt  <  1:480818‑M1/247‑136 (MQ=255)                                                                       |                                                                                                                                                                                                                                               TCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGAGTAGATATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGT  >  NC_000913/20530‑20824

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