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,2940TC95.7% 75.4 / ‑6.2 23S9G (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 (10/12);  minor base A (0/1);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCA  >  NC_000913/3140162‑3140423
                                                                                                                                    |                                                                                                                                 
tgCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGcg                                                                                                                             <  2:1026602/138‑2 (MQ=255)
ggCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCa                                                                                                                             <  1:820080/139‑1 (MQ=255)
               aaCTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATg                                                                                                              >  2:902843/1‑139 (MQ=255)
               aaCTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATg                                                                                                              >  1:964884/1‑139 (MQ=255)
                 cTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCa                                                                                                            <  1:230925/139‑1 (MQ=255)
                         gacggCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGcc                                                                                                    <  1:432655/139‑1 (MQ=255)
                          acggCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCAt                                                                                                               >  2:980750/1‑127 (MQ=255)
                          acggCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCACTGGCAAAGAGACAACTCAt                                                                                                               <  1:980750/127‑1 (MQ=255)
                                 tctcCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGc                                                                                            <  2:964884/139‑1 (MQ=255)
                                               gagaGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGccct                                                                                                                                 >  1:1028898/1‑88 (MQ=255)
                                               gagaGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCa                                                                                                                             <  1:826149/92‑1 (MQ=255)
                                               gagaGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCa                                                                                                                             >  2:826149/1‑92 (MQ=255)
                                                    cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGcgcg                                                                         >  1:386561/1‑139 (MQ=255)
                                                                                     cTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGc                                                                                                                              >  2:1066229/1‑53 (MQ=255)
                                                                                     cTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGc                                                                                                                              <  1:1066229/53‑1 (MQ=255)
                                                                                                gctgctCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGaaa                             >  2:182029/1‑139 (MQ=255)
                                                                                                    ctCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCaa                                                                                                                            <  1:97713/40‑1 (MQ=255)
                                                                                                    ctCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCaa                                                                                                                            >  2:97713/1‑40 (MQ=39)
                                                                                                    ctCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGAt                         <  1:182029/139‑1 (MQ=255)
                                                                                                                 tGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAg            >  2:822186/1‑139 (MQ=255)
                                                                                                                    aTAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCa         <  2:386561/139‑1 (MQ=255)
                                                                                                                           cAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCa  <  1:902843/139‑1 (MQ=255)
                                                                                                                             aCGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCt                                                                  >  2:1523754/1‑73 (MQ=255)
                                                                                                                             aCGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCt                                                                  <  1:1523754/73‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                 
GGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCA  >  NC_000913/3140162‑3140423

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