Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,358,859	T→C	Y110C (TAT→TGT)	puuP ←	putrescine:H(+) symporter PuuP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,358,859	0	T	C	94.1%	48.9 / NA	17	Y110C (TAT→TGT)	puuP	putrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (8/8);  minor base A (1/0);  total (9/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.71e-01

CTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATC  >  NC_000913/1358629‑1359101                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   cTTCTCCTTTATGCAGTCCCTGCACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGGTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGAcc                                                                                                                                                                                                                                  >  1:314259/1‑249 (MQ=255)           aTGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGGATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATa                                                                                                                                                                                                                       <  2:398143/250‑1 (MQ=255)                           aCCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTGGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACTAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGCCGTGCGGGtt                                                                                                                                                                                                      <  2:296487/251‑1 (MQ=255)                            ccAGGAAGATAAACACCTCCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCCCCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTa                                                                                                                                                                                                     <  1:49250/251‑1 (MQ=255)                                gAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACCCCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCg                                                                                                                                                                                                 <  2:107473/251‑1 (MQ=255)                                               caTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGTCCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAgg                                                                                                                                                                                  <  1:383692/251‑1 (MQ=255)                                                                ccAGTACAAACAGGGTATTGAAGTTAGCGCCCAGGCTGCCGCCCTTCAGATTCGCGGCGGTTAAAATGGCGACGCAGGTTACCACCCACACCCAAGGCGGCCCTTCCCGGAAGCGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAAAGACAAACCAGCCGCGAGGACCCGCCGACCATAAATACACCATGCGGGTTAATCGACTTTTGCGCGTCGGCATAGGCCGAACCGGcc                                                                                                                                                                 >  2:50877/1‑251 (MQ=255)                                                                     aCAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCg                                                                                                                                                             <  1:359826/250‑1 (MQ=255)                                                                                  tGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAAcc                                                                                                                                               >  1:261975/1‑251 (MQ=255)                                                                                                                       ggcggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTg                                                                                                                                         <  1:22882/220‑1 (MQ=255)                                                                                                                                         gaAGGTTACCACCCCCACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAgg                                                                                        >  1:107052/1‑251 (MQ=255)                                                                                                                                                                             gAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCCAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCAGACGTGCGGGATAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACa                                                     >  1:201088/1‑250 (MQ=255)                                                                                                                                                                                       gagagaTAGCTTTTCGACAACAAGACGTTAATCATCGGCAAAAAGAGACCATACAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGCGACAATACCAAAAGt                                          >  2:216780/1‑251 (MQ=255)                                                                                                                                                                                                atTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATAc                                 >  1:411326/1‑251 (MQ=255)                                                                                                                                                                                                                     atcaccgaaaaaaaGAGACAACCCAGCAGCGATGACCAGCGGACCATAAATCCGACGTGCGGTTTAACCGCTTTTTGCGCGTAGGTATAGGCCGAACCGGCTCCCGGAAACTGGCGAACCAGTTTGCCTTAGCTGATAGCGGTAAACAGCCCCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGc            <  1:86711/242‑1 (MQ=255)                                                                                                                                                                                                                       catcGGCCAAAAGAGAAAATCCAGCAGCCATCACCAGCCGACCATACATCCGCCGTGCGGGTTAATCGACTTTTGCGCGTAGATATAAGCCGCACCGGCCTCCGGAAACTGGCGCACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCATCGCCATCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTCTCAAATACGGTCATCGGCGTgagagagacc           >  2:42965/1‑244 (MQ=255)                                                                                                                                                                                                                               aaaaGAGACAATCCAGCAGCGATACCCAGCCGGCCCTAAATCCGACGGGCGGGCTAATCCCCTTTTGCGCGCAGGTATAGGCCGCCCAGGCCTCCGGCAACTGGCGAACCAGTTTGACGTAGCTGATAGAGGCCAACAGCACACCCGCCCGCGCCAGCACATAAGCCGCCGGAACGTGACCGTCGCTAGTGCCGGCGACAATACCAAAAGTATCAAATACGGTCATCGGCGGGTGATAGGCCCGACCcatc  >  2:148577/1‑251 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   CTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATC  >  NC_000913/1358629‑1359101

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