Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,603,828 G→A A671T (GCC→ACC)  hyfB → hydrogenase 4 component B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,603,8280GA95.0% 49.0 / ‑7.1 20A671T (GCC→ACC) hyfBhydrogenase 4 component B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (8/11);  minor base T (0/1);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATACGCGCGCGTATGCACTCCCGCCGCGGGCCGGGGATCTGGCAGGATTATCGCGATATCCACAAACTGTTTAAACGCCAGGAAGTTGCGCCGACATCTTCAGGTCTGATGTTCCGCCTGAT  >  NC_000913/2603629‑2604070
                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
gCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGNc                                                                                                                                                                                                 <  1:395340/251‑1 (MQ=255)
  ngAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACCATg                                                                                                                                                                                                                               >  1:300558/2‑219 (MQ=255)
                                                                gATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTATCACAGGCCGTGACTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTccn                                                                                                                                 <  1:104887/251‑2 (MQ=255)
                                                                       nnnnnntCNNNNNCATNGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACCATg                                                                                                                                                                                                                               <  2:300558/144‑1 (MQ=255)
                                                                                                 aCCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACGGGTTGTTCGGCTAATCGCTATTTCCGTCTAAg                                                                                                                                                                                                                                           <  1:406951/112‑1 (MQ=255)
                                                                                                         ncTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACAATGAGACAAACTCttt                                                                                                                                                                                                                  >  1:299268/2‑129 (MQ=255)
                                                                                                             gCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGc                                                                                                                              <  1:140282/209‑1 (MQ=255)
                                                                                                                        aCAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCCCTTTTTACGGGTATTTCCCGGCATATACGCGCGCGTATGCACTCCCGCCGCGGGCCGGGGATCTGGCAGGATTa                                                                          <  1:205071/250‑1 (MQ=255)
                                                                                                                           gCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTGTGTGGTCGCCGCACTGGTTGTGATGCTAATCACTATTACCGTCTAAGGAAATCACCATg                                                                                                                                                                                                                               >  1:472555/1‑98 (MQ=255)
                                                                                                                                         gcgACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACCAGGAGACAAACTCTTTGCGACGGATAGCTGGTCATGGTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAAccc                                                                                                                                                      <  1:185708/157‑1 (MQ=255)
                                                                                                                                           gACTTTCGTCTCTATTGCCTGGATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGt                                                                                                                                                                                                  <  1:240135/111‑1 (MQ=255)
                                                                                                                                                  gTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACCATGAGACAAACTCTTTNNGACgg                                                                                                                                                                                                           >  2:13100/1‑95 (MQ=255)
                                                                                                                                                  gTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAAcc                                                                                                                                                       <  1:13100/147‑1 (MQ=255)
                                                                                                                                                             cTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATAcgcgcg                                                                                                                   <  1:163395/172‑1 (MQ=255)
                                                                                                                                                                tATGTGGTAGACGCACTGGTTGTGCTNCTAATCNNTCNTACCGTCTAGGAAATAACCNTGNNNCNNNCNCTTTGNGANNGATATCTGGNCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTANGGGTATTTCNNGGCAGATACGCGCGCGTATGCNCTCNNGNNGNGGGNCg                                                                                           >  2:106874/1‑192 (MQ=255)
                                                                                                                                                                      gTCGCCGCACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCGCCATGAGACAAACTCTTTGCGACGGATATCTGGGCATGTGTGCGTTAGCACAGGCCGTGATTctgct                                                                                                                                                                 <  1:206113/117‑1 (MQ=255)
                                                                                                                                                                             cACTGGTTGTGCTGCTAATCGCTATTACCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTNNGcac                                                                                                                                                                                 >  2:7559/1‑94 (MQ=255)
                                                                                                                                                                                          gctAATCGCTATTACCGTCTAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCGGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATACGCGCGCGTATGCACTCCCGCCGCGGGCCGGGGATCTGGCAGGATTATCGCGATATCCACAAACTGTTTAAACGCCAGGAAGTTGCGCCGACATCTTCAGGTCTGATGTTCCg        <  1:106874/249‑1 (MQ=255)
                                                                                                                                                                                               tCGCTATTACCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGAAGGATATCTGGTCATTTTTACGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATACGCGCGCGTATGCACTCCCGCCGCGGGCAGGGGATCTGGAAGGATTATCGCGATATCCACAAACTGTTTAAACGCCAGGCAGTTGCGCCGACATCTTCAGGTCTGATGTTCCGCCTGAt  >  1:3047/1‑251 (MQ=255)
                                                                                                                                                                                                  cTATTACCGTCTAAGGAAATCACCATGAGACAAAATCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATACGCGCGCGTATGCACTCCCGCCGCGGGCCGGGGATCTGGCAGGATTATCGCGATATCCACAAACTGTTTAAACGCCAGGAAGTTGCGCCGACATCTTCAGGTCTGa               >  1:174041/1‑235 (MQ=255)
                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
GCTGAGGCTAAATAAAGGTCTTGCGCACGTCACCGCCAGGGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATACGCGCGCGTATGCACTCCCGCCGCGGGCCGGGGATCTGGCAGGATTATCGCGATATCCACAAACTGTTTAAACGCCAGGAAGTTGCGCCGACATCTTCAGGTCTGATGTTCCGCCTGAT  >  NC_000913/2603629‑2604070

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