Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,000,230	A→G	noncoding (468/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,230	0	A	G	95.2%	95.6 / ‑0.0	83	noncoding (468/2926 nt)	SAUSA300_RS10040	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (75/4);  total (76/7)
Fisher's exact test for biased strand distribution p-value = 1.47e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

TTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCA‑CTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTA  >  NC_007793/2000171‑2000323                                                             |                                                                                              ttccgatcTTTCCCTCACGGTACTGGTTCACTATCGGTCA‑CTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACgtgtg                                                       <  1:3823929/93‑2 (MQ=2)  tttCGCCTTTCCCTCACGGTACTGGCTTACTATCGGTCA‑GTCAGGAGTATTTAGCCTCGGAGGATGGTCCCCCCATATTCAGACAGGATACCACgtgtcc                                                      >  1:590397/1‑98 (MQ=2)  cctcgccTTTCCCTCACGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACgtgtcc                                                      >  1:3021361/6‑98 (MQ=2)    tcgccTTTCCCTCACGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCcg                                                    >  1:1160797/4‑100 (MQ=2)    gccccTTTCCCTCACGGTACTGGATTACTATCGGTCA‑GTCTGGAGTAATTAGCCTTGGAGGAAGGTCCCCCCATATTCAGACAGGATACCACGCGTCCcg                                                    >  1:2090595/4‑100 (MQ=0)     cgccTTTCCCTCACGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCcgt                                                   >  1:2651954/3‑100 (MQ=2)     cACCTTTCCCTCACAGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCcgc                                                   >  1:2441902/1‑99 (MQ=11)      aCCTTTCCCTCACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCTGACGGAATTTCACGTGTTCcgcc                                                  >  1:2891429/1‑98 (MQ=11)         tgtCCCTCACGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCcgccctc                                               <  1:114578/98‑6 (MQ=2)          ttCCCTCACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTAc                                              >  1:4917542/1‑100 (MQ=11)           tCCCTCACGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTAcc                                             >  1:468391/1‑99 (MQ=2)           ccccTCACGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCcgccccact                                             >  1:4114781/2‑93 (MQ=2)            gatctgtttcaCTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCCGCCCTACTc                                            <  1:988709/90‑1 (MQ=2)            cccTCACGGTACTGGTTCACTATCGGTCA‑GTCAGGAGTATTTAGCCTTGGAGGATGGTCCCCCCATATTCAGACAGGATACCACGTGTCCcgcccaacta                                            >  1:5957671/1‑92 (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:6019464/1‑96 (MQ=11)                                   cGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagatcgg                     >  1:1203417/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:3315205/1‑92 (MQ=11)                                       cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtt                 >  1:1573898/1‑92 (MQ=11)                                       cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtt                 >  1:6352830/1‑92 (MQ=11)                                       cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtt                 >  1:558801/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: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:3726397/1‑96 (MQ=11)                                                       aGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAgggctg  >  1:2429999/1‑96 (MQ=11)                                                       aGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGCGTTCCGCCGTACTCAGGATCCACTCAGGAGAGAACGAAGTTTTGACTACAggggca  >  1:2423555/1‑96 (MQ=11)                                                             |                                                                                              TTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCA‑CTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTA  >  NC_007793/2000171‑2000323

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