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	95.8%	41.4 / ‑2.7	49	noncoding (424/2926 nt)	SAUSA300_RS10040	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (45/1);  total (47/2)
Fisher's exact test for biased strand distribution p-value = 8.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

TCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATGTCG  >  NC_007793/2000181‑2000367                                                                                                |                                                                                              ccccTCACGGTACTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACt                                                                                         >  1:435423/2‑100 (MQ=2)    tccACGGTACTGGTTCACTATCGGTCACT‑AGGGAGTATTTAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAg                                                                                      >  1:2421894/3‑100 (MQ=11)    cTCACGGTACTGGTTCACTATCAGTC‑‑‑‑AG‑GAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTAcagatcgga                                                                                  >  1:3323828/1‑92 (MQ=2)     tCACGGTACTGGTTCACTATCGGTCACT‑AGGGAGTATTTAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAgg                                                                                     >  1:2821083/1‑100 (MQ=11)      taCGGTACTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTTGGAT‑GATGGTCCCCCCATATTGAGACAGGATAGCACGTGTCCCGCCCTACTCatcg                                                                                    >  1:2292046/2‑97 (MQ=2)      cACGGTACTGGTTCACTATCGGTCAGTCAG‑GAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTAcaagatc                                                                                    >  1:305452/1‑94 (MQ=2)      cACGGTACTGGTTCACTATCGGTCAGTCAG‑GAGTATTTAGCCTT‑GGAGGATGGTCCCCTCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcg                                                                                    >  1:4201089/1‑97 (MQ=2)      cACGGTACTGGTTCACTATCGGTCAGTCAG‑GAGTATTTAGCCTT‑GGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGTCCTAcagatcg                                                                                    >  1:2734698/1‑94 (MQ=2)      cACGGTACTGGTTCACTATCGGTCAGTCAG‑GAGTATTTAGCCTT‑GGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTAcaagatc                                                                                    >  1:5038482/1‑94 (MQ=2)      cACGGTACTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcg                                                                                    >  1:3619821/1‑97 (MQ=2)      cACGGTACTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcg                                                                                    >  1:5395560/1‑97 (MQ=2)      cACGGTACTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCACCCTACTCatcg                                                                                    >  1:755000/1‑97 (MQ=2)      cACGGTACTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTT‑GGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTAcacaacg                                                                                    >  1:5427497/1‑97 (MQ=2)      cACGGTACTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTT‑GGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACtagatc                                                                                    >  1:2102161/1‑95 (MQ=2)      cACGGTACTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTT‑GGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACtagatc                                                                                    >  1:2921799/1‑95 (MQ=2)      cACGGTACTGGTTCACTATCGGTCAGT‑CAGGAGAATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcg                                                                                    >  1:503572/1‑97 (MQ=2)      cACGGTACTGGTTCACTATCGGTCACT‑AGGGAGTATTTAGCCTTGGGA‑GATGGTCCTCCCTGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAggg                                                                                    >  1:492853/1‑99 (MQ=11)      cACGGTACTGGTTCACTATCGGTCACT‑AGGGAGTATTTAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGa                                                                                    >  1:4165811/1‑100 (MQ=11)      cACGGTACTGGTTCACTATCGGTCACT‑AGGGAGTATTTAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGa                                                                                    >  1:1209161/1‑100 (MQ=11)      cACGGTACTGGTTCACTATCGGTCACT‑AGGGAGTATTTAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGa                                                                                    >  1:3508281/1‑100 (MQ=11)      cACGGTACTGGTTCACTATCGGTCACT‑AGGGAGTATTTAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGa                                                                                    >  1:4237151/1‑100 (MQ=11)      cACGGTACTGGTTCACTATCGGTCACT‑AGGGAGTATTTAGCCTTGGGA‑GATGGTCCTCCCGGATACCGACGGAATTTCACGTGTTCCGGCGTACTCAGGa                                                                                    >  1:2365851/1‑100 (MQ=11)       aCGGTACTGGTTCACTATCGGTCACT‑AGGGAGTATTTAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGAt                                                                                   >  1:2161045/1‑100 (MQ=11)        aggTACTGGTTCACTATCGGTCAGTCAG‑GAGTATTTAGCCTT‑GGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgag                                                                                  >  1:4398470/2‑95 (MQ=2)         cgatcTGGTTTACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTAATCAGGATcc                                                                                 <  1:1913881/96‑1 (MQ=255)         agTACTGGTTCACTATCGGTCAGTCAG‑GAGTATTTAGCCTT‑GGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagc                                                                                 >  1:3481676/2‑94 (MQ=2)          gTACTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagca                                                                                >  1:4082565/1‑93 (MQ=2)          ctACTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagct                                                                                >  1:3120111/2‑93 (MQ=2)          ccaCTGGTTCACTATCGGTCAGTCAG‑GAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagct                                                                                >  1:3353365/3‑93 (MQ=2)           tACTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagctc                                                                               >  1:730249/1‑92 (MQ=2)           tACTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagctc                                                                               >  1:2384307/1‑92 (MQ=2)            gcTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTTGGAG‑GATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagctca                                                                              >  1:4640889/2‑91 (MQ=2)            aCTGGTTCACTATCGGTCAGT‑CAGGAGTATTTAGCCTT‑GGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagctca                                                                              >  1:3099935/1‑91 (MQ=2)                  tCACTATCGGTCACT‑AGGGAGTATTTAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTTAggag                                                                        <  1:1255614/100‑1 (MQ=11)                         cGGTCACT‑AGGGAGTATTTAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCGGgagagatcgg                                                                 >  1:559528/1‑96 (MQ=11)                            ccACT‑AGGGAGTATTTAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagt                                                              >  1:1302916/2‑93 (MQ=11)                                         aTTTAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGACCCACTCAGGAGAGAACGAAGTTTTGACTACAgg                                                  >  1:3966395/1‑100 (MQ=11)                                           ttAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggc                                                >  1:2052484/1‑98 (MQ=11)                                           ttAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggc                                                >  1:5207377/1‑98 (MQ=11)                                           ttAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggc                                                >  1:1188375/1‑98 (MQ=11)                                            tAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggct                                               >  1:1987959/1‑97 (MQ=11)                                            tAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggct                                               >  1:1831766/1‑97 (MQ=11)                                            tAGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggct                                               >  1:587029/1‑97 (MQ=11)                                             cgCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggctg                                              >  1:4134688/2‑96 (MQ=11)                                             aGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggctg                                              >  1:2431906/1‑96 (MQ=11)                                             aGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggctg                                              >  1:4082830/1‑96 (MQ=11)                                             aGCCTTGGGA‑GATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggcag                                              >  1:859698/1‑96 (MQ=11)                                             aGCCTTGGGA‑GATGGGCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggctg                                              >  1:1710462/1‑96 (MQ=11)                                                                                          gtTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAGGGCTGTTACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATgtcg  >  1:4939747/1‑100 (MQ=2)                                                                                                |                                                                                              TCCCTCACGGTACTGGTTCACTATCGGTCACT‑AGAGAGTATTTAGCCTTAGGA‑GATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATGTCG  >  NC_007793/2000181‑2000367

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑