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	96.0%	90.0 / ‑2.6	76	noncoding (424/2926 nt)	SAUSA300_RS10040	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (69/4);  total (69/7)
Fisher's exact test for biased strand distribution p-value = 4.98e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

CACCTTTCCCTCACGGTACTGGTTCACTATCGGTCA‑CTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATGTC  >  NC_007793/2000175‑2000366                                                                                                     |                                                                                             cgccTTTCCCTCACGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCcgt                                                                                              >  1:2651954/3‑100 (MQ=2)      ttCCCTCACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTAc                                                                                         >  1:4917542/1‑100 (MQ=11)       tCCCTCACGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTAcc                                                                                        >  1:468391/1‑99 (MQ=2)        gatctgtttcaCTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTc                                                                                       <  1:988709/90‑1 (MQ=2)        cccTCACGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTc                                                                                       >  1:4800911/1‑100 (MQ=2)        cccTCACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTc                                                                                       >  1:853752/1‑100 (MQ=11)         ccTCACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCa                                                                                      >  1:5197207/1‑100 (MQ=11)          cTCACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTATTCCGCCGTACTCAg                                                                                     >  1:1926016/1‑100 (MQ=11)           gcACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACCTGTTCCGCCGTACTCAgg                                                                                    >  1:1346521/2‑100 (MQ=11)           ccACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAgg                                                                                    >  1:372028/2‑100 (MQ=11)            cACGGTACTGGTTCACTATCGGTCAGCCAG‑GAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcg                                                                                   >  1:2795481/1‑97 (MQ=2)            cACGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCccgccgctagatcgg                                                                                   >  1:216811/1‑91 (MQ=2)            cACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCTTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGa                                                                                   >  1:3391326/1‑100 (MQ=11)            cACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGa                                                                                   >  1:6062109/1‑100 (MQ=11)            cACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGa                                                                                   >  1:2953805/1‑100 (MQ=11)            cACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGAATTTAGCCTTGGGAGATGGTTCTCCCGGAGTCCGACGGAATGTCACGTGTTCCGCCGTACTCAGGa                                                                                   >  1:4794273/1‑100 (MQ=11)            cACGGTACTGGTTCACTATCGGGCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTAAgat                                                                                   >  1:4665685/1‑98 (MQ=2)              cGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTAcacatcgag                                                                                 >  1:1335003/1‑95 (MQ=2)              cGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgaa                                                                                 <  1:2705486/100‑6 (MQ=2)                gTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagct                                                                               >  1:6420362/1‑93 (MQ=2)                 gaCTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCAc                                                                              >  1:358303/2‑100 (MQ=11)                  aCTGGTTCACTATCGGTCA‑GTCCGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagctca                                                                             >  1:2630865/1‑91 (MQ=2)                  aCTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagctca                                                                             >  1:1304660/1‑91 (MQ=2)                  aCTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagccag                                                                             >  1:162266/1‑91 (MQ=2)                  aCTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGCGTCCCGCCCTACTCatcgagccag                                                                             >  1:2330284/1‑91 (MQ=2)                  aCTGGTTCACTATCGGTCA‑CTAGGGAGTCTTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCaga                                                                             >  1:2690411/1‑98 (MQ=11)                   cTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTCatcgagctcgc                                                                            <  1:5634390/100‑11 (MQ=2)                     ggTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCTACTCag                                                                          >  1:4911710/1‑99 (MQ=11)                       tttgcTATCGGTCA‑CTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAagg                                                                        <  1:3194600/96‑2 (MQ=17)                        tCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTgaggag                                                                       >  1:1178917/1‑100 (MQ=11)                        tCACTATCGGTCA‑CTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAagag                                                                       <  1:2582828/100‑1 (MQ=32)                         cACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTGCTCAGGATCCTCTCagatcg                                                                      >  1:3990049/1‑95 (MQ=11)                         cACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTcggagag                                                                      >  1:3235794/1‑94 (MQ=11)                         cACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAggaga                                                                      >  1:2892153/1‑100 (MQ=11)                             aggggTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaac                                                                  >  1:5349494/4‑98 (MQ=11)                               cGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAggagatcggaa                                                                >  1:3891825/1‑94 (MQ=11)                               cGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagatcgg                                                                >  1:1203417/1‑96 (MQ=11)                               cGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagatcgg                                                                >  1:6019464/1‑96 (MQ=11)                               cGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgacagatcgg                                                                >  1:3708010/1‑96 (MQ=11)                               cGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGACCCACTCGGgagagatcgg                                                                >  1:3717218/1‑96 (MQ=11)                               cGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGCCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTGAGgagagatcgg                                                                >  1:1275240/1‑96 (MQ=11)                               cGGTCA‑CTAGGGAGTATTTAGCCTTGCGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCTGgagagatcgg                                                                >  1:190647/1‑96 (MQ=11)                               cGGTCA‑CTAGAGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCCGgagagatcgg                                                                >  1:3092343/1‑96 (MQ=11)                               atctCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacga                                                                >  1:373623/4‑96 (MQ=11)                                tgTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagatcgga                                                               >  1:3728990/2‑95 (MQ=11)                                 gccA‑CTAGGGAATATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaag                                                              >  1:4182283/3‑94 (MQ=11)                                 gTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCTGgagagatcggaa                                                              >  1:3778917/1‑94 (MQ=11)                                 gTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagatcggaa                                                              >  1:1241669/1‑94 (MQ=11)                                 gTCA‑CTAGGGAGTATTTAGCCTTAGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCTGgagagatcggac                                                              >  1:4303112/1‑94 (MQ=11)                                  ccA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagt                                                             >  1:4460257/2‑93 (MQ=11)                                   cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgcgagatcggaaga                                                            >  1:6339887/1‑92 (MQ=11)                                   cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtt                                                            >  1:6352830/1‑92 (MQ=11)                                   cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtt                                                            >  1:1573898/1‑92 (MQ=11)                                   cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtt                                                            >  1:558801/1‑92 (MQ=11)                                   cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtt                                                            >  1:3315205/1‑92 (MQ=11)                                   cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCT‑‑CGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtttt                                                          >  1:1906961/1‑90 (MQ=11)                                      cTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtttt                                                          >  1:2151714/1‑90 (MQ=11)                                               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:2359160/1‑97 (MQ=11)                                                  tAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggct                                              >  1:1070083/1‑97 (MQ=11)                                                  tAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggct                                              >  1:2566322/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:3726397/1‑96 (MQ=11)                                                   aGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggctg                                             >  1:2429999/1‑96 (MQ=11)                                                   aGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGCGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAggggca                                             >  1:2423555/1‑96 (MQ=11)                                                              aTGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAGGGCTGTTAcctcctat                                  >  1:5285011/1‑98 (MQ=2)                                                                   cctccCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAGGGCTGTTACCTCCTATggcgg                             >  1:1439884/1‑96 (MQ=2)                                                                   cctccCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAGGGCTGTTACCTCCTATggcgg                             >  1:2924857/1‑96 (MQ=2)                                                                      cccGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAGGGCTGTTACCTCCTATggcgggcc                          <  1:1450552/100‑8 (MQ=2)                                                                                              tgtTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATgtc  <  1:305466/100‑1 (MQ=32)                                                                                                     |                                                                                             CACCTTTCCCTCACGGTACTGGTTCACTATCGGTCA‑CTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATTACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATGTC  >  NC_007793/2000175‑2000366

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