Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,979,666 T→G 8.9% D322A (GAC→GCC) ‡ cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,979,6660TG 8.9% 132.4 / 4.5 45D322A (GAC→GCC) ‡cobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base G (2/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

GCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGG  >  NC_002947/3979526‑3979814
                                                                                                                                            |                                                                                                                                                     
gCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑cc                                                                                                                                              <  8:118096/149‑1 (MQ=255)
 cTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCATACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCg                                                                                                                                             >  5:188720/1‑149 (MQ=255)
  tCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCgg                                                                                                                                            <  6:285540/149‑1 (MQ=255)
   cGACCCATCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCTGGGTGTGCTCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAATCTGACCGGGCTCGTGATAATGCGCCCGTCCAGCT‑CCGgc                                                                                                                                           <  7:50928/149‑1 (MQ=255)
    gccccAGCGGCGGGCCAGCTCGGCGACAAACTCCATGCGCTCGGGGTCGGCCCCGGAGCAGACAACGACGGACCGGCACCGATCGCTGCGCCAGGCCATGTCCATGAAGCTGACCGGGCGCGTGATACTGCGCCCGTACAGCG‑CCGgc                                                                                                                                           >  8:284613/3‑148 (MQ=255)
      ccccccggccggcCAGCGCGCCGACCACCGCCATCCGCTCGGCGGGGGCGCGGTAGCAGACAACGTCGGCCTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGgcaat                                                                                                                                        <  4:43924/144‑1 (MQ=255)
           cggcggTCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGgcaatgcaat                                                                                                                                   <  7:153193/149‑1 (MQ=255)
                       tCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAg                                                                                                                       >  1:328013/1‑149 (MQ=255)
                          gggACATAGTCGATGGGCTGGGGGGGGGGGCGGTAGCAGTCAACGTCGGAGCGGCTGCGCTCGCTGCGGCAGGCCATGTCCTTGAAGGTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCATTGTGCAAGGCCAGGTc                                                                                                                    <  7:284613/147‑1 (MQ=255)
                             aCAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTcgcg                                                                                                                 >  5:167543/1‑149 (MQ=255)
                                 aGTCCATGCGCTCGGGGTGGGTGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTcgcgcgcg                                                                                                             <  5:67409/149‑1 (MQ=255)
                                   tCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGcc                                                                                                           >  7:237146/1‑149 (MQ=255)
                                       tGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAgg                                                                                                       >  5:205339/1‑149 (MQ=255)
                                       tGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAgg                                                                                                       <  6:205339/149‑1 (MQ=255)
                                           cTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCt                                                                                                   <  8:237146/149‑1 (MQ=255)
                                           cTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCt                                                                                                   <  2:222500/149‑1 (MQ=255)
                                                ggTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTc                                                                                               <  1:192220/148‑1 (MQ=255)
                                                ggTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTc                                                                                               >  2:192220/1‑148 (MQ=255)
                                                        cgGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTc                                                                                      <  8:100181/149‑1 (MQ=255)
                                                                gACAACGTCGGACTGGCTGCGCTCGCTGGGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagccag                                                                              >  1:186845/1‑149 (MQ=255)
                                                                 aCAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAgg                                                                                                                      <  2:128517/108‑1 (MQ=255)
                                                                 aCAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAgg                                                                                                                      >  1:128517/1‑108 (MQ=255)
                                                                   aaCGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCTTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCAAGGTCGCGCGCGCGCAGGCGCTGCTCGCTGGCTTACCAGCCAGGCt                                                                           >  6:154669/1‑149 (MQ=255)
                                                                               gCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTgcg                                                               >  6:14120/1‑149 (MQ=255)
                                                                                  gcgcTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGcc                                                                                                     >  5:104531/1‑108 (MQ=255)
                                                                                  gcgcTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGcc                                                                                                     <  6:104531/108‑1 (MQ=255)
                                                                                      tcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGTGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGAt                                                        >  5:310966/1‑149 (MQ=255)
                                                                                      tcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTccc                                                                                    <  2:189747/121‑1 (MQ=255)
                                                                                      tcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTccc                                                                                    >  1:189747/1‑121 (MQ=255)
                                                                                      tcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGAt                                                        >  6:137196/1‑149 (MQ=255)
                                                                                                     aTGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCgg                                         >  3:308250/1‑149 (MQ=255)
                                                                                                        tCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGTGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGa                                       <  6:310966/148‑1 (MQ=255)
                                                                                                         ccTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATg                                     >  6:324715/1‑149 (MQ=255)
                                                                                                          cTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCATGCAGCACCGGGAGGt                                    >  4:184725/1‑149 (MQ=255)
                                                                                                           ttGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTgg                                                                         >  4:142763/1‑111 (MQ=255)
                                                                                                           ttGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTgg                                                                         <  3:142763/111‑1 (MQ=255)
                                                                                                            tGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTcg                                  >  4:149652/1‑149 (MQ=255)
                                                                                                                 cTGACCGGGCGCGTGAT‑ATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGa                            <  3:275781/149‑1 (MQ=255)
                                                                                                                  tggacGGGCGCATTATCACGCGCCCGGTCAGCTGCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCgcggc                              >  1:211220/5‑148 (MQ=255)
                                                                                                                  tggacGGGCGCATTATCACGCGCCCGGTCAGCTGCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCgcggc                              >  1:206697/5‑148 (MQ=255)
                                                                                                                         gcgcGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGgc                     <  6:167543/149‑1 (MQ=255)
                                                                                                                                 tcacgcgcCCGGTCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTg             <  2:211220/145‑1 (MQ=255)
                                                                                                                                 tcacgcgcCCGGTCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTg             <  2:206697/145‑1 (MQ=255)
                                                                                                                                    tGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTgggg          >  7:119215/1‑149 (MQ=255)
                                                                                                                                           gTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCgg   >  6:253389/1‑149 (MQ=255)
                                                                                                                                            tCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCggg  >  5:311485/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
GCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCT‑CCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGG  >  NC_002947/3979526‑3979814

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