Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,044,642	T→C	intergenic (+18/‑296)	asnT → / → yeeJ	tRNA‑Asn/inverse autotransporter adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,044,642	0	T	C	93.8%	48.5 / ‑5.6	16	intergenic (+18/‑296)	asnT/yeeJ	tRNA‑Asn/inverse autotransporter adhesin
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (5/10);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTTTCGGCGACGAATGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCAGTATT  >  NC_000913/2044423‑2044876                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           cGTTTTCGGCGACGAATGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGa                                                                                                                                                                                                              <  2:191180/250‑1 (MQ=255)               aaTGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAAc                                                                                                                                                                                               <  1:68595/251‑1 (MQ=255)                             aacaccTTTGCAGCTTAATTAACCAATTGAAATGACTTATGAAATTTCGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCAATGGTTCGAGTCCAGTCAGCGGAGCCAAATTCTAAAAATTCGCTTTTTTTGCGCAATGTCCATGAACTTAGTTGAACATTGTTTTTTACAGGATAGCGGGattt                                                                                                                                                                                 >  2:180536/1‑251 (MQ=25)                                         gATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTa                                                                                                                                                                       >  2:44164/1‑249 (MQ=255)                                              attaaCCAATCGAAATGACTTATGAAATTTAGCGTTGACAGACAAGGTACCGATAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGACCGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATATCGGGTTTTTAACATCTTAAGCGccc                                                                                                                                                                >  2:237256/1‑251 (MQ=255)                                                taaCCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGAGTTTTTAACATCTTAAGCGCCCt                                                                                                                                                               <  2:39136/250‑1 (MQ=255)                                                               aCTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATCTTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCTGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTg                                                                                                                                                >  1:121725/1‑250 (MQ=255)                                                                                                    aaGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTAGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTcccc                                                                                                           <  1:167927/250‑1 (MQ=255)                                                                                                                           cGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAGAATTCGCTCTTTTAGCGCAATGTCACTGACCTTATTTGAACATTGTTGTTTAACGGATAGCGTTTTTTTAACATCTTAAGCGGCCTTGTCCGCTAGGGGTGAGGGCGTTTTTCTTTGAGCTCTGTCTGCATTTTCCCCTCGATTATAAAACTTGAGTTAtcc                                                                                   <  1:100576/251‑3 (MQ=21)                                                                                                                                             gTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGc                                                                 >  1:131145/1‑251 (MQ=255)                                                                                                                                                                   aaTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTCCAACTCACACCCAAAACTGCCTa                                           <  2:11935/251‑1 (MQ=255)                                                                                                                                                                                cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGttatt                              >  1:186836/1‑251 (MQ=255)                                                                                                                                                                                  ggTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCGCCATTTGCCCCTCGATTATAAAGCTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGttatta                             <  1:237256/250‑1 (MQ=255)                                                                                                                                                                                                cAGAGGAGCCAAATTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTGATTAATTTTCAATCatat               <  2:111464/250‑1 (MQ=255)                                                                                                                                                                                                           aaTTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCAGTAt   <  2:126406/251‑1 (MQ=255)                                                                                                                                                                                                            aTTCTAAAAATTCGCTCTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCAGTAtt  <  2:202142/251‑1 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           CGTTTTCGGCGACGAATGTTCATTAATTAACAACTTTGCAGATTAATTAACCAATTGAAATGACTTATGAAATTTAGTGTTGACAGACAAGGTACCGCTAAGTAATATTCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCTAAAAATTCGCTTTTTTAGCGCAATGTCACTGACCTTAGTTGAACATTGTTTTTTAACGGATAGCGGGTTTTTAACATCTTAAGCGCCCTCGACCTTTATGGTTGAGGGCGTTTTGCTATGAACGCCATCACCATTTTCCCCTCGATTATAAAACTTGAGTTATTCAGTAGTCTCCCCTCTTGCAACTCACACCCAAAACTGCCTAACGAAAAGTTATTAATTTTCAATCATATTGCTATCAGTATT  >  NC_000913/2044423‑2044876

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