Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,719,734	G→A	E8K (GAA→AAA)	SAUSA300_RS13970 →	amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,719,734	0	G	A	100.0%	136.0 / NA	52	E8K (GAA→AAA)	SAUSA300_RS13970	amidohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (51/1);  total (51/1)

CTGCAATAGATGCAATAAACATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAGAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATCAG‑CAGATCCTAAAGGTGTA  >  NC_007793/2719666‑2719802                                                                     |                                                                      cTGCAATAGATGCAATAAACATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGa                                        >  1:2342192/1‑100 (MQ=255)   gCAATAGATGCAATAAACATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGaaa                                      >  1:1859164/1‑100 (MQ=255)   gCAATAGATGCAATAAACATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATACTCAAAAAGaaa                                      >  1:3374343/1‑100 (MQ=255)   ccAATAGATGCAATAAACATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATAAAAAAAGatc                                      >  1:4707986/2‑98 (MQ=255)      aTAGATGCAATAAACATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGaagatc                                   >  1:4869356/1‑96 (MQ=255)           tGCAATAAACATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAAt                              >  1:2317989/1‑100 (MQ=255)              aaTAAACATGCAATACAATATGGAGGCGAAGTAAATGTAAAGTATTACGTTCGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCg                           >  1:1417271/1‑100 (MQ=255)               aTAAACATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGa                          >  1:859206/1‑100 (MQ=255)                gccaCATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGat                         >  1:1447761/4‑100 (MQ=255)                 aaaCATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGata                        >  1:2578218/1‑100 (MQ=255)                  taCATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGatat                       >  1:5656116/2‑100 (MQ=255)                  gccATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGatat                       >  1:6280869/3‑100 (MQ=255)                   gtaTGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATc                      >  1:2316289/3‑100 (MQ=255)                   aCATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGatacc                      >  1:6645332/1‑98 (MQ=255)                   aCATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATc                      <  1:2784334/100‑1 (MQ=255)                    cctGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATAAca                     >  1:4746207/3‑100 (MQ=255)                    cATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGatatta                     >  1:1674108/1‑98 (MQ=255)                     gtGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATACcag                    >  1:2392333/2‑100 (MQ=255)                     ctGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATGTcag                    >  1:5927871/2‑100 (MQ=255)                     ctGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATCTcag                    >  1:16986/2‑100 (MQ=255)                     ctGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATcag                    >  1:4945612/2‑100 (MQ=255)                     aTGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATcag                    >  1:1102698/1‑100 (MQ=255)                       gCAATACAATATGGAGGCGAAGTAAATGAAAGGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATAgcag‑ca                 >  1:6641226/1‑100 (MQ=255)                       gCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTGCGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATcag‑ca                 >  1:3000767/1‑100 (MQ=255)                       gCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATGCcag‑ca                 >  1:2380674/1‑100 (MQ=255)                       gCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATcagtc                  >  1:6208143/1‑98 (MQ=255)                       gCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATcag‑ca                 >  1:143171/1‑100 (MQ=255)                       gCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATAGcagat                  >  1:2097463/1‑98 (MQ=255)                       gCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAACGCGAGATcggaa                  >  1:4038703/1‑96 (MQ=255)                       gCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATCGAAGATATTCAAAAAGAAACGATGAATGCGATATcag‑ca                 >  1:3502089/1‑100 (MQ=255)                        ctatacaataTGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATGTcag‑cag                >  1:679981/3‑100 (MQ=255)                         tatacaataTGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATAACAG‑CAGa               >  1:3195914/2‑100 (MQ=255)                         aatacaataTGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATCAG‑CAGa               >  1:1618072/1‑100 (MQ=255)                          atacaataTGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATCAG‑CAGAt              >  1:5427358/1‑100 (MQ=255)                           tacaataTGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATACCAG‑CAGATc             >  1:896187/1‑100 (MQ=255)                            acaataTGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATGTCAG‑CAGATcg            >  1:997351/1‑99 (MQ=255)                            acaataTGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATGTCAG‑CAGATcg            >  1:1336369/1‑99 (MQ=255)                            acaataTGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATCAG‑CAGATcc            >  1:3074280/1‑100 (MQ=255)                            acaataTGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATACCAG‑CAGATcg            >  1:1129596/1‑99 (MQ=255)                              aataTGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATCAG‑CAGAttgga          >  1:5487853/1‑96 (MQ=255)                               ctttGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATGCCAG‑CAGATcggaa         >  1:454325/4‑96 (MQ=255)                               ctatGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATACCAG‑CAGATcggaa         >  1:241748/2‑96 (MQ=255)                                 atGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATGTCAG‑CAGATcggaaga       >  1:544736/1‑94 (MQ=255)                                 atGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATGCCAG‑CAGATcggaaga       >  1:5136209/1‑94 (MQ=255)                                 atGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATTAG‑CAGATcggaaga       >  1:1521809/1‑94 (MQ=255)                                 atGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATACCAG‑CAGATcggaaga       >  1:1575071/1‑94 (MQ=255)                                 acggAGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATAAG‑CAGATcggaaga       >  1:2136424/3‑94 (MQ=255)                                  tccaGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATGTCAG‑CAGATcggaagag      >  1:2184425/4‑93 (MQ=255)                                   tcaGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGCTGAATGCGATACAAG‑TAGATcggaagagc     >  1:1134561/3‑92 (MQ=255)                                     aGGCGAAGTAAATGAAAAGTATTACGTTTGAAAAGCATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATCAA‑CAGATcggaagagcac   >  1:2880180/1‑90 (MQ=255)                                      ggCGAAGTAAATGAAAAGTATTACGTTTGAAAAACATTATGTCATTGAAGATATTCGAAAAGAAACGATGAATGCGATATCAG‑CAGATCCTAAAGGTGTa  >  1:2079420/1‑100 (MQ=255)                                      ggCGAAGTAAATGAAAAGCATTACGTTTGAAAAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATCAGTCAGATcggaagagcac   >  1:1003751/1‑90 (MQ=255)                                                                     |                                                                      CTGCAATAGATGCAATAAACATGCAATACAATATGGAGGCGAAGTAAATGAAAAGTATTACGTTTGAAGAACATTATGTCATTGAAGATATTCAAAAAGAAACGATGAATGCGATATCAG‑CAGATCCTAAAGGTGTA  >  NC_007793/2719666‑2719802

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