Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,586,131	C→T	E149K (GAA→AAA)	SAUSA300_RS13290 ←	peptidase M28

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,586,131	0	C	T	100.0%	69.1 / NA	25	E149K (GAA→AAA)	SAUSA300_RS13290	peptidase M28
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (5/20);  total (5/20)

TTCGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTCAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTAAGACATGTGGCGG  >  NC_007793/2586046‑2586224                                                                                      |                                                                                              ttCGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTccgggt                                                                                 <  1:3800459/100‑5 (MQ=255)   cGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGcta                                                                               >  1:4056755/1‑100 (MQ=255)   cGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGcac                                                                               <  1:4580247/100‑3 (MQ=255)   cGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTAGcac                                                                               <  1:949209/100‑3 (MQ=255)    gCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGctac                                                                              <  1:1003647/100‑1 (MQ=255)     cAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCtcggcatg                                                                             <  1:5313596/100‑8 (MQ=255)     cAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGctact                                                                             <  1:2785794/100‑1 (MQ=255)     cAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTAGcacac                                                                             <  1:4288919/100‑5 (MQ=255)      aGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGctactg                                                                            <  1:2517538/100‑2 (MQ=255)      aGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGctacta                                                                            <  1:540830/100‑1 (MQ=255)      aGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGctacta                                                                            <  1:2039349/100‑1 (MQ=255)        ttAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAac                                                                          <  1:2126058/100‑1 (MQ=255)         tAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAaca                                                                         <  1:1427087/100‑1 (MQ=255)                  aTCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCt                                                                <  1:2539278/100‑1 (MQ=255)                            aGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACatat                                                      <  1:3284327/100‑1 (MQ=255)                                     aaTTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGAttt                                             >  1:3554679/1‑100 (MQ=255)                                              tatGGCGTAACCATATTGCCTTCGTCAACGCAAGCTTCTTTAGCGTCTTCCTTGCTCCTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCt                                    <  1:139208/100‑1 (MQ=255)                                                 ggCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCcttt                                 <  1:5524911/100‑1 (MQ=255)                                                        ccATATTGCTTACTTCAACGCCAGCTACTTTAGCTTCTTCCTTGCTTCTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGtta                          <  1:778552/100‑2 (MQ=255)                                                              tGCCTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAgg                    >  1:5269052/1‑100 (MQ=255)                                                                 cTACTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCgg                 <  1:920210/100‑2 (MQ=255)                                                                   tctccggcGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTa               <  1:2397789/93‑1 (MQ=255)                                                                   aCTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTa               >  1:1348074/1‑100 (MQ=255)                                                                   aCTTCAACGCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTa               <  1:51279/100‑1 (MQ=255)                                                                     ttgggggCCAGCTTCTTTAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTAAg             >  1:1951181/7‑100 (MQ=255)                                                                               gcttctttagcttcttcCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTAAGACATGTGGca   <  1:3568215/100‑2 (MQ=255)                                                                                cttctttagcttcttcCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTAAGACATGTGGCgg  <  1:4152649/100‑1 (MQ=255)                                                                                      |                                                                                              TTCGCAGTTAAATATTTATCATTTGCAAGCACTTCAAATTCACTATATGGCGTAACCATATTGCCTACTTCAACGCCAGCTTCTTCAGCTTCTTCCTTGCTACTAACACCAATATCTATAAACATATTTTTGATTTCCATTGGCTTTTTACGTTCTTCAGGCGTTAAGACATGTGGCGG  >  NC_007793/2586046‑2586224

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