Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	877,476	A→G	K241E (AAA→GAA)	SAUSA300_RS04280 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	877,476	0	A	G	100.0%	97.6 / NA	36	K241E (AAA→GAA)	SAUSA300_RS04280	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/13);  total (23/13)

GTAACACATTTACTGCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGAAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATGGGGTGATTTAGGTATTGCCAATATCTT  >  NC_007793/877384‑877567                                                                                             |                                                                                            gtAACACATTTACTGCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTa                                                                                      >  1:2699645/1‑100 (MQ=255)   aaCACATTTACTGCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATg                                                                                    >  1:5533867/1‑100 (MQ=255)     cacaTTTACTGCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAATTAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGAt                                                                                  <  1:4817650/100‑1 (MQ=255)       caTTTACTGCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATag                                                                                <  1:1282969/100‑1 (MQ=255)          ttACTGCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGt                                                                             >  1:4769314/1‑100 (MQ=255)           tACTGCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTc                                                                            >  1:631233/1‑100 (MQ=255)             cTGCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCtt                                                                          <  1:4704860/100‑1 (MQ=255)              ggCTAAGGAACCATTTAATCAACTTGGAGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTa                                                                         >  1:2693805/2‑100 (MQ=255)               gCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCGAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTAc                                                                        >  1:902832/1‑100 (MQ=255)               gCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTAc                                                                        >  1:589892/1‑100 (MQ=255)               gCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTAc                                                                        >  1:1808532/1‑100 (MQ=255)               gCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGCCTTAc                                                                        >  1:400049/1‑100 (MQ=255)                  aaGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGc                                                                     <  1:5216714/100‑1 (MQ=255)                  aaGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGc                                                                     <  1:296795/100‑1 (MQ=255)                              aaTCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCACGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTa                                                         >  1:353081/1‑100 (MQ=255)                                  aaCTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTGCAag                                                     <  1:2745540/100‑1 (MQ=255)                                  aaCTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAag                                                     <  1:4538542/100‑1 (MQ=255)                                       ggtggtTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATg                                                >  1:1297988/1‑100 (MQ=255)                                         tggtTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGt                                              <  1:1099570/100‑1 (MQ=255)                                          ggtTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTa                                             >  1:2850304/1‑100 (MQ=255)                                             tATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAAtt                                          >  1:5579912/1‑100 (MQ=255)                                             tATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAAtt                                          >  1:5181019/1‑100 (MQ=255)                                                          aCCAAATAGGGCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATg                             >  1:4779286/1‑100 (MQ=255)                                                             aaaTAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATgggg                          <  1:3517022/100‑1 (MQ=255)                                                               aTAGATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATGGGGta                        <  1:405816/100‑2 (MQ=255)                                                                gggATCCAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAATAGATGGTAATTCTTCGATTATATggggaag                       >  1:4039564/3‑97 (MQ=255)                                                                     ccAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATggggcagatcgg                  >  1:1860498/1‑92 (MQ=255)                                                                     ccAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATggggaaagatcg                  >  1:5112066/1‑92 (MQ=255)                                                                      cAAGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATGGGGTGATTTAgg                 >  1:1353269/1‑100 (MQ=255)                                                                       aaGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATGGGGTGATTTAGGt                >  1:4957893/1‑100 (MQ=255)                                                                       aaGAACGAACGATCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATGGGGAGAtcggaag                >  1:2229924/1‑94 (MQ=255)                                                                              ttccgatctgAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATGGGGTGATTTAGGTATTGCCa         <  1:3029453/91‑1 (MQ=255)                                                                              aacgaTCAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATGGGGTGATTTAGGTATTGCCa         <  1:1535970/100‑1 (MQ=255)                                                                                    cAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATGGGGTGATTTAGGTATTGCCAATATCt   >  1:5137154/1‑100 (MQ=255)                                                                                    cAAGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATGGGGTGATTTAGGTATTGCCAATATCt   >  1:440462/1‑100 (MQ=255)                                                                                     aaGAACTGGAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATGGGGTGATTTAGGTATTGCCAATATCtt  <  1:700004/100‑1 (MQ=255)                                                                                             |                                                                                            GTAACACATTTACTGCTAAGGAACCATTTAATCAACTTGGTGGTTATCCATATTTTGACCAAATAGATCCAAGAACGAACGATCAAGAACTGAAAATGTATGATAGAGTCTTACTGCAAATTGATTCTACAAGAGATGGTAATTCTTCGATTATATGGGGTGATTTAGGTATTGCCAATATCTT  >  NC_007793/877384‑877567

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: ‑