Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	16,981	G→A	noncoding (30/59 nt)	sokC →	small regulatory RNA antitoxin SokC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	16,981	0	G	A	93.3%	33.0 / ‑5.8	16	noncoding (30/59 nt)	sokC	small regulatory RNA antitoxin SokC
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (8/6);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCC‑GTCTGGCCAGTTCGGATGTGAACC‑TCACAGAGGT‑CTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGA‑TGGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAAC  >  NC_000913/16736‑17225                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      cccccGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCCCC‑GTCTGGCCAGTTCGGATGTGAACC‑TCACAGAGGT‑CTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑aaa                                                                                                                                                                                                                                                 <  1:180545/250‑1 (MQ=255)            gccctcatacgcgtgcttAAGCCGTGAAAACAGCAGCCCCCTGGC‑GGCCAGTCCGTATGTCATCCGTCA‑AGAGGT‑TTTTCCTCTCTACCAGCTCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGCCAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCCCCTTGCCCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATAGAGGTGGCCTCGGGTTGA‑TTGT‑AAAATTTCACTCgg                                                                                                                                                                                                                                      <  2:353164/234‑1 (MQ=255)                      ggATTCGTAAGCCGTGAAAACATCACCCTCC‑GTCTGGCCAGGTCGGATGTGAACC‑TCACAGAGTT‑CTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTtctc                                                                                                                                                                                                                            <  1:213644/251‑1 (MQ=255)                                   gTGAAAACAGAAACCTCC‑GTCTGGCCAGTTCGGATGTGAACC‑TCACAGAGGT‑CTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAAACCTCGCCTTATGCTGCTTCATTTCTCTCTTCTCCTTGACCTTACCGTCAGTAAGAGGCACTCTAACTGTGTTTAGCATATCGGGGGCCTCGGGTTGA‑‑GATGAAAAAATACCTCGCGGCTTTTCTCTATCTGCCGTTcc                                                                                                                                                                                                               >  2:115968/1‑250 (MQ=255)                                         aCAGCAACCTCC‑GTCGGGCCAGTTCGGATGTGACCC‑TCACAGAGGT‑CTTTTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAAt                                                                                                                                                                                                         <  1:200015/251‑1 (MQ=255)                                                   cccttcgggCCAGTTCGGATGGGAACC‑CCA‑AGGGGTTCTTTTCTCTTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGagac                                                                                                                                                                                                <  2:550978/244‑1 (MQ=255)                                                    ccctCTGTCCAGGTCGGATGTGAACC‑CCACAGAGTT‑CTTTGCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGagaca                                                                                                                                                                                               <  1:352979/248‑1 (MQ=255)                                                                                                     aCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTAAGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGGTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAACACCCGCCGCTATTATATCGCTCTCTTTAAcccct                                                                                                                                                >  1:6746/1‑249 (MQ=255)                                                                                                                   cgcgggggtaaatcggATGACCATCAGGGCTACATTCACCTCCTTCTGTTGCTTCATTCCCCTCTTCCCCTTCATCTCCCGGTCGGAAAGGGCCCCTCTACATGTTTTCAGTAGATGGGGGGCCTCGG‑TTGAATAGA‑AAAATATCCCCCGGGTTTTTTTTCTATCTCCCGTTCAGATAATCCCTAAGACATATAGCCAAAACCTCCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGAt                                                                                                                                  <  1:522452/236‑1 (MQ=255)                                                                                                                      ggcggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCt                                                                                                                               >  1:204489/1‑251 (MQ=255)                                                                                                                                                  aCAATCATCGCCTTATGCTGCTTCATTGATATCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGACAACGCCTCGCAtttttttg                                                                                                   >  2:395651/1‑248 (MQ=255)                                                                                                                                                                         cATTGCTCTCTTCTACGTGAACTTACGGTCAGTAAGAGGCACTCTACCTGTGTTCAGCATATAGAGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTACTCCTGAAGACGCATCGCATATCTGAGGCGTAATTTTTTGATGATTTAat                                                                            >  2:470459/1‑251 (MQ=255)                                                                                                                                                                                ctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTa                                                                      >  1:300693/1‑250 (MQ=255)                                                                                                                                                                                        cTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGctctct                                                                                                                                                         >  1:62723/1‑159 (MQ=255)                                                                                                                                                                                                            gagGCACTCTTCATGTGGTCAGCATATAGGGGGCCTCGGGTTGA‑TAGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCATGAGCCAGATAGCCTCAAGCACCCGCCGCTATTATCTCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGCCGCCGCATCGAATTCTTGTGGCGAATTTTTTTACTGCTTTAATTATTTAACTTTATTTTATATCtttatagtaatttt                                         >  2:150418/1‑239 (MQ=255)                                                                                                                                                                                                                                                    ttGA‑TAGT‑AAAAAATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTGTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCATAGCATTTTTGTGGCGTACTTTTTTCATGGTTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGaaaacttgtcttaac  >  2:249995/1‑238 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      TCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCC‑GTCTGGCCAGTTCGGATGTGAACC‑TCACAGAGGT‑CTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGA‑TGGT‑AAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAAC  >  NC_000913/16736‑17225

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