Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,589,996	G→A	*42* (TAG→TAA)	SAUSA300_RS07845 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,589,996	0	G	A	100.0%	71.0 / NA	26	*42* (TAG→TAA)	SAUSA300_RS07845	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/16);  total (10/16)

ACTCGTAGTAACCCATTAAACAAAAGTTTTGAAAAGTCAGGCGCCAGCGAAAAATTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAGTATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGCGTATGTTTAATATACAATTCATCGTTTGACGGTAAATC  >  NC_007793/1589901‑1590088                                                                                                |                                                                                             actCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATatt                                                                                          >  1:2700871/1‑100 (MQ=255)  ctCGTAGTAACCCAATAAACAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATattc                                                                                         <  1:1716201/100‑1 (MQ=255)                  aaaCAAAAGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATaaa                                                                         >  1:3159360/1‑100 (MQ=255)                        aaGTTTTGAAAAGCCGGGTGCCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCaag                                                                   <  1:1511918/100‑2 (MQ=255)                                 gtgtgcgggtgcCAGCGAAAACTTAAAAAGCACTTTATCAGAAATAGCTAAGCAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATttgt                                                          >  1:4814450/11‑100 (MQ=255)                                     ccgggtgcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTgttgtt                                                      <  1:2796794/94‑1 (MQ=255)                                     ccgggtgcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTgttgtt                                                      >  1:4454824/7‑100 (MQ=255)                                      cgggtgcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGttt                                                     <  1:2144643/95‑1 (MQ=255)                                       gggtgtcAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTa                                                    <  1:1604523/94‑1 (MQ=255)                                       gggtgcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTa                                                    <  1:5193890/96‑1 (MQ=255)                                           gcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGTTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGt                                                >  1:576476/1‑100 (MQ=255)                                           gcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGt                                                <  1:4842213/100‑1 (MQ=255)                                           gcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGt                                                >  1:4010836/1‑100 (MQ=255)                                           gcCAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGt                                                >  1:4050606/1‑100 (MQ=255)                                            ccAGCGAAAACTTAAAAAGCACTTTAGCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTc                                               >  1:6050457/1‑100 (MQ=255)                                             cAGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCt                                              <  1:566346/100‑1 (MQ=255)                                              aGCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTa                                             <  1:3789678/100‑1 (MQ=255)                                               gCGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTAc                                            <  1:2217506/100‑1 (MQ=255)                                                cGAAAACTTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCtgaa                                           <  1:1576743/100‑4 (MQ=255)                                                       ttAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGTGTAt                                    <  1:2864464/100‑1 (MQ=255)                                                       ttAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAATATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGTGTAt                                    <  1:5641317/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)                                                                                                |                                                                                             ACTCGTAGTAACCCATTAAACAAAAGTTTTGAAAAGTCAGGCGCCAGCGAAAAATTAAAAAGCACTTTATCAGAAAAAGCTAAGAAAAAAGATTAGTATTCATTCATTAAATATAAATCCAATTTAATTTGTTGTTTAAGGTCTACAAGCGTATGTTTAATATACAATTCATCGTTTGACGGTAAATC  >  NC_007793/1589901‑1590088

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