Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	716,517	G→A	intergenic (+23/+142)	SAUSA300_RS03435 → / ← SAUSA300_RS03440	lipase LipA/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	716,517	0	G	A	100.0%	98.4 / NA	36	intergenic (+23/+142)	SAUSA300_RS03435/SAUSA300_RS03440	lipase LipA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/6);  total (30/6)

AGTGGTATTCAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTGCTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAA  >  NC_007793/716420‑716603                                                                                                  |                                                                                       aGTGGTATTCAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACt                                                                                      >  1:4124990/1‑100 (MQ=255) aGTGGTATTCAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACt                                                                                      >  1:2527012/1‑100 (MQ=255)      tATTCAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGc                                                                                 >  1:2803619/1‑100 (MQ=255)       aTTCAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACCGTAGTAATTTATTACTTAAGCa                                                                                >  1:2901649/1‑100 (MQ=255)        gccAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCaa                                                                               >  1:3136802/3‑100 (MQ=255)        cccAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCaa                                                                               >  1:2938614/3‑100 (MQ=255)         ttaaaCTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAAc                                                                              >  1:1351987/3‑100 (MQ=255)         tCAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAAc                                                                              <  1:2478748/100‑1 (MQ=255)         tCAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAAc                                                                              >  1:2695654/1‑100 (MQ=255)         tCAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAAc                                                                              <  1:4197851/100‑1 (MQ=255)          cAAACTGAAGAGAAACCGCAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACa                                                                             >  1:404064/1‑100 (MQ=255)          cAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCGACa                                                                             >  1:2359536/1‑100 (MQ=255)          cAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACa                                                                             >  1:2445954/1‑100 (MQ=255)          cAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACa                                                                             >  1:1437334/1‑100 (MQ=255)              cTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTa                                                                         <  1:1694045/100‑1 (MQ=255)                gAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAg                                                                       >  1:2863134/1‑100 (MQ=255)                gAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAg                                                                       >  1:3802207/1‑100 (MQ=255)                    agaAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGAttt                                                                   >  1:4562386/1‑100 (MQ=255)                     taaaCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGAtttt                                                                  >  1:5633780/2‑100 (MQ=255)                     gaAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGAtttt                                                                  >  1:1899198/1‑100 (MQ=255)                      aaaCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTc                                                                 >  1:148011/1‑100 (MQ=255)                      aaaCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTc                                                                 >  1:4205069/1‑100 (MQ=255)                      aaaCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTc                                                                 >  1:4273190/1‑100 (MQ=255)                       aaCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCa                                                                >  1:4649834/1‑100 (MQ=255)                        gccACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCATTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCaa                                                               >  1:3735733/2‑100 (MQ=255)                        aCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCaa                                                               >  1:4851781/1‑100 (MQ=255)                        aCCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCaa                                                               >  1:2326173/1‑100 (MQ=255)                         ccACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCaaa                                                              >  1:581326/1‑100 (MQ=255)                                                           tAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTaa                            >  1:123255/1‑100 (MQ=255)                                                            ccttCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAAt                           >  1:2872696/4‑100 (MQ=255)                                                            aaaTCCTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTagc                           <  1:5353264/100‑3 (MQ=255)                                                                ccTTTAAACTAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATgtgt                       <  1:3305884/100‑1 (MQ=255)                                                                         tAATGATAAACAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGGATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTg              <  1:5521008/100‑1 (MQ=255)                                                                                ctaCAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGatcgg       >  1:915334/3‑97 (MQ=255)                                                                                aaaCAGTAGTAATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGACTGAGatcgg       >  1:3745616/1‑97 (MQ=255)                                                                                     gtagtaATTTATTACTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATaaa  >  1:3252445/1‑100 (MQ=255)                                                                                                  |                                                                                       AGTGGTATTCAAACTGAAGAGAAACCACAAATAGAAAATCCGAGTAATGAATTACCGTTAAATCCTTTAAACTAATGATAAACAGTAGTAATTTATTGCTTAAGCAACATTTAAGATTTTCAAATTAAAAACGAAGAATTTAAAACATGTGGTGCTAATGTGTAAGAGTCTGAGATATAATAAA  >  NC_007793/716420‑716603

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