Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,140,294 T→C S9G (AGC→GGC)  hybF ← protein involved with the maturation of hydrogenases 1 and 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,140,2940TC96.3% 90.4 / ‑4.9 27S9G (AGC→GGC) hybFprotein involved with the maturation of hydrogenases 1 and 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (13/13);  minor base A (0/1);  total (13/14)
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.97e-01

AATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGT  >  NC_000913/3140169‑3140421
                                                                                                                             |                                                                                                                               
aaTTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAagagtc                                                                                                                    <  1:1001137/139‑3 (MQ=255)
aaTTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGAc                                                                                                                    >  2:889004/1‑139 (MQ=255)
           tAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCAt                                                                                                         <  1:889004/139‑1 (MQ=255)
                  gacggCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGcc                                                                                                  <  1:196932/139‑1 (MQ=255)
                              ctcAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGAcc                                                                                      <  2:340685/139‑1 (MQ=255)
                                cAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCaa                                                                                    <  2:295953/139‑1 (MQ=255)
                                        gagaGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAgg                                                                            >  1:17601/1‑139 (MQ=255)
                                        gagaGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAgg                                                                            >  2:858216/1‑139 (MQ=255)
                                        gagaGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAgg                                                                            >  2:15692/1‑139 (MQ=255)
                                             cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGAc                                                                                                                    <  2:581614/94‑1 (MQ=255)
                                             cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGAc                                                                                                                    >  1:581614/1‑94 (MQ=255)
                                             cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGcgcg                                                                       >  2:941395/1‑139 (MQ=255)
                                                  cAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGcc                                                                  <  2:17601/139‑1 (MQ=255)
                                                  cAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGcc                                                                  <  2:854155/139‑1 (MQ=255)
                                                          aGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCgtgt                                                          >  1:566019/1‑139 (MQ=255)
                                                                acGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGc                                                                                                                            >  1:792932/1‑67 (MQ=255)
                                                                acGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGc                                                                                                                            <  2:792932/67‑1 (MQ=255)
                                                                      gTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGcc                                                                  >  2:658270/1‑119 (MQ=255)
                                                                      gTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCACTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGcc                                                                  <  1:658270/119‑1 (MQ=255)
                                                                              cTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGAc                                      >  2:240335/1‑139 (MQ=255)
                                                                                 aaCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTgg                                   <  1:15692/139‑1 (MQ=255)
                                                                                                cGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCAt                    <  1:579003/139‑1 (MQ=255)
                                                                                                  ccTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACATCAGGGCGCGCCGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTc                  <  1:746281/139‑1 (MQ=255)
                                                                                                          tGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGa                                              >  2:341529/1‑103 (MQ=255)
                                                                                                          tGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGa                                              <  1:341529/103‑1 (MQ=255)
                                                                                                                 tttCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCg   >  1:54859/1‑139 (MQ=255)
                                                                                                                 tttCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGt  >  1:530302/1‑139 (MQ=255)
                                                                                                                             |                                                                                                                               
AATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGT  >  NC_000913/3140169‑3140421

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