Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,983,365 C→T I136I (ATC→ATT aslB → putative anaerobic sulfatase maturation enzyme AslB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,983,3650CT100.0% 66.2 / NA 24I136I (ATC→ATTaslBputative anaerobic sulfatase maturation enzyme AslB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/14);  total (10/14)

GCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATGAAACAACTGCCCGCGATGGAC  >  NC_000913/3983178‑3983600
                                                                                                                                                                                           |                                                                                                                                                                                                                                           
gngAACCAACGCAACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTg                                                                                                                                                                              >  1:393042/1‑251 (MQ=255)
                                 gccgtgccgtAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATTAACAATCAATAGCGCGGGACTAAAGGGGGCAGACCCGCGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGAc                                                                                                                                                   <  1:450733/245‑1 (MQ=255)
                                                                         gAAGAGAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGAGGCGCATCATTCTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGCTTCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGt                                                                                                     >  1:485132/1‑251 (MQ=255)
                                                                                 gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGt                                                                                              >  1:377408/1‑250 (MQ=255)
                                                                                           ccAGACTAACGGCGTGCTGCTTGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCGTGAGATTCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGAGGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGAttttt                                                                                   <  1:439217/251‑1 (MQ=255)
                                                                                                             gctgGATGACGAATGGTGCTCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGTggn                                                                                                                                                                                                               <  1:48170/109‑2 (MQ=255)
                                                                                                             gctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGCGATTCCCAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTaa                                                                                                                            >  1:109232/1‑192 (MQ=255)
                                                                                                             gcngGANGNNGNNTGGNGCGCGGTTCTCGCGGAGCATCNNTTTCTTGTTNGTTTATCGCTGGATGGCCCGCCTGAGATTCACAATCAATATCGCGTGACTAAANGTGGCATACNNGCNCATAANCNNNTNNNGCNTGCCCTGACGCTCCTGCAANAANNTCATGTNGACTATAACGTGCTGGTCTGCGTTaa                                                                                                                            <  2:109232/192‑1 (MQ=255)
                                                                                                                ggATGACGTATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGTgg                                                                                                                                                                                                                <  1:40267/105‑1 (MQ=255)
                                                                                                                ggATGACGAATGGTGAGCGTTTCTCGCGGGGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCATGAGATTCACAATCAATATNNCGTGACNna                                                                                                                                                                                                                      >  2:40267/1‑99 (MQ=255)
                                                                                                                      cGAAGGGTGCGCGTGTCTCGCGGAGCATCATTTTCTTGGTGGTTTATCGCTGGATGGCGCGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGTggn                                                                                                                                                                                                               <  1:51204/100‑2 (MQ=255)
                                                                                                                                     tctcGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATTCACAATCAATATcgcn                                                                                                                                                                                                                             <  1:52823/71‑2 (MQ=255)
                                                                                                                                        cgcgGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGTGGCa                                                                                                                                                                                                              <  2:41001/83‑1 (MQ=255)
                                                                                                                                        cgcgGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGTGGCa                                                                                                                                                                                                              >  1:41001/1‑83 (MQ=255)
                                                                                                                                         gcgGANCANCNTTNNCNTGTTNGTTTATCGCTGGATGGCCCGCCTGAGATTCAcaatc                                                                                                                                                                                                                                      >  2:226937/1‑58 (MQ=255)
                                                                                                                                               catcatTTTCTTGTTGGTGTATCGCTGGATGGCCCGCCGGAGATTCACAATCAATATCGCGTGACTAAAGGTGGCATACCCACGCATAAGCTGGTGATGCTTGCCCTGACGATCCTGCAAAAACATCATGTCGACTATAACg                                                                                                                                            <  1:490932/142‑1 (MQ=255)
                                                                                                                                                           gttgGTTTATCGCTGGATGGCACGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCg                                                                                                                                                     >  1:136265/1‑121 (MQ=255)
                                                                                                                                                                   aTCGCTGGATGGCCCGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGta                                                                                             <  1:197094/169‑1 (MQ=255)
                                                                                                                                                                   aTCGCNGGATCNCNNGCCTGAGATTCACAATCAANATCGCGTGACTAAAGNTGGCNGACCCNCgn                                                                                                                                                                                                     >  2:197094/1‑64 (MQ=255)
                                                                                                                                                                         ggATGGNCCGCCTGAGANTCANNATCANTNTNGCGTGACTAAAGGTGGCAGACCCACGCATNAGCTGGTGATGCGTGCCCTGACGNNCCNNNAANNAcatc                                                                                                                                                           <  2:367746/101‑1 (MQ=255)
                                                                                                                                                                         ggATGGCCCGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAAcatc                                                                                                                                                           >  1:367746/1‑101 (MQ=255)
                                                                                                                                                                             ggCCCGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCTGGAGTCGAATTCATCCAGTTTATTCCGGTGGTNg                                      <  1:395704/214‑1 (MQ=255)
                                                                                                                                                                             ggCCCGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGGTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATGAAACAACTGCCCGCGATGGAc  <  1:391922/250‑1 (MQ=255)
                                                                                                                                                                              gCCCGCTTTAGATTCACAATCAATATCGCGTGCCTAAAGGTGGCAGTCCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTcc                                                                                                                                                                       <  1:186290/84‑1 (MQ=255)
                                                                                                                                                                                           |                                                                                                                                                                                                                                           
GCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGACGCTCCTGCAAAAACATCATGTCGACTATAACGTGCTGGTCTGCGTTAATCGCACCAGCGCGCAGCAACCGTTGCAGGTATATGATTTTTTGTGCGATGCGGGAGTCGAATTCATCCAGTTTATTCCGGTGGTCGAGCGCCTGGCTGATGAAACAACTGCCCGCGATGGAC  >  NC_000913/3983178‑3983600

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