Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,602,017 (C)7→8 coding (200/2019 nt) hyfB → hydrogenase 4 component B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,602,0101.C91.2% 104.3 / 3.6 34coding (193/2019 nt)hyfBhydrogenase 4 component B
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base C (13/18);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 9.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGC  >  NC_000913/2601883‑2602129
                                                                                                                                |                                                                                                                         
gTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGccccc                                                                                                                    <  2:363349/133‑1 (MQ=255)
               ggACAGACTGGCTATTAATCTTTCCGGCATGACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGc                                                                                                      <  1:160127/133‑1 (MQ=255)
                     aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATTATCAGCGGAATTACGCAATTACATGCTGGTGTAACTATAGTCGCCGGTTTTGCTCCCCCCTTGTGAATTTGCCGAttt                                                                                                >  2:183912/1‑133 (MQ=255)
                     aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGAttt                                                                                                >  1:106194/1‑133 (MQ=255)
                     aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGAttt                                                                                                >  1:469171/1‑133 (MQ=255)
                     aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGAttt                                                                                                >  1:373853/1‑133 (MQ=255)
                     aCTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGGTTTGCCCACCCGTTGTGAATTAGCCGAttt                                                                                                >  1:238961/1‑133 (MQ=255)
                      cTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTa                                                                                               >  2:440798/1‑133 (MQ=255)
                                       cGGCATCACATCGCTGGTGGGCGGCGTGATGGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAg                                                                              <  1:203104/133‑1 (MQ=255)
                                          cATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCct                                                                           <  1:312865/133‑1 (MQ=255)
                                               caTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCgg                                                                      >  1:329306/1‑133 (MQ=255)
                                                aTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGc                                                                     <  1:14774/133‑1 (MQ=255)
                                                       gtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATg                                                              <  1:440798/133‑1 (MQ=255)
                                                          ggcggcGTGATTGGCTTCCTCACCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTg                                                           <  1:158975/133‑1 (MQ=255)
                                                               cgTGATTGGCATCATCAGCGAAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGtt                                                     >  1:93174/1‑133 (MQ=255)
                                                                       gCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCt                                              <  1:62332/133‑1 (MQ=255)
                                                                          tcatcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGc                                           <  2:418271/133‑1 (MQ=255)
                                                                             tcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtgg                                        >  2:483558/1‑133 (MQ=255)
                                                                             tcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGACATTTATAGTGATGGTTATCTACTTGCtgg                                        >  1:160648/1‑133 (MQ=255)
                                                                               accGGAATTTCGCAACTACATGCTTGTGTAACTTTAGTCGCGCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggtg                                      <  2:238961/131‑1 (MQ=255)
                                                                               aGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggtg                                      <  2:64539/133‑1 (MQ=255)
                                                                               aGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggtg                                      <  1:106778/133‑1 (MQ=255)
                                                                                    aaTTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTggtggtggt                                 <  2:106194/133‑1 (MQ=255)
                                                                                       tACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTg                              <  2:270798/133‑1 (MQ=255)
                                                                                        aCGCAATTTCATGCTGGTCTAACTTTAGTCCCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGt                             <  1:183912/133‑1 (MQ=255)
                                                                                         cGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGtt                            >  2:32815/1‑133 (MQ=255)
                                                                                          tcAATTACATGCTGGTATATCTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtct                          <  2:93174/132‑1 (MQ=255)
                                                                                          gCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtc                           <  1:32815/133‑1 (MQ=255)
                                                                                          gCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTGGCCGATTTAACCCTGCGTATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGCGGTGGTTTGTtc                           <  2:373853/133‑1 (MQ=255)
                                                                                          gCAATTACATGCTGGCGTAACTTTAGTCTCCCCTTTTGCCCTCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtc                           <  2:160648/133‑1 (MQ=255)
                                                                                           cAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtct                          <  2:397306/133‑1 (MQ=255)
                                                                                              ttACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtctctc                       >  2:481230/1‑133 (MQ=255)
                                                                                              ttACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTAGTTATCTCCTTGCTGGTGGTGGTTTGTtctctc                       <  1:117532/133‑1 (MQ=255)
                                                                                                     cTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCAt                <  1:481230/133‑1 (MQ=255)
                                                                                                        gtgtAACTTTAGTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGa             >  2:17912/1‑133 (MQ=255)
                                                                                                                ttAGTCGCCCGTTTTGCCCCCCCCTTTTGTGTTTGCCGATTTAACGCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGGTATCTCGTTGCTGGTGGTGGATTGTTCTCTCTATTCATTGACTTATATg     >  2:4115/1‑133 (MQ=255)
                                                                                                                   gTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATgcgc  >  1:428387/1‑133 (MQ=255)
                                                                                                                   gTCGCCCGTTTTGCCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATgcgc  >  2:89325/1‑133 (MQ=255)
                                                                                                                                |                                                                                                                         
GTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGC  >  NC_000913/2601883‑2602129

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: