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,2940TC98.2% 190.7 / ‑3.7 54S9G (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 (27/26);  minor base G (0/1);  total (27/27)
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.83e-01

CCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGT  >  NC_000913/3140158‑3140428
                                                                                                                                        |                                                                                                                                      
ccGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTAACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGccct                                                                                                                                      >  2:857351/1‑139 (MQ=255)
   tGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGc                                                                                                                                   <  2:1391674/139‑1 (MQ=255)
          cAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAgaga                                                                                                                            <  2:1374982/139‑1 (MQ=255)
                              acggCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCa                                                                                                        >  2:296283/1‑139 (MQ=255)
                                            aaCGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCaaa                                                                                          >  2:950064/1‑139 (MQ=255)
                                                 aGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAAcagca                                                                                     >  2:1602481/1‑139 (MQ=255)
                                                 aGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAAcagca                                                                                     >  2:1310843/1‑139 (MQ=255)
                                                        cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCa                                                                                                                                  >  2:188187/1‑87 (MQ=255)
                                                        cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCa                                                                                                                                  <  1:188187/87‑1 (MQ=255)
                                                        cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAag                                                                                                                               <  1:1476531/90‑1 (MQ=255)
                                                        cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAag                                                                                                                               >  2:1476531/1‑90 (MQ=255)
                                                        cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGcgcg                                                                              <  2:560757/139‑1 (MQ=255)
                                                        cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGcgcg                                                                              >  1:560757/1‑139 (MQ=255)
                                                            ccAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGc                                                                          <  2:982080/139‑1 (MQ=255)
                                                            ccAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGc                                                                          <  1:296283/139‑1 (MQ=255)
                                                             cAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACt                                                                                                                       >  1:1404720/1‑93 (MQ=255)
                                                             cAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACt                                                                                                                       <  2:1404720/93‑1 (MQ=255)
                                                             cAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGcc                                                                         <  2:630073/139‑1 (MQ=255)
                                                               aTTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCt                                                                       <  1:950064/139‑1 (MQ=255)
                                                               aTTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCt                                                                       <  1:1602481/139‑1 (MQ=255)
                                                                  tCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCg                                                                    <  1:1565485/139‑1 (MQ=255)
                                                                          cacGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGcc                                                                         >  1:819778/1‑126 (MQ=255)
                                                                          cacGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGcc                                                                         <  2:819778/126‑1 (MQ=255)
                                                                           acGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATg                                                                                                                   >  1:1404059/1‑83 (MQ=255)
                                                                           acGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATg                                                                                                                   <  2:1404059/83‑1 (MQ=255)
                                                                                 gTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAgaga                                                                                                                            >  2:537454/1‑68 (MQ=255)
                                                                                 gTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAgaga                                                                                                                            <  1:537454/68‑1 (MQ=255)
                                                                                  tGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCa                                                                                                                                  >  2:725261/1‑61 (MQ=255)
                                                                                  tGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCa                                                                                                                                  <  1:725261/61‑1 (MQ=255)
                                                                                    aCGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGccc                                                                        >  1:1510752/1‑117 (MQ=255)
                                                                                    aCGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGccc                                                                        <  2:1510752/117‑1 (MQ=255)
                                                                                     cgcgCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGa                                                                                               >  1:1574879/1‑93 (MQ=255)
                                                                                     cgcgCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGa                                                                                               <  2:1574879/93‑1 (MQ=255)
                                                                                       cgcTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTg                                               >  2:1431750/1‑139 (MQ=255)
                                                                                         cTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAaga                                                                                                                              <  2:1435077/58‑1 (MQ=255)
                                                                                         cTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAaga                                                                                                                              >  1:1435077/1‑58 (MQ=255)
                                                                                         cTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCt                                                                       >  1:1212489/1‑113 (MQ=255)
                                                                                         cTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCt                                                                       <  2:1212489/113‑1 (MQ=255)
                                                                                           tAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGAGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTg                                           <  2:1018020/139‑1 (MQ=255)
                                                                                            aaCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCt                                                                       <  2:1603700/110‑1 (MQ=255)
                                                                                            aaCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCt                                                                       >  1:1603700/1‑110 (MQ=255)
                                                                                              cATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCa                                        <  1:857351/139‑1 (MQ=255)
                                                                                                 cGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCg                                     >  2:1548252/1‑139 (MQ=255)
                                                                                                   tgctgcTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTgg                                          <  1:865027/132‑1 (MQ=255)
                                                                                                   tgctgcTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTgg                                          >  2:865027/1‑132 (MQ=255)
                                                                                                        ctCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGAt                              >  2:106188/1‑139 (MQ=255)
                                                                                                        ctCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGACAGCGAAAGGAt                              >  2:467570/1‑139 (MQ=255)
                                                                                                                 ccGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTgc                     <  2:946112/139‑1 (MQ=255)
                                                                                                                       gATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCCGTc               >  1:287535/1‑139 (MQ=255)
                                                                                                                            tttCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGa                                                                                               <  1:1280652/54‑1 (MQ=255)
                                                                                                                            tttCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGa                                                                                               >  2:1280652/1‑54 (MQ=255)
                                                                                                                            tttCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCg          >  2:80103/1‑139 (MQ=255)
                                                                                                                                 aCGGCGCGCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTgcgc                   <  2:950771/125‑1 (MQ=255)
                                                                                                                                 aCGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGCACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTgcgc                   >  1:950771/1‑125 (MQ=255)
                                                                                                                                    gcgcCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGt  >  1:897091/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                      
CCGTGGCAGACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGT  >  NC_000913/3140158‑3140428

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