Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,000,306	T→G	noncoding (392/2926 nt)	SAUSA300_RS10040 ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,000,306	0	T	G	95.2%	19.9 / ‑3.8	21	noncoding (392/2926 nt)	SAUSA300_RS10040	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (19/1);  total (19/2)
Fisher's exact test for biased strand distribution p-value = 9.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

ATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATGTC  >  NC_007793/2000220‑2000366                                                                                       |                                                             aTTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgg                                                 >  1:2042896/1‑100 (MQ=11) aTTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAga                                                 >  1:1718859/1‑99 (MQ=11) aTTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACATGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgg                                                 >  1:2371059/1‑100 (MQ=11)  ttcaGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAggg                                                >  1:2840595/4‑99 (MQ=11)   ttAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCATGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggc                                               >  1:2641407/1‑98 (MQ=11)   ttAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTTAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggc                                               >  1:6303581/1‑98 (MQ=11)   ttAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACCCAGGAGAGAACGAAGTTTTGACTACAgggc                                               >  1:531210/1‑98 (MQ=11)    tAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggct                                              >  1:2566322/1‑97 (MQ=11)    tAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggct                                              >  1:1070083/1‑97 (MQ=11)    tAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggct                                              >  1:2359160/1‑97 (MQ=11)    tAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggca                                              >  1:4383481/1‑97 (MQ=11)    tAGCCTTGGAAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggct                                              >  1:33654/1‑97 (MQ=11)     aGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGGGAACGAAGTTTTGACTACAgggctg                                             >  1:2268379/1‑96 (MQ=11)     aGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggctg                                             >  1:2429999/1‑96 (MQ=11)     aGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggctg                                             >  1:3726397/1‑96 (MQ=11)     aGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGCGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAggggca                                             >  1:2423555/1‑96 (MQ=11)                aTGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAGGGCTGTTAcctcctat                                  >  1:5285011/1‑98 (MQ=2)                     cctccCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAGGGCTGTTACCTCCTATggcgg                             >  1:2924857/1‑96 (MQ=2)                     cctccCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAGGGCTGTTACCTCCTATggcgg                             >  1:1439884/1‑96 (MQ=2)                        cccGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAGGGCTGTTACCTCCTATggcgggcc                          <  1:1450552/100‑8 (MQ=2)                                                tgtTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATgtc  <  1:305466/100‑1 (MQ=32)                                                                                       |                                                             ATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATGTC  >  NC_007793/2000220‑2000366

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