Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,589,916	T→A	L16I (TTA→ATA)	SAUSA300_RS07845 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,589,916	0	T	A	100.0%	169.6 / NA	60	L16I (TTA→ATA)	SAUSA300_RS07845	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (52/8);  total (52/8)

AATCAAAAAACGCCTACTAGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGTAGATAAAAACAAAAAACAAGAAACTACTCGTAGTAACCCATTAAACAAAAGTTTTGAAAAGTCAGGCGCCAGCGAAAAATTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAGTATTC  >  NC_007793/1589822‑1590001                                                                                               |                                                                                      aaTCAAAAAACGCCTACAAGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAAc                                                                                  >  1:3584357/1‑100 (MQ=25)               tacaaGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGaaaa                                                                    >  1:2046978/1‑100 (MQ=38)                acaaGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACCAAAGTTTTGAAAAg                                                                   >  1:4948690/4‑100 (MQ=38)                 caaGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAgc                                                                  >  1:3996684/3‑99 (MQ=38)                 aaaGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAAAAAAAAACAAGAAGCTACTCGTAGGAACCCAATAAACAAAAGTTTTGAAAAgc                                                                  >  1:298778/3‑99 (MQ=38)                     tgtAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAAAAAAAGTTTTGCAAAgccggg                                                              >  1:3375461/1‑95 (MQ=255)                     cgtAGACGTTGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAgccggg                                                              >  1:4670977/2‑95 (MQ=255)                           tccgatctggtggTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAg                                                        <  1:4220980/93‑1 (MQ=255)                             ctGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCg                                                      >  1:3347942/2‑100 (MQ=255)                             ccgatctgtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTGGAAAAGCCGGGTGCCAGCg                                                      <  1:3142452/93‑1 (MQ=255)                              tGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGa                                                     >  1:1680375/1‑100 (MQ=255)                              tGAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAAACCAATAAACAAAAGTTTTGAAAAGCCGGGCGCCAGCGa                                                     >  1:1104808/1‑100 (MQ=255)                               gAATGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaa                                                    >  1:5777351/1‑100 (MQ=255)                                tctggtggTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaa                                                   >  1:3452206/3‑100 (MQ=255)                                taTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaa                                                   >  1:5444817/2‑100 (MQ=255)                                gctggtggTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaa                                                   >  1:3312694/3‑100 (MQ=255)                                gaTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTATTTGTAGCAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaag                                                   >  1:953029/2‑99 (MQ=255)                                aaTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaa                                                   >  1:828730/1‑100 (MQ=255)                                aaTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGGAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaa                                                   >  1:3914006/1‑100 (MQ=255)                                 ttggtggTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAGAGCCGGGTGCCAGCGaaaa                                                  >  1:2617275/2‑100 (MQ=255)                                 aTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaaa                                                  >  1:4588470/1‑100 (MQ=255)                                 aTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaaa                                                  >  1:4740688/1‑100 (MQ=255)                                 aTGGTGGTGAGAATTTTATGGCGGATAAAAACAAAAAACAAAAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGaaaa                                                  >  1:626814/1‑100 (MQ=255)                                  gtgtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaaac                                                 >  1:2295226/3‑99 (MQ=255)                                  gcctggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGaaaac                                                 >  1:2287713/4‑99 (MQ=255)                                  ctgtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAACTACTCGTAATAACCCAATAAACAAAAGCTTTGAAAAACCGGGGGGCAGCGaagat                                                 >  1:445744/3‑97 (MQ=255)                                  acatggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGGAACTACTTGTAGTAAACCAATAAACAAAAGGGTTGAAAAGCCGGGGGCCAGCGgaaac                                                 >  1:6321703/4‑99 (MQ=255)                                   ggtggtGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTTAAAAGCCGGGTGCCAGCGaaaact                                                >  1:5815377/1‑98 (MQ=255)                                   ggtggtGAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGCAGTAACCCAATAAACAAAAGCTTTGAAAAGCCGGGTGCCAGCGaaaact                                                >  1:5568292/1‑98 (MQ=255)                                    gtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGGAGCCACTCGTAGTAATCCAATAAACAAAAGCTTTGAAAAGCCCGGCGCCAgatcggaaga                                               >  1:496662/1‑91 (MQ=255)                                    gtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAAACCAATAAACAAAAGCTTTGAAAAGCCGGGTGCCAGCGAAAACtt                                               >  1:2350372/1‑100 (MQ=255)                                    gtggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGATACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGAGCCAGCGAAAACtt                                               >  1:4786197/1‑100 (MQ=255)                                      ggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGAAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGAAAACTTaa                                             >  1:5107228/1‑100 (MQ=255)                                      ggtgAGAATTTTATGGCGGATAAAAACAAAAAACAAGAACCTCCTCGTTGTAACCCCATTAACAAAAGTTTTGAAAAGCCGGGAGCCAGCGAAAACTTaa                                             >  1:3111504/1‑100 (MQ=255)                                                ttATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGGGCCAGCGAAAACTTAAAAAGCACttt                                   >  1:3428327/1‑100 (MQ=255)                                                ctATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACttt                                   >  1:3021768/2‑100 (MQ=255)                                                 tgtggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTa                                  >  1:1009260/8‑100 (MQ=255)                                                 tATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTa                                  >  1:4147905/1‑100 (MQ=255)                                                 tATGGCGGATAAAAACAAAAAACAAGAAGCTACTCGCAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTa                                  >  1:4847659/1‑100 (MQ=255)                                                  gtggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                 >  1:1486363/7‑100 (MQ=255)                                                  ctggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                 >  1:4542906/7‑100 (MQ=255)                                                  ctggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                 >  1:3590381/7‑100 (MQ=255)                                                  atggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                 >  1:6167256/7‑100 (MQ=255)                                                  atggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                 >  1:847219/7‑100 (MQ=255)                                                  atggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTAt                                 >  1:5030204/7‑100 (MQ=255)                                                   tggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc                                >  1:6307574/6‑100 (MQ=255)                                                   tggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc                                <  1:2546694/95‑1 (MQ=255)                                                   tggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGATTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc                                >  1:5653315/6‑100 (MQ=255)                                                   tGGCAGATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATc                                <  1:1395034/100‑1 (MQ=255)                                                    ggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCa                               >  1:227284/5‑100 (MQ=255)                                                    ggcggATAAAAACAAAAAACAAGAAGCTACTCGTAGTAACCCAATAAACAAAAGCTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCa                               >  1:141162/5‑100 (MQ=255)                                                              aacaaaaaacaaGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTa                     <  1:6159713/100‑1 (MQ=255)                                                              aacaaaaaacaaGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTCTATCAGAAAAAGCTa                     <  1:4160601/100‑1 (MQ=255)                                                                   aaaacaaGAAGCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAGAGCTAAGatc                >  1:5303032/1‑98 (MQ=255)                                                                        aagaagCTACTCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGa           >  1:365721/1‑100 (MQ=255)                                                                           tccgatctCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTa        <  1:1631534/94‑1 (MQ=255)                                                                             gctactCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTaat      >  1:2743385/2‑98 (MQ=255)                                                                             agtactCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTaat      >  1:5917729/3‑98 (MQ=255)                                                                                actCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATatt   >  1:2700871/1‑100 (MQ=255)                                                                                 ctCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATattc  <  1:1716201/100‑1 (MQ=255)                                                                                               |                                                                                      AATCAAAAAACGCCTACTAGTGTAGACGTTGAATGGTGGTGAGAATTTTATGGTAGATAAAAACAAAAAACAAGAAACTACTCGTAGTAACCCATTAAACAAAAGTTTTGAAAAGTCAGGCGCCAGCGAAAAATTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAGTATTC  >  NC_007793/1589822‑1590001

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