Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,983,666	T→C	F237L (TTC→CTC)	aslB →	putative anaerobic sulfatase maturation enzyme AslB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,983,666	0	T	C	100.0%	148.1 / NA	40	F237L (TTC→CTC)	aslB	putative anaerobic sulfatase maturation enzyme AslB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/19);  total (21/19)

TGATATTCAGGGTGAGCTAACGGAATGGTCGGTGCGCCCCGAGGAGTTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCAAACGCG  >  NC_000913/3983620‑3983711                                               |                                              tGATATTCAGGGTGAGCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTc                                             <  1:564229/49‑1 (MQ=255)   atatTCAGGGTGAGCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTCgg                                           >  1:1917753/1‑49 (MQ=255)    tatTCAGGGTGAGCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTCGGt                                          <  1:78568/49‑1 (MQ=255)     atTCAGGGTGAGCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTCGGt                                          >  1:1102191/1‑48 (MQ=255)      ttCAGGGTGAGCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTCGGt                                          >  1:1037755/1‑47 (MQ=255)        cAGGGTGAGCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAg                                       <  1:121263/48‑1 (MQ=255)         aGGGTGAGCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGtt                                     >  1:1030133/1‑49 (MQ=255)         aGGGTGAGCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGt                                      <  1:1539950/48‑1 (MQ=255)             tGAGCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCtg                                 >  1:2289855/1‑49 (MQ=255)              gAGCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCtgg                                >  1:1117986/1‑49 (MQ=255)               aGCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCtggt                               >  1:1056317/1‑49 (MQ=255)               aGCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCtggt                               <  1:1203452/49‑1 (MQ=255)                gCTAACGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCtggtg                              <  1:1995900/49‑1 (MQ=255)                  tAACGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCTGGTGGc                            >  1:1479022/1‑49 (MQ=255)                   aaCGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCTGGTGGc                            >  1:740683/1‑48 (MQ=255)                    aCGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCTGGTGGCGa                          >  1:2257157/1‑49 (MQ=255)                     cGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCTGGTGGCGat                         <  1:779735/49‑1 (MQ=255)                     cGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCTGGTGGCGa                          >  1:231985/1‑48 (MQ=255)                     cGGAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCTGGTGGCGa                          >  1:1098267/1‑48 (MQ=255)                      ggAATGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCTGGTGGCGata                        <  1:1979770/49‑1 (MQ=255)                         aTGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTc                     <  1:155932/49‑1 (MQ=255)                          tGGTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTc                     >  1:1178803/1‑48 (MQ=255)                            gTCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTCGAc                  >  1:236568/1‑49 (MQ=255)                             tCGGTGCGCCCCGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTCGAcc                 <  1:1936757/49‑1 (MQ=255)                               ggTGCGCCCCGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTCGACCAc               <  1:435451/49‑1 (MQ=255)                                    gcCCCGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGAt          >  1:90370/1‑49 (MQ=255)                                     ccccGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATc         <  1:1597609/49‑1 (MQ=255)                                      cccGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATc         >  1:2007844/1‑48 (MQ=255)                                      cccGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCa        <  1:1808849/49‑1 (MQ=255)                                       ccGAGGAGCTCTGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCa        >  1:33133/1‑48 (MQ=39)                                       ccGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCa        >  1:1198679/1‑48 (MQ=255)                                        cGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCaaa      <  1:1915865/49‑1 (MQ=255)                                        cGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCaaa      <  1:607883/49‑1 (MQ=255)                                        cGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCaaa      <  1:1255751/49‑1 (MQ=255)                                        cGAGGAGCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCaaa      <  1:886170/49‑1 (MQ=255)                                         gaggagCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCAAAc     >  1:2049680/1‑49 (MQ=255)                                          aggagCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCAAAcg    <  1:336435/49‑1 (MQ=255)                                          aggagCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCAAAcg    >  1:1564877/1‑49 (MQ=255)                                            gagCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCAAAcgcg  <  1:457580/49‑1 (MQ=255)                                            gagCTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCAAAcgc   >  1:1006824/1‑48 (MQ=255)                                               |                                              TGATATTCAGGGTGAGCTAACGGAATGGTCGGTGCGCCCCGAGGAGTTCGGTGAGTTTCTGGTGGCGATATTCGACCACTGGATCAAACGCG  >  NC_000913/3983620‑3983711

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