Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,983,380	C→A	R141R (CGC→CGA)	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,380	0	C	A	100.0%	192.5 / NA	53	R141R (CGC→CGA)	aslB	putative anaerobic sulfatase maturation enzyme AslB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (22/31);  total (22/31)

GTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGAC  >  NC_000913/3983302‑3983455                                                                               |                                                                            gTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGt                                                                           <  1:797120/81‑1 (MQ=255)    cgcgTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGAc                                                                        <  2:797121/81‑1 (MQ=255)     gcgTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGAc                                                                        >  2:797180/1‑80 (MQ=255)     gcgTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACt                                                                       <  1:797123/81‑1 (MQ=255)     gcgTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACt                                                                       <  2:797122/81‑1 (MQ=255)      cgTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTa                                                                      <  2:797124/81‑1 (MQ=255)         ttCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAg                                                                   >  2:797176/1‑81 (MQ=255)         ttCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAg                                                                   >  1:797179/1‑81 (MQ=255)         ttCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAg                                                                   >  1:797177/1‑81 (MQ=255)              gcgGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTg                                                                <  1:797125/79‑1 (MQ=255)                  aGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGAc                                                          >  1:797193/1‑81 (MQ=255)                    catcatTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGAcc                                                         >  1:797195/1‑80 (MQ=255)                     atcatTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCa                                                       >  1:797206/1‑81 (MQ=255)                      tcatTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCAc                                                      <  2:797126/81‑1 (MQ=255)                      tcatTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCAc                                                      >  1:797173/1‑81 (MQ=255)                        atTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGc                                                    <  1:797127/81‑1 (MQ=255)                         ttttCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCa                                                   <  2:797128/81‑1 (MQ=255)                           ttCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCAt                                                  <  1:797129/80‑1 (MQ=255)                                gttgGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCt                                             <  2:797130/80‑1 (MQ=255)                                   ggTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGt                                          <  2:797131/80‑1 (MQ=255)                                    gTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa                                        >  1:797178/1‑81 (MQ=255)                                     tttATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGAt                                       <  2:797132/81‑1 (MQ=255)                                       tATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGc                                     >  1:797183/1‑81 (MQ=255)                                        aTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGc                                     <  1:797134/80‑1 (MQ=255)                                        aTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGc                                     <  2:797133/80‑1 (MQ=255)                                        aTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGc                                     <  2:797135/80‑1 (MQ=255)                                        aTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGc                                     <  2:797136/80‑1 (MQ=255)                                        aTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCg                                    >  1:797166/1‑81 (MQ=255)                                         tCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGt                                   <  2:797137/81‑1 (MQ=255)                                         tCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGt                                   >  1:797169/1‑81 (MQ=255)                                          cGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTg                                  >  1:797187/1‑81 (MQ=255)                                           gCTGGATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGc                                 >  1:797192/1‑81 (MQ=255)                                              ggATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCt                              >  1:797205/1‑81 (MQ=255)                                               gATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGgcc                               <  2:797140/79‑1 (MQ=255)                                               gATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGccc                               <  1:797141/79‑1 (MQ=255)                                               gATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCt                              <  1:797138/80‑1 (MQ=255)                                               gATGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTg                             <  2:797139/81‑1 (MQ=255)                                                 tGGCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTAGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAc                           <  1:797142/81‑1 (MQ=255)                                                   gCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGc                         >  2:797186/1‑81 (MQ=255)                                                   gCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGc                         <  1:797144/81‑1 (MQ=255)                                                   gCCCGCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGc                         <  1:797143/81‑1 (MQ=255)                                                       gCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCt                     >  2:797204/1‑81 (MQ=255)                                                       gCCTGAGATCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCt                     <  2:797145/81‑1 (MQ=255)                                                            agaTCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCaa                 <  2:797146/80‑1 (MQ=255)                                                              aTCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCaaaaa              >  2:797214/1‑81 (MQ=255)                                                              aTCCACAATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCCGGTGATGCGTGCCCTGACGCTCCTGCaaaaa              >  2:797190/1‑81 (MQ=255)                                                                  acaATCAATATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAAcat           >  1:797189/1‑80 (MQ=255)                                                                    aatcaatATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAAcatcat        <  1:797147/81‑1 (MQ=255)                                                                       caatATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTc     >  2:797198/1‑81 (MQ=255)                                                                        aatATCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTc     <  1:797148/80‑1 (MQ=255)                                                                          tatCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGAc  <  2:797149/81‑1 (MQ=255)                                                                          tatCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGAc  <  1:797150/81‑1 (MQ=255)                                                                           atCGAGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGAc  >  1:797229/1‑80 (MQ=255)                                                                               |                                                                            GTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGAC  >  NC_000913/3983302‑3983455

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑