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,2940TC97.9% 165.2 / NA 47S9G (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/19);  minor base A (0/1);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCA  >  NC_000913/3140166‑3140423
                                                                                                                                |                                                                                                                                 
gACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAaga                                                                                                                         >  2:1055326/1‑139 (MQ=255)
 aCAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAagag                                                                                                                        <  1:1055326/139‑1 (MQ=255)
 aCAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAagag                                                                                                                        >  2:590094/1‑139 (MQ=255)
        cAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTc                                                                                                                 >  2:799803/1‑139 (MQ=255)
               aaaaCGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCAtt                                                                                                          <  2:520174/139‑1 (MQ=255)
                                    aaCGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCACTGGCAAAGAGCCAACTCATGCATTTTCGCCACTCCTGCGACCaaa                                                                                     <  1:590094/139‑1 (MQ=255)
                                         aGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCAtt                                                                                                          >  2:718448/1‑113 (MQ=255)
                                         aGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCAtt                                                                                                          <  1:718448/113‑1 (MQ=255)
                                           gagaGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGa                                                                                          >  2:610397/1‑127 (MQ=255)
                                           gagaGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGa                                                                                          <  1:610397/127‑1 (MQ=255)
                                                cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGcgcg                                                                         >  2:1174387/1‑139 (MQ=255)
                                                     cAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGcc                                                                    >  2:460946/1‑139 (MQ=255)
                                                        tttCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACaa                                                                                                                    <  2:924548/88‑1 (MQ=255)
                                                        tttCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACaa                                                                                                                    >  1:924548/1‑88 (MQ=255)
                                                                   acGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGcc                                                                    >  1:1203469/1‑125 (MQ=255)
                                                                   acGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGcc                                                                    <  2:1203469/125‑1 (MQ=255)
                                                                   acGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACa                                                      >  1:549749/1‑139 (MQ=255)
                                                                    cggcggTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCgtg                                                             <  1:708702/131‑1 (MQ=255)
                                                                    cggcggTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCgtg                                                             >  2:708702/1‑131 (MQ=255)
                                                                          tGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGATACATCCGGAt                                               <  1:13376/139‑1 (MQ=255)
                                                                                 cTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGt                                                         >  2:1346397/1‑122 (MQ=255)
                                                                                 cTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGt                                                         <  1:1346397/122‑1 (MQ=255)
                                                                                 cTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGAc                                        >  2:206129/1‑139 (MQ=255)
                                                                                    aaCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTgg                                     >  1:669/1‑139 (MQ=255)
                                                                                     aCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGc                                    <  1:1174387/139‑1 (MQ=255)
                                                                                     aCATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGc                                    >  1:753931/1‑139 (MQ=255)
                                                                                           tgctgcTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGAc                                                                                                                      >  1:467484/1‑51 (MQ=255)
                                                                                           tgctgcTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGAc                                                                                                                      <  2:467484/51‑1 (MQ=255)
                                                                                           tgctgcTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAAcagca                                                                                <  2:76585/89‑1 (MQ=255)
                                                                                           tgctgcTCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAAcagca                                                                                >  1:76585/1‑89 (MQ=255)
                                                                                              tgctCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAgg                           >  2:1123417/1‑139 (MQ=255)
                                                                                                ctCCGCCTGCCGTTGGATAATTTCAACGGCGccct                                                                                                                                 <  2:394882/35‑1 (MQ=255)
                                                                                                ctCCGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGAt                         >  2:595266/1‑139 (MQ=255)
                                                                                                   cGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCAt                                                                                                           >  1:520052/1‑54 (MQ=255)
                                                                                                   cGCCTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCAt                                                                                                           <  2:520052/54‑1 (MQ=255)
                                                                                                     ccTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGAc                                                                                                                      >  2:1333938/1‑41 (MQ=255)
                                                                                                     ccTGCCGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGAc                                                                                                                      <  1:1333938/41‑1 (MQ=255)
                                                                                                          cGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGcgcg                                                                         >  1:757610/1‑81 (MQ=255)
                                                                                                          cGTTGGATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGcgcg                                                                         <  2:757610/81‑1 (MQ=255)
                                                                                                               gATAATTTCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTc          >  1:258135/1‑139 (MQ=255)
                                                                                                                    tttCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCt                                                                  <  2:550798/78‑1 (MQ=255)
                                                                                                                    tttCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCt                                                                  >  1:550798/1‑78 (MQ=255)
                                                                                                                      tCAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTc   >  2:356101/1‑139 (MQ=255)
                                                                                                                       cAACGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCa  <  2:549749/139‑1 (MQ=255)
                                                                                                                         aCGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTgcgc              >  1:439062/1‑125 (MQ=255)
                                                                                                                         aCGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTgcgc              <  2:439062/125‑1 (MQ=255)
                                                                                                                          cGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGa                                                                                          <  1:961683/48‑1 (MQ=255)
                                                                                                                          cGGCGCCCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGa                                                                                          >  2:961683/1‑48 (MQ=255)
                                                                                                                                |                                                                                                                                 
GACAATTTCAAAACTAAAACGGACGGCGCTCTCCTCAACGCAGGAGAGCGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCA  >  NC_000913/3140166‑3140423

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