Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,602,118 A→G T101A (ACT→GCT)  hyfB → hydrogenase 4 component B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,602,1180AG93.8% 39.1 / ‑4.6 17T101A (ACT→GCT) hyfBhydrogenase 4 component B
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base G (8/7);  minor base T (0/1);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

GCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAACGTCGAT  >  NC_000913/2601880‑2602299
                                                                                                                                                                                                                                              |                                                                                                                                                                                     
gCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCCCATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTCGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGGCTTATATGcgcg                                                                                                                                                                           >  1:354624/1‑251 (MQ=255)
gCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGCGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGCGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGGCTTCTATGcgcg                                                                                                                                                                           >  1:354601/1‑251 (MQ=255)
             ngTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGCCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTtata                                                                                                                                                                                 >  1:71067/2‑232 (MQ=255)
                                   aGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCGCTCTATTCATTGGCTTATATGcgcg                                                                                                                                                                           <  1:105564/216‑1 (MQ=255)
                                            caatcacgccgcccaccagcgatgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGAGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTCACGTTTTACCCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGGGGTGGTTTGTTCTCTCTATTCATTGGTTTATATGAGCAAATACGAGGGCAAGGGCGCGGCGGCGATGGCCTTTTTTATGAat                                                                                                                               >  1:158269/24‑251 (MQ=255)
                                                       ctgttggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGGCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGAGAGCCTCTCGGCATTTATGGTTCTGTTTATCTCCTTGCTGGTGGTTGTTTGTTCTCTCTATTCATTGTCTTATATGcgcg                                                                                                                                                                           <  1:252916/192‑1 (MQ=255)
                                                                   ntGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGGCTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGc                                                                                                        >  1:212423/2‑251 (MQ=255)
                                                                          gCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTCTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGGCTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATTGTTGCCctgctg                                                                                                 <  1:151479/251‑1 (MQ=255)
                                                                                                          ggTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGGCTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTAATCGTGc                                                                  <  1:25871/250‑1 (MQ=255)
                                                                                                          ggTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGATGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGGCTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCt                                                                 >  1:317626/1‑251 (MQ=255)
                                                                                                                         gCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAAACCTGCGAATGGATCGACACTTGGACTTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGGCTTATATGCGCGCATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTc                                                  >  1:306776/1‑251 (MQ=255)
                                                                                                                                   cccccccTTTTGAAGTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGGCTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGt                                                                           <  1:25587/216‑1 (MQ=255)
                                                                                                                                                     ccGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTTGCTTATATGCGCTAATACGAGGGCAAAGGCGCGGCGGCTATGGGCTTCGTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCttttn                                                                              <  1:36545/195‑2 (MQ=255)
                                                                                                                                                     ccGATGTAACCCTGCGAATGGATAGCCTCGCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGGCTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTg                                                            <  1:243882/213‑1 (MQ=255)
                                                                                                                                                                         gATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGATGGTGGTGGTTTATTCTCTCTATTCATTGGCTTATATGAGCGCATAAGAGGGAAAAGGCGCGGCGGCGATGGGCTTCTTTATGACTATTTTCTTAGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAAAGTCGAt  >  1:57392/1‑251 (MQ=255)
                                                                                                                                                                                                gCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTTATTGGCTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTg                                      <  1:151971/192‑1 (MQ=255)
                                                                                                                                                                                                gCTGGNTATCTNCNNGCTGGTGGTGGTTTGTTCTNTCTATTTATTGGCTTNTATGNGCGAANAcn                                                                                                                                                                     >  2:151971/1‑64 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                     
GCTGTTTTTACTCGGTCTGGACAGACTGGCTATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGCCCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTCTCTCTATTCATTGACTTATATGCGCGAATACGAGGGCAAAGGCGCGGCGGCGATGGGCTTCTTTATGAATATTTTCATCGCATCGATGGTTGCCCTGCTGGTGATGGACAACGCTTTTTGGTTCATCGTGCTGTTTGAAATGATGTCGCTGTCTTCCTGGTTTCTGGTCATTGCCAGGCAGGATAAAACGTCGAT  >  NC_000913/2601880‑2602299

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑