Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,110,894	C→T	M33I (ATG→ATA)	SAUSA300_RS10670 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,110,894	0	C	T	100.0%	61.1 / NA	24	M33I (ATG→ATA)	SAUSA300_RS10670	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (4/20);  total (4/20)

ATCGTTATTCATTCGCTTTAATTCAATCCTGTATTCTTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTCATATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTGATACAAACTAATAGCAATATAGCAAATACAATAGTCATAGTAATCCA  >  NC_007793/2110804‑2110989                                                                                           |                                                                                                atCGTTATTCATTCGCTTTAATTCAATCCTGTATTCTTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTg                                                                                        <  1:1342989/100‑1 (MQ=255)          cattcGCTTTAATTCAATCCTGTATTCTTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGt                                                                               <  1:3110993/100‑1 (MQ=255)                     tccgatctcTGTATTCTTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGc                                                                    <  1:2632011/92‑1 (MQ=255)                     aTTCAATCCTGTATTCTTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGc                                                                    <  1:4224847/100‑1 (MQ=255)                      ttCAATCCTGTATTCCTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCtgcgt                                                                   <  1:1209832/100‑5 (MQ=255)                      tccgaTCTTGTATTCTTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCtggtg                                                                   <  1:1875327/96‑5 (MQ=255)                      ccgatctatGTATTCTTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCa                                                                   <  1:2913284/92‑1 (MQ=255)                        ctttCCTGTATTCTTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACg                                                                 <  1:2272903/97‑1 (MQ=255)                         cgatcTGTATTCTTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACgg                                                                <  1:1452890/96‑2 (MQ=255)                          tctcTGTACTCTTTCAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGAt                                                               <  1:4938848/97‑1 (MQ=255)                          aTCCTGTATTCTTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGag                                                               <  1:801277/100‑2 (MQ=255)                             cTGTATTCTTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCagatcgga                                                            >  1:6069515/1‑93 (MQ=255)                                 tccgatctatcccGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTAtt                                                        <  1:1908823/90‑1 (MQ=255)                                        tAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCATGAAGATAATCATTTATATATCTGCGCGCTTGTACCTCCCTTGCACGATCACTATTAATACTg                                                 >  1:5061750/1‑100 (MQ=255)                                        tAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTg                                                 >  1:2755447/1‑100 (MQ=255)                                        tAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTg                                                 <  1:549568/100‑1 (MQ=255)                                                 gTATCCTTTAGTTTTAACTACTTCATCAAGTAGTTAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTGATACTAACt                                        >  1:386454/1‑100 (MQ=255)                                                         tAGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTGATACAAACTAATAGCaa                                <  1:3922034/100‑1 (MQ=255)                                                          aGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTGATACAAACTAATAGCAat                               <  1:2334033/100‑1 (MQ=255)                                                          aGTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTGATACAAACTAATAGCAat                               <  1:5641397/100‑1 (MQ=255)                                                           gTTTTAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTGATACAAACTAATAGCagct                              <  1:3849419/100‑4 (MQ=255)                                                               tAACTACTTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGATCTCTATTAATACTGATACAAACTAATAGCAATATAGc                          <  1:3681751/100‑1 (MQ=255)                                                                 actactTCATCAAGTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTGATACAAACTAATAGCAATATAGCaa                        <  1:6164914/100‑1 (MQ=255)                                                                              gTAGATAATCATTTATATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTGATACAAACTAATAGCAATATAGCAAATACAATAGTCAt           <  1:1844074/100‑1 (MQ=255)                                                                                      tcatttatATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTGATACAAACTAATAGCAATATAGCAAATACAATAGTCATAGTAATcc   <  1:3339293/100‑1 (MQ=255)                                                                                       catttatATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTGATACAAACTAATAGCAATATAGCAAATACAATAGTCATAGTAATCca  >  1:4718736/1‑100 (MQ=255)                                                                                           |                                                                                                ATCGTTATTCATTCGCTTTAATTCAATCCTGTATTCTTTTAACCCGTTGTATCCTTTAGTTTTAACTACTTCATCAAGTAGATAATCATTCATATATCTGAGCGCTTGTATCTCCCTTGCACGATCACTATTAATACTGATACAAACTAATAGCAATATAGCAAATACAATAGTCATAGTAATCCA  >  NC_007793/2110804‑2110989

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