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	93.9%	68.4 / 0.3	67	noncoding (483/2926 nt)	SAUSA300_RS10040	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (58/4);  total (60/7)
Fisher's exact test for biased strand distribution p-value = 2.91e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

TTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCA‑CTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTC  >  NC_007793/2000171‑2000310                                              |                                                                                                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:2953805/1‑100 (MQ=11)                cACGGTACTGGTTCACTATCGGTCA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGa                           >  1:6062109/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:558801/1‑92 (MQ=11)                                       cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtt    >  1:3315205/1‑92 (MQ=11)                                       cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtt    >  1:6352830/1‑92 (MQ=11)                                       cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtt    >  1:1573898/1‑92 (MQ=11)                                       cA‑CTAGGGAGTATTTAGCCTTGGGAGATGGTCCT‑‑CGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtttt  >  1:1906961/1‑90 (MQ=11)                                          cTAGGGAGTATTTAGCCTTGGGAGATGGTCCTCCCGGATTCCGACGGAATTTCACGTGTTCCGCCGTACTCAGGATCCACTCAGgagagaacgaagtttt  >  1:2151714/1‑90 (MQ=11)                                              |                                                                                                TTTTCACCTTTCCCTCACGGTACTGGTTCACTATCGGTCA‑CTAGAGAGTATTTAGCCTTAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGTTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTC  >  NC_007793/2000171‑2000310

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