Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,873,937 G→C 12.5% G69A (GGC→GCC)  cobU → nicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,873,9370GC12.5% 92.7 / 4.7 32G69A (GGC→GCC) cobUnicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (14/14);  new base C (2/2);  total (16/16)
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

CGCCCGTCAGCAGCAACTGACCAAACCTGCTGGCTCGCTCGGCCAACTTGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCAGGCCAGCCTGGACGTGGTCGACCTGGGCACCATT  >  NC_002947/1873797‑1874079
                                                                                                                                            |                                                                                                                                              
cgcCCGTCAGCAGCAACTGACCAAACCTGCTGGCTCGCTCGGCCAACTTGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGCCGTGGTCgc                                                                                                                                       <  1:71778/150‑2 (MQ=255)
  cccGTCAGCAGCAACTGACCAAACCTGCTGGCTCGCTCGGCCAACTTGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGCCGTGGTCgccg                                                                                                                                     <  2:84702/150‑4 (MQ=255)
   ccGTCAGCAGCAACTGACCAAACCTGCTGGCTCGCTCGGCCAACTTGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAaga                                                                                                                                    <  1:3674/150‑1 (MQ=255)
         gcagcaACTGACCAAACCTGCTGGCTCGCTCGGCCAACTTGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCAt                                                                                                                              >  1:790323/1‑150 (MQ=255)
               aCTGACCAAACCTGCTGGCTCGCTCGGCCAACTTGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGc                                                                                                                        <  1:277350/150‑1 (MQ=255)
                    ccAAACCTGCTGGCTCGCTCGGCCAACTTGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATc                                                                                                                   >  2:595071/1‑150 (MQ=255)
                                gctcgctcGGCCAACTTGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCCATCCGCAGGCGGTaa                                                                                                       <  1:700975/150‑1 (MQ=255)
                                    gctcgGCCAACTTGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCgg                                                                                                   >  1:22046/1‑150 (MQ=255)
                                    gctcgGCCAACTTGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCgg                                                                                                   >  1:400553/1‑150 (MQ=255)
                                              cTTGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTa                                                                                         <  1:135072/150‑1 (MQ=255)
                                               ttGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTAc                                                                                        <  1:529161/150‑1 (MQ=255)
                                                    gcgcTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAg                                                                                                              <  1:460590/123‑1 (MQ=255)
                                                    gcgcTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAg                                                                                                              >  2:460590/1‑123 (MQ=255)
                                                       cTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTtc                                                                                >  1:596586/1‑150 (MQ=255)
                                                                        tggcCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCgcggc                                                               >  1:71534/1‑150 (MQ=255)
                                                                        tggcCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCgcggc                                                               >  1:570116/1‑150 (MQ=255)
                                                                             ttCAGGGGCTTGAACGACCCGCGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCTGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCa                                                          <  2:238239/150‑1 (MQ=255)
                                                                               cAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAgt                                                        <  2:400553/150‑1 (MQ=255)
                                                                               cAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAgt                                                        <  2:596586/150‑1 (MQ=255)
                                                                                   ggCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGc                                                    <  1:595071/150‑1 (MQ=255)
                                                                                          aCGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGcgcg                                             >  1:582174/1‑150 (MQ=255)
                                                                                                   gCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCa                                    >  2:162035/1‑150 (MQ=255)
                                                                                                        aTCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCAGGCCa                               >  1:483679/1‑150 (MQ=255)
                                                                                                              tCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCAg                                   <  1:659525/140‑1 (MQ=255)
                                                                                                              tCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCAg                                   >  2:659525/1‑140 (MQ=255)
                                                                                                                 ccATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCAGGCCAGCCTGGACg                      >  2:234510/1‑150 (MQ=255)
                                                                                                                      aCCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCAGGCCAGCCTGGACGTGGTc                 <  2:790323/150‑1 (MQ=255)
                                                                                                                          tttttGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCAGGCCAGCCTGGACGTGGTCGAcc             <  2:22046/150‑1 (MQ=255)
                                                                                                                              cctcttcGACCACGCCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCAGGCCAGCCTGGACGTGGTCGACCTggg         >  2:352881/7‑150 (MQ=255)
                                                                                                                               ctcttcGACCACGCCGTGGTCGAAGAGGGCATCACGTCATATCTGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCAGGCCAGCCTGGACGTGGTCGACCTGGGc        >  1:153963/6‑150 (MQ=255)
                                                                                                                                    ggACCACGGCGGGGTCGAAGGGGGCGGCGCGGCGTAGGCGCAGGGGGGAACCGGGCAGATGCGACGCAACGTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCAGGCCAGCCTGGACGTGGTCGACCTGGGCACCAt   <  2:483679/149‑1 (MQ=255)
                                                                                                                                     gACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCAGGCCAGCCTGGACGTGGTCGACCTGGGCACCAtt  >  1:676189/1‑150 (MQ=255)
                                                                                                                                            gCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTtcgt                                                                              >  2:596628/1‑67 (MQ=255)
                                                                                                                                            gCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTtcgt                                                                              <  1:596628/67‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
CGCCCGTCAGCAGCAACTGACCAAACCTGCTGGCTCGCTCGGCCAACTTGAAGCGCTGGCGATTCAGCTGGCTGGCCTTCAGGGGCTTGAACGACCCACGCTCGATCAGGTCGCCATCACCATTTTTGCCGGCGACCACGGCGTGGTCGAAGAGGGCATCTCGGCCTATCCGCAGGCGGTAACCGGGCAGATGCTACGCAACTTCGTCGGCGGGGGCGCGGCGATCAGTGTGCTGGCGCGGCAGCTGCAGGCCAGCCTGGACGTGGTCGACCTGGGCACCATT  >  NC_002947/1873797‑1874079

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: