Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,607,273 C→T F119F (TTC→TTT hyfF → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,607,2730CT100.0% 103.3 / NA 30F119F (TTC→TTThyfFhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/15);  total (15/15)

GTCTGGGCGGTTTGTTCCTCGCCATTCTTGGTGTGATTGGCTTTCTCACCGGTATTTACTCGATTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGC  >  NC_000913/2607137‑2607407
                                                                                                                                        |                                                                                                                                      
gTCTGGGCGGTTTGTTCCTCGCCATTCTTGGTGTGATTGGCTTTCTCACCGGTATTTACTCGATTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCa                                                                                                                                      <  1:89062/139‑1 (MQ=255)
  cTGGGCGGTTTGTTCCTCGCCATTCTTGGTGTGATTGGCTTTCTCACCGGTATTTACTCGATTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATc                                                                                                                                    >  1:523916/1‑139 (MQ=255)
    gggCGGTTTGTTCCTCGCCATTCTTGGTGTGATTGGCTTTCTCACCGGTATTTACTCGATTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTg                                                                                                                                  >  1:274784/1‑139 (MQ=255)
     ggCGGTTTGTTCCTCGCCATTCTTGGTGTGATTGGCTTTCTCACCGGTATTTACTCGATTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGt                                                                                                                                 <  2:523916/139‑1 (MQ=255)
     ggCGGTTTGTTCCTCGCCATTCTTGGTGTGATTGGCTTTCTCACCGGTATTTACTCGATTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCAACTGt                                                                                                                                 >  2:382418/1‑139 (MQ=255)
                    gCCATTCTTGGTGTGATTGGCTTTCTCACCGGTATTTACTCGATTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCAcc                                                                                                                      >  2:528637/1‑135 (MQ=255)
                    gCCATTCTTGGTGTGATTGGCTTTCTCACCGGTATTTACTCGATTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCAcc                                                                                                                      <  1:528637/135‑1 (MQ=255)
                          cTTGGTGTGATTGGCTTTCTCACCGGTATTTACTCGATTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTgg                                                                                                            <  2:274784/139‑1 (MQ=255)
                             ggTGTGATTGGCTTTCTCACCGGTATTTACTCGATTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGgttg                                                                                                         <  2:138078/139‑1 (MQ=255)
                                       gCTTTCTCACCGGTATTTACTCGATTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCaa                                                                                               <  2:446756/139‑1 (MQ=255)
                                                              aTTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGggcgg                                                                        >  1:78518/1‑139 (MQ=255)
                                                              aTTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGAGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGggcag                                                                        >  1:180009/1‑137 (MQ=255)
                                                                  gCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTAc                                                                                                     >  2:220293/1‑106 (MQ=255)
                                                                  gCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTAc                                                                                                     <  1:220293/106‑1 (MQ=255)
                                                                  gCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCa                                                                                                   >  1:270990/1‑108 (MQ=255)
                                                                  gCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCa                                                                                                   <  2:270990/108‑1 (MQ=255)
                                                                                gAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAAcc                                                                                          <  2:22439/103‑1 (MQ=255)
                                                                                gAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAAcc                                                                                          >  1:22439/1‑103 (MQ=255)
                                                                                      gCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCt                                                >  2:505104/1‑139 (MQ=255)
                                                                                       cacaCGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTc                                               <  2:78518/139‑1 (MQ=255)
                                                                                            ggCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATaa                                                                                            <  1:198257/89‑1 (MQ=255)
                                                                                            ggCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATaa                                                                                            >  2:198257/1‑89 (MQ=255)
                                                                                                      cACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCa                                                                                                                        <  1:525755/51‑1 (MQ=255)
                                                                                                      cACCCGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCa                                                                                                                        >  2:525755/1‑51 (MQ=255)
                                                                                                          cGTTACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCAtt                            <  1:505104/139‑1 (MQ=255)
                                                                                                            ttACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCAtgctgc                                                                                                               >  2:164310/1‑54 (MQ=255)
                                                                                                            ttACGCTGTGCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCAtgctgc                                                                                                               <  1:164310/54‑1 (MQ=255)
                                                                                                                   gtgCGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAg                   >  2:93788/1‑139 (MQ=255)
                                                                                                                      cGATTACTACGGTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGt                >  1:393374/1‑139 (MQ=255)
                                                                                                                                 gTTTCTTTCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGa     <  2:393374/139‑1 (MQ=255)
                                                                                                                                   ttctttCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAg   >  1:219299/1‑139 (MQ=255)
                                                                                                                                    tctttCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGc  <  1:393998/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GTCTGGGCGGTTTGTTCCTCGCCATTCTTGGTGTGATTGGCTTTCTCACCGGTATTTACTCGATTGGCTACATGCGTCATGAAGTGGCACACGGCGAGCTTTCACCCGTTACGCTGTGCGATTACTACGGTTTCTTCCATCTGTTTTTGTTCACCATGCTGCTGGTTGTTACCAGCAATAACCTGATTGTGATGTGGGCGGCGATCGAAGCCACCACCTTAAGCTCGGCGTTTCTGGTAGGCATTTACGGTCAGCGTTCATCGCTGGAAGC  >  NC_000913/2607137‑2607407

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