Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,590,197	T→A	I137I (ATA→ATT)	SAUSA300_RS07850 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,590,197	0	T	A	100.0%	196.3 / NA	70	I137I (ATA→ATT)	SAUSA300_RS07850	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (32/38);  total (32/38)

CAACTTCTAGTAAATCGAAATCGCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTTATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAATTTCATAAACTTTCTGGTTAT  >  NC_007793/1590112‑1590292                                                                                      |                                                                                                cAACCTCTAGTAAATCGAAATCGCTACCAGCTGAATTATAGGTTTTAAGTTCAGCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATAc                                                                                   >  1:1618829/1‑100 (MQ=255)   cgatCTAGT‑AATCGAAATCGCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTgg                                                                                <  1:2015704/97‑1 (MQ=255)       ctggcaaatcgaaatcgCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGttt                                                                             <  1:1759375/95‑1 (MQ=255)       cTAGTGGATCGAAATCGCTGCCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGttt                                                                             <  1:1820548/100‑1 (MQ=255)       cTAGTAAATCGAAATCGCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGttt                                                                             >  1:1587458/1‑100 (MQ=255)       cTAGTAAATCGAAAGCGCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGttt                                                                             <  1:2217864/100‑1 (MQ=255)       cTAGCAAATCGAAATCGCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAACTACTGGttt                                                                             >  1:4178658/1‑100 (MQ=255)         aGTAAATCGAAATCGCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTAc                                                                           >  1:2332984/1‑100 (MQ=255)         aGTAAATCGAAATCGCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTAc                                                                           >  1:5688123/1‑100 (MQ=255)         aGTAAATCGAAATCGCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTAc                                                                           >  1:6577050/1‑100 (MQ=255)         aGTAAATCGAAATCGCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTATGTTTTCGAAGTCTTTAATAACTATAAATACTGGttaga                                                                           >  1:2461745/1‑97 (MQ=255)               tcgaaatcgCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCgttgt                                                                     <  1:3638405/100‑1 (MQ=255)                 gaaatcgCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGtta                                                                   <  1:4345202/100‑1 (MQ=255)                      cgCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTaaa                                                              >  1:3058771/1‑100 (MQ=255)                      cgCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTaaa                                                              >  1:1294503/1‑100 (MQ=255)                         tACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAaacaa                                                           >  1:471883/1‑100 (MQ=255)                          gccAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAaacaac                                                          >  1:3562918/2‑100 (MQ=255)                          aCCAGCTGAATTATAGGTTTTATGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAaacaac                                                          >  1:6637331/1‑100 (MQ=255)                          aCCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAaacaac                                                          >  1:3243529/1‑100 (MQ=255)                          aCCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAaacaac                                                          >  1:2679318/1‑100 (MQ=255)                          aCCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAaacaac                                                          >  1:5666154/1‑100 (MQ=255)                          aCCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAaacaac                                                          >  1:5142399/1‑100 (MQ=255)                          aCCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAaacaac                                                          >  1:1751400/1‑100 (MQ=255)                          aCCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAaacaac                                                          >  1:1312946/1‑100 (MQ=255)                             cgCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTg                                                       >  1:4650812/2‑100 (MQ=255)                             cgCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTg                                                       >  1:5033159/2‑100 (MQ=255)                              gCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGt                                                      >  1:1489819/1‑100 (MQ=255)                                     tataGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCtttt                                               >  1:863531/1‑100 (MQ=255)                                       taGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGt                                             >  1:6059824/1‑100 (MQ=255)                                         ggTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCt                                           <  1:4638032/100‑1 (MQ=255)                                         ggTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCt                                           <  1:4317655/100‑1 (MQ=255)                                         ggTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCt                                           <  1:2337150/100‑1 (MQ=255)                                         ggTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCt                                           <  1:2283675/100‑1 (MQ=255)                                         ggTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCt                                           <  1:3004134/100‑1 (MQ=255)                                         ggTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCt                                           <  1:2215236/100‑1 (MQ=255)                                          gTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCta                                          <  1:5640449/100‑1 (MQ=255)                                            tttAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACTGTTGTTATTAAACAACTTGTCTCTTTTGTCtaat                                        <  1:4065244/100‑1 (MQ=255)                                            tttAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCtagg                                        <  1:733647/100‑3 (MQ=255)                                            tttAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCtaat                                        <  1:433446/100‑1 (MQ=255)                                            tttAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCtaat                                        <  1:3465434/100‑1 (MQ=255)                                            tttAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCtaat                                        <  1:4430891/100‑1 (MQ=255)                                            tttAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCtaat                                        <  1:4469270/100‑1 (MQ=255)                                            tttAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCtaag                                        <  1:6288786/100‑2 (MQ=255)                                             ttAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCtaata                                       <  1:4406178/100‑1 (MQ=255)                                             ttAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCtaata                                       <  1:4593013/100‑1 (MQ=255)                                             ttAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCtaata                                       <  1:5386685/100‑1 (MQ=255)                                             ttAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCtaata                                       <  1:2370743/100‑1 (MQ=255)                                               aaGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAg                                     >  1:5004088/1‑100 (MQ=255)                                                  ttCCCCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATGCTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCtt                                  <  1:3000275/100‑1 (MQ=255)                                                  ttCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCtt                                  <  1:5042159/100‑1 (MQ=255)                                                  tccgatctctTCGATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCtt                                  <  1:4396353/94‑1 (MQ=255)                                                  ccgaCCCTCTTCAATGATTCTGTTTTCACAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCtt                                  >  1:2065595/4‑100 (MQ=255)                                                  atctcccTTTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGcat                                  <  1:6313446/96‑3 (MQ=255)                                                    tcgccTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTtat                                <  1:4244069/97‑1 (MQ=255)                                                    cACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATTCTGGTTTACCGTTGTTATTAAACAACTTGTTTCTTTTGTCTAATAAGCTtat                                >  1:1533480/1‑100 (MQ=255)                                                    cACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTtat                                <  1:4403436/100‑1 (MQ=255)                                                    cACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTtag                                <  1:6422831/100‑2 (MQ=255)                                                    cACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTATCTAATAAGCTtat                                >  1:4998168/1‑100 (MQ=255)                                                     aCCCTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTtata                               >  1:5614957/1‑100 (MQ=255)                                                      cccTCTTCAATGATTCTGTTTTCAAATTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATAc                              >  1:2561560/1‑100 (MQ=255)                                                      cccTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATAc                              <  1:2320798/100‑1 (MQ=255)                                                      accTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATAc                              >  1:5374467/2‑100 (MQ=255)                                                       ccTCTTCAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACa                             <  1:5631356/100‑1 (MQ=255)                                                             cAATGATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCaa                       >  1:3596080/1‑100 (MQ=255)                                                                 gATTCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTc                   >  1:6477392/1‑100 (MQ=255)                                                                   ttCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCAt                 <  1:459804/100‑1 (MQ=255)                                                                   ttCTGTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCAt                 <  1:1578514/100‑1 (MQ=255)                                                                       gTTTTCAAAGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAAc             <  1:1311993/100‑1 (MQ=255)                                                                             aaaGTCTTTAATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCtt       <  1:2803566/100‑2 (MQ=255)                                                                                  cTTTAATAACTATAAATACTGGTTTACCGTTGTTACTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAAATTCATAAACTTTCTTGtttc  <  1:4373173/100‑3 (MQ=255)                                                                                      |                                                                                                CAACTTCTAGTAAATCGAAATCGCTACCAGCTGAATTATAGGTTTTAAGTTCACCCTCTTCAATGATTCTGTTTTCAAAGTCTTTTATAACTATAAATACTGGTTTACCGTTGTTATTAAACAACTTGTCTCTTTTGTCTAATAAGCTTATACAATCCAATTTCATAAACTTTCTGGTTAT  >  NC_007793/1590112‑1590292

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