Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,578,251	A→C	L50V (TTG→GTG)	SAUSA300_RS13260 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,578,251	0	A	C	100.0%	92.7 / NA	36	L50V (TTG→GTG)	SAUSA300_RS13260	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (1/35);  total (1/35)

TGTATGATATTTATCCTCCAAAATGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCAACGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCATTAAAGATTAACGCTCTCAAATATTTATATTCTTCTTGCCATTG  >  NC_007793/2578156‑2578344                                                                                                |                                                                                              tgtATGATATTTATCCTCCAAAATGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACA‑TTTTATTATTCACCgtgt                                                                                          <  1:5168473/100‑1 (MQ=255)  gtATGATATTTATCCTCCAAAATGGCTCCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTCCTACATTTTTATTATTCACCgtgt                                                                                          <  1:1854204/100‑1 (MQ=255)  gtATGATATTTATCCTCCAAAATGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCgtgt                                                                                          <  1:2243424/100‑1 (MQ=255)   tATGATATTTATCCTCCAAAATGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTa                                                                                         <  1:4082687/100‑1 (MQ=255)   tATGATATTTATCCTCCAAAATGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTa                                                                                         <  1:1655509/100‑1 (MQ=255)      gATATTTATCCTCCAAAATGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACa                                                                                      <  1:888633/100‑1 (MQ=255)      gATATTTATCCTCCAAAATGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACa                                                                                      <  1:1850948/100‑1 (MQ=255)       atatTTATCCTCCAAAATGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACaa                                                                                     <  1:2563316/100‑1 (MQ=255)        tatTTATCCTCCAAAATGGCACCTTTATGAAATAACAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAg                                                                                    <  1:1928057/100‑1 (MQ=255)              tcctccAAAATGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATgtt                                                                              <  1:5514827/100‑3 (MQ=255)              tcctccAAAATGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATgag                                                                              <  1:3338591/100‑3 (MQ=255)                ctccAAAATGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATc                                                                            <  1:4691644/100‑1 (MQ=255)                    aaaaTGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCAtt                                                                        <  1:3608363/100‑1 (MQ=255)                                 tATGAAATAATAGTTCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCctc                                                           >  1:1555232/1‑100 (MQ=255)                                  aTGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCctct                                                          <  1:5362923/100‑1 (MQ=255)                                         aataGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTg                                                   <  1:1369141/100‑1 (MQ=255)                                         aataGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTg                                                   <  1:1597018/100‑1 (MQ=255)                                          ataGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTga                                                  <  1:856649/100‑2 (MQ=255)                                          ataGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGt                                                  <  1:5701037/100‑1 (MQ=255)                                          ataGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGt                                                  <  1:1047741/100‑1 (MQ=255)                                          ataGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGt                                                  <  1:4265772/100‑1 (MQ=255)                                            aGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTggca                                                <  1:5691616/100‑4 (MQ=255)                                             gTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAGTTCTGTTTc                                               <  1:2857706/100‑1 (MQ=255)                                             gTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTc                                               <  1:4252463/100‑1 (MQ=255)                                                  aCATAATTTTTGAAGCCATGGATAAGTACTGCATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTTTTTCATTaa                                          <  1:3717130/100‑1 (MQ=255)                                                  aCATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCATTaa                                          <  1:4919702/100‑1 (MQ=255)                                                  aCATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCATTaa                                          <  1:3285573/100‑1 (MQ=255)                                                  aCATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCATTaa                                          <  1:3057696/100‑1 (MQ=255)                                                  aCATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCATTaa                                          <  1:5731414/100‑1 (MQ=255)                                                  aCATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCATTaa                                          <  1:710617/100‑1 (MQ=255)                                                  aCATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCATTaa                                          <  1:2122791/100‑1 (MQ=255)                                                     tAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCAttgcttt                                       <  1:3449606/100‑6 (MQ=255)                                                     tAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCAttgcttt                                       <  1:3000628/100‑6 (MQ=255)                                                                  cATGGATAAGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCATTAAAGATTAACGCTCTCaa                          <  1:1476161/100‑1 (MQ=255)                                                                          aGTACTACATTTTTATTATTCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCATTAAAGATTAACGCTCTCAAATATTtat                  <  1:1726014/100‑1 (MQ=255)                                                                                          tattCACCGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCATTAAAGATTAACGCTCTCAAATATTTATATTCTTCTTGCCATTg  <  1:2203234/100‑1 (MQ=255)                                                                                                |                                                                                              TGTATGATATTTATCCTCCAAAATGGCACCTTTATGAAATAATAGTGCAACATAATTTTTGAAGCCATGGATAAGTACTACATTTTTATTATTCAACGTGTAACAAGGATGCATCCATTTATACGCTTCCTCTAATTCTGTTTCATTAAAGATTAACGCTCTCAAATATTTATATTCTTCTTGCCATTG  >  NC_007793/2578156‑2578344

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