Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	46,867	G→A	F218F (TTC→TTT)	SAUSA300_RS00200 ←	recombinase RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	46,867	0	G	A	98.4%	180.8 / ‑3.8	62	coding (654/1629 nt)	SAUSA300_RS00200	recombinase RecB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (4/57);  minor base . (0/1);  total (4/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

TTTTTCACTCCAATCTTTGTATTTCGCAAATTGAATTTTGCCGATATAGAAAGGGTTAGCTAAGATGTATGTGATTGAACTAATACTGAAAGGTTTCCCCTTTTTAGTGACATATCCTTTGTGATTCAATGCATTGGCAATTTTACGATAGCCATGTCCTTTGGCATAGCACTCGAAT  >  NC_007793/46780‑46957                                                                                        |                                                                                           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:5045592/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt             <  1:2753798/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt             <  1:538591/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt             <  1:3583442/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt             <  1:538974/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt             <  1:1776020/100‑1 (MQ=255)                                                                    tatgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAt             <  1:5660355/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:3306724/100‑2 (MQ=255)                                                                     atgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCAtt            <  1:2955066/100‑2 (MQ=255)                                                                     atgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATa            <  1:3764844/100‑1 (MQ=255)                                                                     atgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATa            <  1:4025662/100‑1 (MQ=255)                                                                     atgtGATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATa            <  1:157734/100‑1 (MQ=255)                                                                       gtgATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAgg          <  1:1124748/100‑2 (MQ=255)                                                                       gtgATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAgg          <  1:5666495/100‑2 (MQ=255)                                                                       gtgATAGAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAgg          <  1:2551880/100‑2 (MQ=255)                                                                            agAACTAATACTAAAAGGTTTCCCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAGGACTCg     <  1:369472/99‑1 (MQ=255)                                                                            agAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAGGACTCg     <  1:936344/99‑1 (MQ=255)                                                                            agAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAGGACTCg     <  1:1582789/99‑1 (MQ=255)                                                                            agAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAGGACTCg     <  1:2520688/99‑1 (MQ=255)                                                                            agAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAGGACTCg     <  1:5372755/99‑1 (MQ=255)                                                                            agAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAGGACTCg     <  1:5279726/99‑1 (MQ=255)                                                                            agAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAGGACTCg     <  1:5207442/99‑1 (MQ=255)                                                                            agAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAGGACTCg     <  1:265076/99‑1 (MQ=255)                                                                            agAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTAGTGATTCAATGCATTGGCTAATTTACGATAGCCATGTCCTTTGGCATAGGACTCg     <  1:3720532/99‑1 (MQ=255)                                                                             gAACTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAGGACTcat    <  1:480023/100‑3 (MQ=255)                                                                               aCTAATACTAAAAGGTTTACCCTTTTTAGTGACATAGCCTTTGTGATTCAATGCATTGGCTATTTTACGATAGCCATGTCCTTTGGCATAGGACTCGagt  <  1:4507590/100‑3 (MQ=255)                                                                                        |                                                                                           TTTTTCACTCCAATCTTTGTATTTCGCAAATTGAATTTTGCCGATATAGAAAGGGTTAGCTAAGATGTATGTGATTGAACTAATACTGAAAGGTTTCCCCTTTTTAGTGACATATCCTTTGTGATTCAATGCATTGGCAATTTTACGATAGCCATGTCCTTTGGCATAGCACTCGAAT  >  NC_007793/46780‑46957

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