Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007790	1,967	C→G	intergenic (+200/+126)	SAUSA300_RS14710 → / ← SAUSA300_RS14715	membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007790	1,967	0	C	G	85.0%	40.3 / 2.7	21	intergenic (+200/+126)	SAUSA300_RS14710/SAUSA300_RS14715	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base G (13/4);  total (17/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

TTTGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTG‑ATGTTGAGAAGAACCCTTAACTAAACTTGCAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACAAGTTTTGGGATTGTTAAGGGTTC  >  NC_007790/1873‑2052                                                                                                |                                                                                      tttGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCtt                                                                                   >  1:4612656/1‑100 (MQ=255) tttGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCtt                                                                                   >  1:3154310/1‑100 (MQ=255) cttGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCtt                                                                                   <  1:3636334/99‑1 (MQ=255)  ttGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTa                                                                                  >  1:265200/1‑100 (MQ=255)   tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTaa                                                                                 >  1:1658630/1‑100 (MQ=255)   tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTaa                                                                                 >  1:730542/1‑100 (MQ=255)   tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTaa                                                                                 >  1:2248526/1‑100 (MQ=255)   tGACGAGATGGAAGAAAAATATGATTATCATGAAGATAAT‑AAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAAc                                                                                >  1:734945/1‑100 (MQ=255)    gACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAAc                                                                                >  1:2660365/1‑100 (MQ=255)     ccGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATTCCGAGAAAATTTATTGTGCGTTGAGAAGAAACCTTAACt                                                                               >  1:5054260/2‑100 (MQ=255)      cGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTa                                                                              >  1:2168409/1‑100 (MQ=255)      cGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTa                                                                              >  1:5200938/1‑100 (MQ=255)            ggAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTAATCTTg                                                                        >  1:5982382/1‑100 (MQ=255)             gaagaaAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTAAACTTgg                                                                       <  1:2294083/100‑2 (MQ=255)             gaagaaAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTAAACTTgg                                                                       <  1:3322774/100‑2 (MQ=255)             gaagaaAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTAAACTTgg                                                                       <  1:2090425/100‑2 (MQ=255)                                                    ccgacgCATTTATGCCGAGAAAATTTATTG‑ATGTTGAGAAGAACCCTTAACTAAACTTGAAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCAt                               >  1:390854/2‑100 (MQ=255)                                                          gccccgtGCCGAGAAAATTTATTG‑ATGTTGAGAAGAACCCTTAACTAAACTTGAAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACa                         >  1:1045771/7‑100 (MQ=255)                                                               aTGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTAAACTTGGAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACAAGtt                     >  1:2800528/1‑100 (MQ=255)                                                               aTGCCGAGAAAATTTATTGTGCGTTGAGAAGAAGCCTTAACTAAACTTGGAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACAAGtt                     >  1:301323/1‑100 (MQ=255)                                                                                 g‑ATGTTGAGAAGAACCCTTAACTAAACTTGAAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACAAGTTTTGGGATTGTTAAGGGTTc  >  1:2049674/1‑100 (MQ=255)                                                                                                |                                                                                      TTTGACGAGATGGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTG‑ATGTTGAGAAGAACCCTTAACTAAACTTGCAGACGAATGTCGGCATAGCGTGAGCTATTAAGCCGACCATTCGACAAGTTTTGGGATTGTTAAGGGTTC  >  NC_007790/1873‑2052

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