Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,817,246	G→A	D191N (GAT→AAT)	yicC →	putative RNase adaptor protein YicC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,817,246	0	G	A	84.2%	31.5 / 0.8	19	D191N (GAT→AAT)	yicC	putative RNase adaptor protein YicC
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (9/7);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

GATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGGACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGATCGCCTCGAAGCGCATGTCAAAGAGACCTACAACATTCTGAAGAAAAAAGAAGCGGTTGGTCGTCGTCTGGATTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTG  >  NC_000913/3817020‑3817444                                                                                                                                                                                                                                   |                                                                                                                                                                                                       gATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGGACGGCTTTATTGTCGCGCGGGAAACCGAAGGTCAGGCACTGAAATGATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAAATGCTCAGGTGCAACTGGAAaac                                                                                                                                                                                 <  1:285170/250‑1 (MQ=255)                                      gCTGAAATTCTCGCGGCGCTGGATGGTACGCTGGACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGTTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAAATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTctgc                                                                                                                                                     <  1:156335/241‑1 (MQ=255)                                           aaTTCTCGCGGCGCTGGATGGTACGCNGGACGCNNTTNTTGTCGCGCGCGAAACCGANGGNNNGNNNCNGAAAGNATNNATCGAGCAGNGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAANTCCTGCAATNNCAGCGTGAGCGTCTGGTCGCGNAGCNNGNNANTGCNCan                                                                                                                                                                                                >  2:75772/1‑192 (MQ=255)                                                    ggcgCTGGATGGTACGCTGGACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAAATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCg                                                                                                                             <  1:75772/251‑1 (MQ=255)                                                        cTGGATGGCACGCTGGACGACTTTATTGTCGCGCGCCAAAAAGAAGGTCAGGCACTGAACGCACTGATCGACCAGCGTCTGGACAGCGCCACAGCCGACGTGGTCAAAGTCGGCTCCCATATGCCGGAAGTCCTGACATGGCAGCGTGTGCGTCTGGTCGCGAAGCTGGAAAATGCTCACGTGCAACTGGAAAACAACCGTCTGGAGCTGGAACTGGATCTGCTGGCACAACGAATTGACGTTGCCga                                                                                                                            >  1:322177/1‑248 (MQ=255)                                                            aTGGTACGCTGGACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAAATGCTCAGGTGCAACTGGAAAACATCCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACtn                                                                                                                     <  1:61669/251‑2 (MQ=255)                                                                                   gTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTAGAAAATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACt                                                                                                                                                             >  1:160537/1‑188 (MQ=255)                                                                                                ccGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAAATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGATCGCCTCGAAGCGCATGTCAAAGAGACCTacaac                                                                                 >  1:267121/1‑251 (MQ=255)                                                                                                  naaGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAAATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGATCGCCTCGAAGCGCATGTCAAAGAGACCTACAACAt                                                                               >  1:78741/2‑251 (MQ=255)                                                                                                                          aTCGAGCAGCGTATGGAAGGCGTCAACGCCGAAATGGTCAAAGTCCGCACACATATGCCGGATATCCTGCAAAGGCAGAGTGAGCGTCTGGTCGAGAAGCTGGAAAATGCTCTGGTGCAACTGGACAACAACCGTCTGGAGAAGGAACTGATTCTGCTGGCCCAACGAATTGCCGTTGCCGAAGACCTGGCTCGCCTCGAAGCGCATGTCAAaga                                                                                           >  1:313479/1‑215 (MQ=255)                                                                                                                           tCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAATCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAAATGCTCAGGTGCAACTGGaaa                                                                                                                                                                                   >  1:5210/1‑126 (MQ=255)                                                                                                                                    gTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAAATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTgg                                                                                                                                                  <  1:259288/150‑1 (MQ=255)                                                                                                                                       tGGAAGGCGTCACCGCCGAAGTGGCCAAAGTCCGCTAGAATATGCCGGCACTCCTGATATAGCAGCGCGCGCGTCTGGTCGCGACGCTGCAAAATGCCCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGGTCTGCAGGCAAAACGAAATGACGTTGGCGAAGAACTGGa                                                                                                                   >  1:303587/1‑178 (MQ=255)                                                                                                                                                                           cccATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAAATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGATCGCCTCGAAGCGCATGTCAAAGAGACCTACAACATTCTGAAGAAAAAAGAAGCGGTTGGTCGTCGTCTGGATTTTATGATGCAGGAGTTCAACCGCGAGTCGAACAct      >  1:129109/1‑251 (MQ=255)                                                                                                                                                                              atatGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAAACAACCGtn                                                                                                                                                                          <  1:41152/84‑2 (MQ=255)                                                                                                                                                                                atGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAAATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGATCGCCTCGAAGCGCATGTCAAAGAGACCTACAACATTCTGAAGAAAAAAGAAGCGGTTGGTCGTCGTCTGGATTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCttn  <  1:79190/250‑2 (MQ=255)                                                                                                                                                                                     ggAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAGCGAATTGACGTTGCCGAAGAACTGGATCGCCTCGAAGCGCATGTCAAAGAGACCTACAACATTCTGAAGAAAAAAGAAGCGGTTGGTCGTCGTCTGGATTTTATGATGCAgg                             <  1:15969/218‑1 (MQ=255)                                                                                                                                                                                            cTGCAATGGCAGCGTGCGCGTCTGGTCGCGACGATGGAAGATGCTCAGGTGCAACTGGAAAACAACCGACTGGAGCAGGAACTGGTTct                                                                                                                                                       >  1:265616/1‑89 (MQ=255)                                                                                                                                                                                                                             cTGGAAAATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGATCTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTCGATCGCCTCGAAGCGCATGTCAAAGAGACCTACAACATTCTGAAGAAAAAAGAATCGGTTGGTCGTCGTCTGGATTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCttn  <  1:29299/205‑2 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                       GATGGCAGCCCAGGAGCAGGATCTTGACGCCATTGCCGCTGAAATTCTCGCGGCGCTGGATGGTACGCTGGACGACTTTATTGTCGCGCGCGAAACCGAAGGTCAGGCACTGAAAGCATTGATCGAGCAGCGTCTGGAAGGCGTCACCGCCGAAGTGGTCAAAGTCCGCTCCCATATGCCGGAAATCCTGCAATGGCAGCGTGAGCGTCTGGTCGCGAAGCTGGAAGATGCTCAGGTGCAACTGGAAAACAACCGTCTGGAGCAGGAACTGGTTCTGCTGGCACAACGAATTGACGTTGCCGAAGAACTGGATCGCCTCGAAGCGCATGTCAAAGAGACCTACAACATTCTGAAGAAAAAAGAAGCGGTTGGTCGTCGTCTGGATTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTG  >  NC_000913/3817020‑3817444

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