Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,590,050	C→T	T186T (ACG→ACA)	SAUSA300_RS07850 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,590,050	0	C	T	100.0%	62.4 / NA	24	T186T (ACG→ACA)	SAUSA300_RS07850	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/11);  total (13/11)

TTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAGTATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGCGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACTTCTAGTAAATCGAAATCGCTACCAGCTG  >  NC_007793/1589955‑1590144                                                                                                |                                                                                               ttAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGTGTAt                                                                                            <  1:5641317/100‑1 (MQ=255) ttAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGTGTAt                                                                                            <  1:2864464/100‑1 (MQ=255)                      aaaGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGt                                                                       >  1:5712994/1‑100 (MQ=255)                       aaGCTAAGAAAAAAGATTAATATTCATTCATTAAATATTAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGtt                                                                      >  1:13827/1‑100 (MQ=255)                           tAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGag                                                                  <  1:1439558/100‑2 (MQ=255)                              gAAAAAGGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGt                                                               <  1:5050645/100‑1 (MQ=255)                                   aaGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATc                                                          <  1:459231/100‑1 (MQ=255)                                                      tAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGTGTATGTTTAATATACAATGCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGt                                       <  1:1105772/100‑1 (MQ=255)                                                           ataAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTAAAATCTTGTCGCTc                                  >  1:4814580/1‑100 (MQ=255)                                                               aTCCAATTTAATTTGTTGTTTAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAAcc                              >  1:3977292/1‑99 (MQ=255)                                                                 ccAATTTAATTTGTTGTTTAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCGGATACTTTGAAATCTTGTCGCTCagcctc                            >  1:2792277/1‑95 (MQ=255)                                                                     tttaatttGTTGTTTAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATTCTTTGAAATCTTGTCGCTCAACCTCTAGt                        <  1:5239522/100‑1 (MQ=255)                                                                     tttaatttGTTGTTTAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACCTCTAGt                        <  1:957058/100‑1 (MQ=255)                                                                        aatttGTTGTTTAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAAccccaagatcg                     >  1:174560/1‑90 (MQ=255)                                                                           ttgttgTTTAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACCTCTAGTaaatcg                  >  1:2567347/1‑100 (MQ=255)                                                                                  ctAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACCTCTAGTAAATCGAAATCGc           >  1:4498920/2‑100 (MQ=255)                                                                                  accaGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACCCCTAGGagatcggaagagc           >  1:2048078/4‑96 (MQ=255)                                                                                   tAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACCTCTAGTAAATCGAAATCGCt          >  1:3510625/1‑100 (MQ=255)                                                                                   tAAGGTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACCTCTAGTAAATCGAAATCGCt          >  1:2284352/1‑100 (MQ=255)                                                                                      ggTCTACAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACCTCTAGTAAATCGAAATCGCTAcc       <  1:5344864/100‑1 (MQ=255)                                                                                           aCAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACCTCTAGTAAATCGAAATCGCTACCAGCTg  <  1:5037247/100‑1 (MQ=255)                                                                                           aCAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACCTCTAGTAAATCGAAATCGCTACCAGCTg  >  1:2107007/1‑100 (MQ=255)                                                                                           aCAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACCTCTAGTAAATCGAAATCGCTACCAGCTg  <  1:1529613/100‑1 (MQ=255)                                                                                           aCAAGTGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACCTCTAGTAAATCGAAATCGCCACCAgatc  >  1:3622139/1‑97 (MQ=255)                                                                                                |                                                                                               TTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAGTATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGCGTATGTTTAATATACAATTCATCGTTTGACGGTAAATCAGATACTTTGAAATCTTGTCGCTCAACTTCTAGTAAATCGAAATCGCTACCAGCTG  >  NC_007793/1589955‑1590144

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