Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,000,274	T→C	noncoding (424/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,274	0	T	C	85.2%	16.3 / 2.6	27	noncoding (424/2926 nt)	SAUSA300_RS10040	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base C (23/0);  total (23/4)
Fisher's exact test for biased strand distribution p-value = 5.70e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

CCTCACGGTACTGGTTCACTATCGGTCACTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTA  >  NC_007793/2000183‑2000323                                                                                            |                                                  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:4028551/1‑100 (MQ=11)     aCGGTACTGGTTCACTATCGGTCACTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt                                       >  1:3634507/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)                                         tAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggct   >  1:1060006/1‑97 (MQ=11)                                         tAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggct   >  1:461287/1‑97 (MQ=11)                                         tAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAAGatc   >  1:91161/1‑99 (MQ=11)                                          aGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggctg  >  1:2904590/1‑96 (MQ=11)                                          aGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggctg  >  1:2636341/1‑96 (MQ=11)                                          aGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggcag  >  1:2359631/1‑96 (MQ=11)                                                                                            |                                                  CCTCACGGTACTGGTTCACTATCGGTCACTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTA  >  NC_007793/2000183‑2000323

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