Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,105,528	G→T	E39* (GAA→TAA)	ymdA →	uncharacterized protein YmdA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,105,528	0	G	T	93.8%	42.8 / NA	16	E39* (GAA→TAA)	ymdA	uncharacterized protein YmdA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (8/7);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGG‑AACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGAAAAAGCTGGCTGTAATGAACATCGAATATAACTAAGTTTCTGAACAACTCACGTTGCTGAGCAGAAAAATGCGATTTA  >  NC_000913/1105289‑1105768                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  cTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGACAGTCTGGATAACACAGCCCTCCAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGAAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGT‑AACCGGCAtgg                                                                                                                                                                                                                                       >  2:153212/1‑250 (MQ=255)                tatCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACATGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGT‑AACCGGCATGTGACGTCAGCAccc                                                                                                                                                                                                                          >  1:297172/1‑249 (MQ=255)                                     ggAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGT‑AACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATg                                                                                                                                                                                                      <  1:224573/248‑1 (MQ=255)                                     ggAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCGCAAGGATGACTAATCATTGTGGAGATAGAATGTATGTTCAGACCTTTTTTAAACTGTCTTAGGCTCGTGAGGTTGTTTTTTCCTTTTATTGTGATTGGTGGAAGCACCGTGCAAGGGGGCGTGATCCAATTTTATGGTCGA‑TTTTGGTAACCGGCTGGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGaa                                                                                                                                                                                                    <  1:302995/250‑1 (MQ=255)                                         agaATGGTCTGCTGCTAATCAGCCAACCTGAAATAATGGATAACTTAACCCTCAAGGATGACTAATCATTGAGTAAATAGAATAAATGTTCAGGCCGTTTTTAAACTCTCTTGGGCTCGTCAGTTTGTTTTTTCTTTTTATTGCCTTTGGTGGATGCTCCGTGCATGGGGGCTTGATCCATTTTTATGGCCAAAGTGTTT‑AACCGGCATTTGGCGTCAACACCCAGTCATCACCCGTAGAAATGAACTGcc                                                                                                                                                                                               <  1:147483/251‑1 (MQ=255)                                                       cGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGT‑AACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCt                                                                                                                                                                                 <  2:297172/251‑1 (MQ=255)                                                                                         cccTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGT‑AACCGGCATGTGACGTCAGCTCCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAGCCTACTCCAGCAAAGCGTTaa                                                                                                                                               <  2:207026/251‑1 (MQ=255)                                                                                                                                         ttttttacattctattctgatCGGCAGTTTGTTTTTTCCATTTATTGTCATTGCTGGAAGCCCTGTGCAAGGGGGCGTGATCCATTTTTATGGCCAACTTGTGT‑AACCGGCATGTGCCGTCCGCAACCAGTACTCCACCGTAGACATGAACTGCCCACAGACTGGTTTTATTCAGGGAAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGACACTAGCCTCAGTAAAAGGGCAGt                                                                                               >  2:315442/21‑251 (MQ=255)                                                                                                                                         ttttttAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGT‑AACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGt                                                                                               >  2:252091/1‑251 (MQ=255)                                                                                                                                               aaaCTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGT‑AACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCtt                                                                                          <  1:148837/250‑1 (MQ=255)                                                                                                                                                                                                     cACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGT‑AACCGGCATGTGACGTCAGCACCCAGTCATCAACCGTAGAAATGCACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGAGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGACAAAGCTGGCTGTAATGAACATCGAATATAACTAAGTTTCTGaac                                   >  2:301550/1‑251 (MQ=255)                                                                                                                                                                                                                      tGATCCATTTTTATGGCCAAATTGTGT‑AACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACTCGCAAATAGCCTCAg                                                                                                           <  2:61367/162‑1 (MQ=255)                                                                                                                                                                                                                      tGATCCATTTTTATGGCCAAATTGTGT‑AACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAg                                                                                                           >  1:61367/1‑162 (MQ=255)                                                                                                                                                                                                                              ttttATGGCCAAATTGTGT‑AACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAACACGCGCAAATAGCCTCCGTaaaa                                                                                                      >  1:308101/1‑159 (MQ=255)                                                                                                                                                                                                                              ttttATGGCCAAAGTGTTT‑AACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACGTCCCACAGAATGGTTCTATTCCTGGCAAAACCTACTCCAGCAAAGCGTTAATGATCGGCAATGTTAGAAACGCGCAAATAGCCTCAGTaaaa                                                                                                      <  2:308101/159‑1 (MQ=255)                                                                                                                                                                                                                                       cAAATTGTGT‑AACCGGCATGTGACGTCAGAACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGAAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGAAAAAGCTGGCTGTAATGAACATCGAATATAACTAAGTTTATGAACAACTAACGTTGCTGAGCAGAAAAATGCGATTTa  >  2:277697/1‑250 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  CTTAATTTGTTGAAATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGG‑AACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGAAAAAGCTGGCTGTAATGAACATCGAATATAACTAAGTTTCTGAACAACTCACGTTGCTGAGCAGAAAAATGCGATTTA  >  NC_000913/1105289‑1105768

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