Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	46,855	T→A	S222S (TCA→TCT)	SAUSA300_RS00200 ←	recombinase RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	46,855	0	T	A	100.0%	148.0 / NA	51	S222S (TCA→TCT)	SAUSA300_RS00200	recombinase RecB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (4/47);  total (4/47)

TTTTTCACTCCAATCTTTGTATTTCGCAAATTGAATTTTGCCGATATAGAAAGGGTTAGCTAAGATGTATGTGATTGAACTAATACTGAAAGGTTTCCCCTTTTTAGTGACATATCCTTTGTGATTCAATGCATTGGCAATTTTACGATAGCCATGTCCTTTGGCATAGC  >  NC_007793/46780‑46949                                                                            |                                                                                               tttttCACTCCAATCTTTGTATTTCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTaccc                                                                        <  1:4330489/100‑1 (MQ=255)             aTCTTTGTATTTCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGAct                                                            <  1:3527879/100‑2 (MQ=255)                 ttGTATTTCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACatagc                                                        >  1:2784622/1‑98 (MQ=255)                 ttGTATTTCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACatagc                                                        >  1:2417541/1‑98 (MQ=255)                   gTATTTCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCt                                                      <  1:3098719/100‑1 (MQ=255)                    tATTTCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGccca                                                     <  1:3563574/100‑3 (MQ=255)                    tATTTCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCtt                                                     <  1:2232301/100‑1 (MQ=255)                    tATTTCGCAAATTGAATATTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGTACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCtt                                                     <  1:1040591/100‑1 (MQ=255)                     aTTTCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATAGATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCtaa                                                    <  1:2682796/100‑3 (MQ=255)                      tttCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTCCCCTTTTTAGTGACATAGCCttta                                                   <  1:2269635/100‑2 (MQ=255)                      tttCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCttag                                                   <  1:3057996/100‑3 (MQ=255)                      tttCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTtg                                                   <  1:4552623/100‑1 (MQ=255)                      tttCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTtg                                                   <  1:2783447/100‑1 (MQ=255)                        tCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTtgtg                                                 <  1:1670947/100‑1 (MQ=255)                        tCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTtgtg                                                 <  1:3248582/100‑1 (MQ=255)                        tCGCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTtgtg                                                 <  1:1524195/100‑1 (MQ=255)                          gCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTtgtgta                                               <  1:365728/100‑3 (MQ=255)                          gCAAATTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGAt                                               <  1:3095528/100‑1 (MQ=255)                            aaaTTGAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTc                                             <  1:5437303/100‑1 (MQ=255)                                 gAATTTTGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGc                                        <  1:3454481/100‑1 (MQ=255)                                    ttttGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCaat                                     <  1:1319408/100‑3 (MQ=255)                                    ttttGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCAtg                                     <  1:3002794/100‑2 (MQ=255)                                       tGCCGATATAGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGc                                  <  1:3934308/100‑1 (MQ=255)                                             tataGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTAtttt                            >  1:3461329/1‑100 (MQ=255)                                                aGAATGGGTTAGCTAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACg                         >  1:1110348/1‑100 (MQ=255)                                                             tAATATATATGTGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCt            <  1:2606793/100‑1 (MQ=255)                                                               atatatatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCttt          <  1:3008391/95‑1 (MQ=255)                                                                  tatatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGc       <  1:610561/98‑1 (MQ=255)                                                                  tatatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGc       <  1:4482214/98‑1 (MQ=255)                                                                  tatatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGc       <  1:2030079/98‑1 (MQ=255)                                                                  tatatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGc       <  1:1258180/98‑1 (MQ=255)                                                                   ataggtgATAGAACTATTACT‑AAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt     <  1:1092827/96‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCag     <  1:5199615/100‑2 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt     <  1:2753798/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt     <  1:5660355/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt     <  1:3583442/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt     <  1:538974/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt     <  1:538591/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt     <  1:1776020/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt     <  1:5045592/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCCTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt     <  1:4797325/100‑1 (MQ=255)                                                                    tatgtGATAAAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt     <  1:4312406/100‑1 (MQ=255)                                                                     atgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCagt    <  1:1869065/100‑3 (MQ=255)                                                                     atgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAtt    <  1:2955066/100‑2 (MQ=255)                                                                     atgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAtt    <  1:3306724/100‑2 (MQ=255)                                                                     atgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATa    <  1:157734/100‑1 (MQ=255)                                                                     atgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATa    <  1:4025662/100‑1 (MQ=255)                                                                     atgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATa    <  1:3764844/100‑1 (MQ=255)                                                                       gtgATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAgg  <  1:1124748/100‑2 (MQ=255)                                                                       gtgATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAgg  <  1:2551880/100‑2 (MQ=255)                                                                       gtgATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAgg  <  1:5666495/100‑2 (MQ=255)                                                                            |                                                                                               TTTTTCACTCCAATCTTTGTATTTCGCAAATTGAATTTTGCCGATATAGAAAGGGTTAGCTAAGATGTATGTGATTGAACTAATACTGAAAGGTTTCCCCTTTTTAGTGACATATCCTTTGTGATTCAATGCATTGGCAATTTTACGATAGCCATGTCCTTTGGCATAGC  >  NC_007793/46780‑46949

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