Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,983,835	A→G	D293G (GAC→GGC)	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,835	0	A	G	100.0%	141.4 / NA	42	D293G (GAC→GGC)	aslB	putative anaerobic sulfatase maturation enzyme AslB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/24);  total (18/24)

CTGGATCAAACGCGACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCG  >  NC_000913/3983698‑3983968                                                                                                                                          |                                                                                                                                      cTGGATCAAACGCGACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGgcggcg                                                                                                                                     >  2:765584/1‑140 (MQ=255)     aTCAAACGCGACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGTGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTa                                                                                                                                 <  1:471537/140‑1 (MQ=255)     aTCAAACGCGACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTa                                                                                                                                 <  1:765584/140‑1 (MQ=255)     aTCAAACGCGACGTCGGCAAGACTTTCGTGACGAACATCTAATGGGCGCTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTa                                                                                                                                 <  2:233909/140‑1 (MQ=255)       cAAACGCGACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACg                                                                                                                               <  1:635267/140‑1 (MQ=255)         aaCGCGACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGcc                                                                                                                             >  1:487986/1‑140 (MQ=255)          aCGCGACGTCGGCAACATTTTCGTGATGATTATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCt                                                                                                                            <  1:360174/140‑1 (MQ=255)            gcgACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCtgt                                                                                                                          >  2:675854/1‑140 (MQ=255)                cGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATc                                                                                                                      >  1:574948/1‑140 (MQ=255)                cGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGATGAATGTGGAGCACAAGGGCGGCGTTTACGCCTGTGATc                                                                                                                      >  2:240842/1‑140 (MQ=255)                 gTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCa                                                                                                                     <  2:532507/140‑1 (MQ=255)                    ggCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTa                                                                                                                  >  2:137231/1‑140 (MQ=255)                            tttCGTGATGAATATGGAATGGGCGGTTGACAATTTTATCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATc                                                                                                          <  1:388373/140‑1 (MQ=255)                             ttCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATcc                                                                                                         <  2:464436/140‑1 (MQ=255)                             ttCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATcc                                                                                                         <  1:112020/140‑1 (MQ=255)                               cGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTTTTGAGCACAACGACAGAGTTTACGCCTGTGGTCACTATGTTTATCCGc                                                                                                       >  2:388373/1‑140 (MQ=255)                               cGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGc                                                                                                       >  1:725506/1‑140 (MQ=255)                                  gatgaATATCGAATGGGGGTTTACCAATTTTTTCGGTGCTCCGGGTGCTGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAat                                                                                                    <  2:286682/140‑1 (MQ=255)                                   atgaATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAata                                                                                                   >  2:88975/1‑140 (MQ=255)                                    tgaATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAatat                                                                                                  <  1:88975/140‑1 (MQ=255)                                    tgaATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAatat                                                                                                  <  1:675854/140‑1 (MQ=255)                                    tgaATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAatat                                                                                                  <  2:475381/140‑1 (MQ=255)                                                gggCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAat                                                                                      >  1:715056/1‑140 (MQ=255)                                                   cGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATg                                                                                   <  2:151813/140‑1 (MQ=255)                                                   cGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATg                                                                                   >  2:155254/1‑140 (MQ=255)                                                           aaTTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCa                                                                           <  2:601860/140‑1 (MQ=255)                                                                   cGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTg                                                                   <  2:34972/140‑1 (MQ=255)                                                                           caggTGCGGTTTGCCATCATCATCCAACCTGGGGGTGCTTGGTGATTGCTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATg                                                           <  2:320430/138‑1 (MQ=255)                                                                            gggTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGACTGGGGAATATACATCAGCAAGCAATTGCAGAAATgg                                                          >  1:398932/1‑139 (MQ=255)                                                                                gCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGa                                                      <  1:155254/140‑1 (MQ=255)                                                                                     ttGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTccc                                                 <  2:487986/140‑1 (MQ=255)                                                                                        ccATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGc                                              <  2:574948/140‑1 (MQ=255)                                                                                         catcatCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTCCAGAAAAGATCGATTccccctg                                             >  1:320430/1‑137 (MQ=255)                                                                                             atcaGCCAACCTGTGGGGGCTCTGTGATTGTCGAGCACAGCGGCGGCGGTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAAcag                                         <  2:253762/140‑1 (MQ=255)                                                                                                aGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAAcagcag                                      <  1:80291/140‑1 (MQ=255)                                                                                                       cTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATGGCAGAAATGATGGAGACCCCGTAACAGCAGGCGTTTg                               >  1:253762/1‑140 (MQ=255)                                                                                                       cTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACAAGCAAACAGTTGGAGAAATGATCGACTCCCCGCAACAGCAGGCGTTTg                               >  2:272925/1‑140 (MQ=255)                                                                                                            ggCGCTCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGaa                          >  2:211561/1‑140 (MQ=255)                                                                                                                 tCGGTGATTGTTGAGCACAACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATaa                     >  2:470348/1‑140 (MQ=255)                                                                                                                                 acaACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCg     >  1:433506/1‑140 (MQ=255)                                                                                                                                  caACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAAttt                 <  1:699757/127‑1 (MQ=255)                                                                                                                                  caACGGCGGCGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAAttt                 >  2:699757/1‑127 (MQ=255)                                                                                                                                    aCGGCGGCGTTTACGTCTGTGATCATTATTTTTATCCCCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGgcg  <  1:240842/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                      CTGGATCAAACGCGACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCG  >  NC_000913/3983698‑3983968

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