Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,065,589 T→C E59G (GAG→GGG)  cobU ← cobinamide kinase and cobinamide phosphate guanylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,065,5890TC97.9% 174.2 / NA 47E59G (GAG→GGG) cobUcobinamide kinase and cobinamide phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (29/17);  minor base A (1/0);  total (30/17)
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

TCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCTCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCT  >  NC_000913/2065443‑2065734
                                                                                                                                                  |                                                                                                                                                 
tCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCAc                                                                                                                                               >  2:488951/1‑151 (MQ=255)
  gCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACtg                                                                                                                                             >  2:61292/1‑151 (MQ=255)
   ccTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACtgt                                                                                                                                            <  2:51844/151‑1 (MQ=255)
            cccATTCATCAGGGTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTg                                                                                                                                   >  2:409549/1‑151 (MQ=255)
              cATTCATCAGGGTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCt                                                                                                                                 <  2:131776/151‑1 (MQ=255)
               aTTCATCAGGGTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTc                                                                                                                                >  2:14694/1‑151 (MQ=255)
                 tcatcaGGGTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCg                                                                                                                              >  1:316403/1‑151 (MQ=255)
                     caGGGTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGggc                                                                                                                           <  2:238575/150‑1 (MQ=255)
                      aGGGTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGggcgg                                                                                                                         >  2:493723/1‑151 (MQ=255)
                      aGGGTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGggcgg                                                                                                                         >  2:333707/1‑151 (MQ=255)
                         gTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCt                                                                                                                      <  2:448817/151‑1 (MQ=255)
                             ttATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGAATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGcc                                                                                                                  >  2:96233/1‑151 (MQ=255)
                                 cgccgcCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCgatg                                                                                                              <  1:409580/151‑1 (MQ=255)
                                 cgccgcCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCgatg                                                                                                              >  2:670354/1‑151 (MQ=255)
                                   ccgccATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCgatgat                                                                                                            <  2:421072/151‑1 (MQ=255)
                                     gccATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGt                                                                                                          >  2:70717/1‑151 (MQ=255)
                                         tAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCt                                                                                                       >  2:433992/1‑150 (MQ=255)
                                          aaTCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTAt                                                                                                     >  1:474667/1‑151 (MQ=255)
                                             cAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCg                                                                                                  >  2:597603/1‑151 (MQ=255)
                                             cAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCg                                                                                                  <  2:33670/151‑1 (MQ=255)
                                                    caGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGcc                                                                                           >  1:1962/1‑151 (MQ=255)
                                                      gATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCACCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCAt                                                                                         >  1:267123/1‑151 (MQ=255)
                                                          aGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCtca                                                                                     <  1:101555/151‑1 (MQ=255)
                                                           gTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCtcat                                                                                    <  1:675770/151‑1 (MQ=255)
                                                             cACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCtcatc                                                                                   <  2:278176/150‑1 (MQ=255)
                                                                 aTTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCtcatcatca                                                                               >  1:1514/1‑150 (MQ=255)
                                                                  ttGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCtcatcatca                                                                               >  1:280617/1‑149 (MQ=255)
                                                                    gtgGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAgg                                                                            >  1:588694/1‑150 (MQ=255)
                                                                    gtgGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGa                                                                           <  2:506115/151‑1 (MQ=255)
                                                                      ggTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGAtt                                                                         <  2:203377/151‑1 (MQ=255)
                                                                      ggTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGAt                                                                          <  2:32903/150‑1 (MQ=255)
                                                                           cGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCg                                                                     >  2:524673/1‑150 (MQ=255)
                                                                                          caACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAg                                                     <  2:11742/151‑1 (MQ=255)
                                                                                           aaCCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCGGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGa                                                    <  1:677204/151‑1 (MQ=255)
                                                                                            aCCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGa                                                    >  1:641094/1‑150 (MQ=255)
                                                                                                 aTTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCtgt                                              >  2:315887/1‑151 (MQ=255)
                                                                                                          aaTGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTcc                                     >  1:7659/1‑151 (MQ=255)
                                                                                                             gTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTcccc                                   >  2:341941/1‑150 (MQ=255)
                                                                                                             gTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCa                                  <  1:293037/151‑1 (MQ=255)
                                                                                                                tGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAAt                                >  1:418199/1‑150 (MQ=255)
                                                                                                                tGCATGAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAAt                                >  2:343062/1‑150 (MQ=255)
                                                                                                                     gAATTAATTCATCAAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAgcgc                          >  2:58851/1‑151 (MQ=255)
                                                                                                                               atcaAGATGTTGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTgg                >  2:565346/1‑151 (MQ=255)
                                                                                                                                         tGCCAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTcc      >  2:344175/1‑151 (MQ=255)
                                                                                                                                           ccAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGTCCTTGCCGATGAGGTTCTATCCGTGCAGCCATCGCATCATCAGGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGGGGCGACGCTTcc      >  1:572400/1‑149 (MQ=255)
                                                                                                                                            cAGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGc   >  1:256697/1‑151 (MQ=255)
                                                                                                                                             aGCGCCCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCt  <  2:177117/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCGCCTGATAATCCCATTCATCAGGGTCTTTATCGCCGCCATAATCAAACAACAGATTAGTCACCATTGTGGTAACGCATTCAAGCAACACAACCTCATTCGGGTTAATGTCTGCATGAATTAATTCATCAAGATGTTGCCAGCGCTCCACTGTGCGCCAGTGCTCCGGGCGGCCTTGCCGATGATGTTCTATCCGTGCAGCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCT  >  NC_000913/2065443‑2065734

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