Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,000,215	A→G	noncoding (483/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,215	0	A	G	83.3%	15.7 / 2.1	24	noncoding (483/2926 nt)	SAUSA300_RS10040	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base G (20/0);  total (20/4)
Fisher's exact test for biased strand distribution p-value = 9.41e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACAT  >  NC_007793/2000173‑2000306                                           |                                                                                            ctCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTcc                                    >  1:3153362/2‑100 (MQ=255)   cACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCcgc                                  >  1:2326834/1‑99 (MQ=255)   cACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCcgc                                  >  1:4677456/1‑99 (MQ=255)           ccTCACGGTACTGGTTCACTATCGGTCAGTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCATACCCa                          >  1:1104828/1‑100 (MQ=2)            cTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAg                         >  1:756687/1‑100 (MQ=11)            cTCACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACCCAg                         >  1:3036603/1‑100 (MQ=11)             tCACGGTACTGGTTCACTATCGGTCAGTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCaag                        >  1:5386454/1‑98 (MQ=2)             gtaCGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAgg                        >  1:219481/3‑100 (MQ=11)              cttgcatcTGGTTCACTATCGGTCACTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAggt                       <  1:4443348/93‑2 (MQ=255)              cACGGTACTGGTTCACTATCGGTCAGTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcg                       >  1:1481229/1‑97 (MQ=2)              cACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGa                       >  1:1335329/1‑100 (MQ=11)              cACGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGa                       >  1:568025/1‑100 (MQ=11)               aCGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt                      >  1:3634507/1‑100 (MQ=11)               aCGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt                      >  1:4028551/1‑100 (MQ=11)                 cttACTGGTTCACTATCGGTCAGTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagc                    >  1:3790264/3‑94 (MQ=2)                   tACTGGTTCACTATCGGTCAGTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagctg                  >  1:1493714/1‑92 (MQ=2)                     cTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTc                >  1:3626209/1‑100 (MQ=11)                       ggTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCag              >  1:326011/1‑99 (MQ=11)                          tccgatctggTCACTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAagag           <  1:5073628/92‑1 (MQ=255)                          tCACTATCGGTCACTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAagag           <  1:3316234/100‑1 (MQ=32)                           cACTATCGGTCACTAGGGAGTATTTAGGCTTAGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAggaga          >  1:3511844/1‑100 (MQ=11)                           cACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAggaga          >  1:4653379/1‑100 (MQ=11)                           cACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGCCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCGggaga          >  1:97870/1‑100 (MQ=11)                                   gTCACTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAgagagagctttg  <  1:3806317/100‑7 (MQ=12)                                           |                                                                                            TTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCACTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACAT  >  NC_007793/2000173‑2000306

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