Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 962,002 G→A 100% intergenic (+13/+22) dmsC → / ← ycaC dimethyl sulfoxide reductase, anaerobic, subunit C/predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606962,0020GA100.0% 54.8 / NA 18intergenic (+13/+22)dmsC/ycaCdimethyl sulfoxide reductase, anaerobic, subunit C/predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/9);  total (9/9)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAATCGGTCGCGGCGTATTCTATGGTTTGCATATGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGGGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGAACAACGTCGCCAGTCCTTCAATATCATTACGCCAGTCGCGATGCAGCTCGCACGC  >  REL606/961894‑962146
                                                                                                            |                                                                                                                                                
gctttATCATTCATTCTGCTTCTGGCAGGGGAATTAATCGGTCGCGGCGTATTCTATGGTTTGCATATGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACg                                                                                                         <  1:194134/146‑1 (MQ=255)
       cattcattCTGCTTCTGGCAGGGGAATTAATCGGTCGCGGCGTATTCTATGGTTTGCATATGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATa                                                                                                  >  1:186394/1‑150 (MQ=255)
          tcattcTGCTTCTGGCAGGGGAATTAATCGGTCGCGGCGTATTCTATGGTTTGCATATGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACt                                                                                               <  2:646609/150‑1 (MQ=255)
                 gCTTCTGGCAGGGGAATTAATCGGTCGCGGCGTATTCTATGGTTTGCATATGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGtcatc                                                                                        >  1:127289/1‑150 (MQ=255)
                               aaTTAATCGGTCGCGGCGTATTCTATGGTTTGCATATGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATc                                                                          >  1:650374/1‑150 (MQ=255)
                                                  aTTCTATGGTTTGCATATGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGaac                                                       >  1:611855/1‑150 (MQ=255)
                                                  aTTCTATGGTTTGCATATGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGaac                                                       >  1:49207/1‑150 (MQ=255)
                                                  aTTCCATGGTTTGCATATGACCGTGGGGATGGCCGTCGCAAGCTCATCATAACAACCGAGGTTCCGGCCCCGGTTCTTTCTTATATCTGGTTCATTCACGTGTTCTAACTGGTCATCAAATTACCATAATCTGGATTCTGGCTAGAGaat                                                       >  1:211794/1‑149 (MQ=255)
                                                              gCATATGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGAACAACGTCGCCAGt                                           >  2:7032/1‑150 (MQ=255)
                                                              gCATATGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGAACAACGTCGCCAGt                                           >  2:160619/1‑150 (MQ=255)
                                                              gCATATGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGGGTCATAACTGGTCATCAAATTACGATAATCAGGACTGTGGTCAGAGAACAACGACGCCaat                                           >  1:123297/1‑148 (MQ=255)
                                                               cttatGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGAACAACGTCGCCAGTc                                          <  1:160171/148‑1 (MQ=255)
                                                                atatGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGAACAACGTCGCCAGTcc                                         <  2:613648/150‑1 (MQ=255)
                                                                     aCCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGAACAACGTCGCCAGTCCTTCaa                                    <  1:10563/150‑1 (MQ=255)
                                                                       cGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGAACAACGTCGCCAGTCCTTCAata                                  <  2:391912/150‑1 (MQ=255)
                                                                        gTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGAACAACGTCGCCAGTCCTTCAatat                                 <  2:294897/150‑1 (MQ=255)
                                                                        gTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGAACAACGTCGCCAGTCCTTCAatat                                 <  2:186930/150‑1 (MQ=255)
                                                                           gggATGGCCGTCGCAAGCTAATCATAACAACCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGAACAACGTCGCCAGTCCTTCAATATCAt                              >  2:391591/1‑150 (MQ=255)
                                                                                                       aacCGAGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGAACAACGTCGCCAGTCCTTCAATATCATTACGCCAGTCGCGATGCAGCTCGcacgc  <  2:549574/150‑1 (MQ=255)
                                                                                                            |                                                                                                                                                
TCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAATCGGTCGCGGCGTATTCTATGGTTTGCATATGACCGTGGGGATGGCCGTCGCAAGCTAATCATAACAACCGGGGTTTCGGCCCCGGTTCTCTTTTATTTCTGCTTCGTTAACGTGTCATAACTGGTCATCAAATTACGATAATCAGGAATGTGGTTAGAGAACAACGTCGCCAGTCCTTCAATATCATTACGCCAGTCGCGATGCAGCTCGCACGC  >  REL606/961894‑962146

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