Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,897,441	A→G	V277V (GTT→GTC)	SAUSA300_RS09370 ←	diaminohydroxyphosphoribosylaminopyrimidine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,897,441	0	A	G	95.0%	42.4 / ‑4.4	20	V277V (GTT→GTC)	SAUSA300_RS09370	diaminohydroxyphosphoribosylaminopyrimidine deaminase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (2/17);  total (2/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

AGATCCGCCAATTAATTTCGGGGCATAATAGAGAATAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCAACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTGAAGAATTGTTGTTAAATCGCAAGACTTCAAGTAAATAATTTCAATATGTGTTTGAT  >  NC_007793/1897348‑1897535                                                                                              |                                                                                               aGATCCGCCAATTAATTTCGGGGCATAATAGAGAATAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCGACTAGc                                                                                          <  1:2560568/100‑1 (MQ=255)  atctcGCCAATTAATTTCGGGGCATAATAGAGAATAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCGACTAGCa                                                                                         <  1:2109940/96‑1 (MQ=255)   tctcGCCATTTAATTTCGGGGCATAATAGAGAATAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCGACTAGCaa                                                                                        <  1:1399194/97‑1 (MQ=255)   aTCTGCCTATTAATTTCGGGGCATAATAGAGAATAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCGACTAGCaa                                                                                        <  1:750616/100‑1 (MQ=255)   aTCTGCCAGTTAATTTCGGGGCATAATAGAGAATAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCGACCAGCaa                                                                                        <  1:3607302/100‑1 (MQ=255)   aTCCGCCAATTAATTTCGGGGCATAATAGAGAATAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCGACTAGCaa                                                                                        <  1:2938581/100‑1 (MQ=255)   aTCCGCCAATTAATTTCGGGGCATAATAGAGAATAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCAACTAGCaa                                                                                        <  1:4837597/100‑1 (MQ=255)    tctgCCGATTAATTTCGGGGCATAATAGGGAATAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTGGTGGTTGGACCTGCCTCGACTAGCaaa                                                                                       <  1:1651174/97‑1 (MQ=255)                                  aaTAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTATTGGTTGGACCTGCCTCGACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTga                                                         <  1:1852702/100‑1 (MQ=255)                                  aaTAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCGACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTga                                                         <  1:759310/100‑1 (MQ=255)                                  aaTAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCGACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTga                                                         <  1:601536/100‑1 (MQ=255)                                   aTAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGACTCGACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTgaa                                                        <  1:860781/100‑1 (MQ=255)                                    tAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCGACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTgaag                                                       <  1:52867/100‑1 (MQ=255)                                              tatataATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCGACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTGAAGAATTgttgtt                                             <  1:1311139/100‑1 (MQ=255)                                              tatataATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCGACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTGAAGAATTgttgtt                                             <  1:1546376/100‑1 (MQ=255)                                                                                     aCCTGCCTCGACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTGAAGAATTGTTGTTAAATCGCAAGACTTCAAGTAAATAATTTCAATATGTGtt      <  1:588032/100‑1 (MQ=255)                                                                                       cTGCCTCGACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTGAAGAATTGTTGTTAAATCGCAAGACTTCAAGTAAATAATTTCAATATGTGTTTg    <  1:5159369/100‑1 (MQ=255)                                                                                       cTGCCTCGACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTGAAGAATTGTTGTTAAATCGCAAGACTTCAAGTAAATAATTTCAATATGTGTTTg    <  1:1031447/100‑1 (MQ=255)                                                                                        tGCCTCGACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTGAAGAATTGTTGTTAAATCGCAAGACTTCAAGTAAATAATTTCAATATGTGTTTGa   >  1:3794031/1‑100 (MQ=255)                                                                                         gCCTCGACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTGAAGAATTGTTGTTAAATCGCAAGACTTCAAGTAAATAATTTCAATATGCGTTAGAt  >  1:3903659/1‑100 (MQ=255)                                                                                              |                                                                                               AGATCCGCCAATTAATTTCGGGGCATAATAGAGAATAAATTCATCTATATAATTAGATTGGAGAAATTCTGAAGTAGTGGTTGGACCTGCCTCAACTAGCAAAGTTCCAACTCCTCTTTTATATAAATTGTGAAGAATTGTTGTTAAATCGCAAGACTTCAAGTAAATAATTTCAATATGTGTTTGAT  >  NC_007793/1897348‑1897535

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