Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,140,351 C→T 12.5% R150H (CGC→CAC) ‡ hybE ← hydrogenase 2‑specific chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,140,3510CT12.5% 131.7 / 10.0 48R150H (CGC→CAC) ‡hybEhydrogenase 2‑specific chaperone
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base T (2/4);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

CGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAA  >  NC_000913/3140214‑3140495
                                                                                                                                          |                                                                                                                                                
cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑TCGCGCCCTGct                                                                                                                                       <  4:236576/149‑2 (MQ=255)
cgcgcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATACATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑TCGCGCCCTGct                                                                                                                                       <  3:325117/149‑2 (MQ=255)
   gcCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCgtgt                                                                                                                                    >  2:333201/1‑149 (MQ=255)
     cAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTgg                                                                                                                                  <  4:235437/149‑1 (MQ=255)
      aaTTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGt                                                                                                                                 >  1:42784/1‑149 (MQ=255)
                          tGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTgg                                                                                                             <  3:122014/149‑1 (MQ=255)
                            aCGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCTAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGGACATCCGGATTCGTGACTGGCa                                                                                                           >  1:408368/1‑149 (MQ=255)
                              gcgcTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGc                                                                                                         >  4:266788/1‑149 (MQ=255)
                                gcTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGa                                                                                                       >  3:403072/1‑149 (MQ=255)
                                gcTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGa                                                                                                       >  2:399923/1‑149 (MQ=255)
                                 cTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGaa                                                                                                      >  2:287478/1‑149 (MQ=255)
                                   tAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAg                                                                                                    <  1:333201/149‑1 (MQ=255)
                                      cATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGCAt                                                                                                 >  3:418920/1‑149 (MQ=255)
                                      cATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCAACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGAt                                                                                                 <  2:196130/149‑1 (MQ=255)
                                            gctgctCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCg                                                                                           >  1:460873/1‑149 (MQ=255)
                                                 tCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGTCG‑CGCCCTGct                                                                                                                                       <  2:352635/100‑2 (MQ=255)
                                                 tCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑TCGCGCCCTGct                                                                                                                                       >  1:352635/1‑99 (MQ=255)
                                                  ccGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCa                                                                                     <  2:300327/149‑1 (MQ=255)
                                                  ccGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCa                                                                                     >  1:291516/1‑149 (MQ=255)
                                                       tGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCAt                                                                                <  2:42784/149‑1 (MQ=255)
                                                                aTAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTgg                                                                                                                                  <  3:387009/90‑1 (MQ=255)
                                                                aTAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTgg                                                                                                                                  >  4:387009/1‑90 (MQ=255)
                                                                 tAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCg                                                                      >  3:205977/1‑149 (MQ=255)
                                                                        aaCGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGccc                                                               >  2:379011/1‑149 (MQ=255)
                                                                           ggCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCctct                                                            <  3:250194/149‑1 (MQ=255)
                                                                                           gaCAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCtgt                                            <  1:287478/149‑1 (MQ=255)
                                                                                             cAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGc                                          <  2:460873/149‑1 (MQ=255)
                                                                                              aaCTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCg                                         <  1:88286/149‑1 (MQ=255)
                                                                                                cTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCAt                                                                                              <  4:321910/94‑1 (MQ=255)
                                                                                                cTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCAt                                                                                              >  3:321910/1‑94 (MQ=255)
                                                                                                 tCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGaaa                                      >  3:445465/1‑149 (MQ=255)
                                                                                                 tCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGaaa                                      >  1:469920/1‑149 (MQ=255)
                                                                                                   aTGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgc                                    <  4:418920/149‑1 (MQ=255)
                                                                                                       aTTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgcggcg                                >  3:348993/1‑149 (MQ=255)
                                                                                                        ttttCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGa                               >  3:94385/1‑149 (MQ=255)
                                                                                                        ttttCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGa                               <  4:445465/149‑1 (MQ=255)
                                                                                                              ccACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATc                           >  2:179471/1‑147 (MQ=255)
                                                                                                              ccACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATc                           <  1:179471/147‑1 (MQ=255)
                                                                                                               cACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGc                        <  2:408368/149‑1 (MQ=255)
                                                                                                                aCTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCg                       <  3:19932/149‑1 (MQ=255)
                                                                                                                aCTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCg                       <  1:314165/149‑1 (MQ=255)
                                                                                                                  tCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGaa                     >  2:273610/1‑149 (MQ=255)
                                                                                                                       cGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGga                <  2:13632/149‑1 (MQ=255)
                                                                                                                            cacacgcagGGCG‑TCGCGCCCTGC‑TGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgcggc                                 <  2:432204/121‑1 (MQ=255)
                                                                                                                              cacgcagGGCG‑TCGCGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGga                >  2:48206/4‑142 (MQ=255)
                                                                                                                              cacgcagGGCG‑TCGCGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGga                <  1:48206/139‑1 (MQ=255)
                                                                                                                                     ggCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGaaa  <  2:469920/149‑1 (MQ=255)
                                                                                                                                     ggCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGaaa  >  2:46387/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                
CGCGCCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCG‑CGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTGCGAAA  >  NC_000913/3140214‑3140495

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: ‑