Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,979,760 C→A 9.8% Q305H (CAG→CAT cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,979,7600CA 9.8% 118.7 / 3.4 41Q305H (CAG→CATcobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base C (23/14);  new base A (2/2);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

GCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGACCAGCACTTCGGCGCCGACCTC  >  NC_002947/3979612‑3979909
                                                                                                                                                    |                                                                                                                                                      
gcTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCt                                                                                                                                                       <  2:281855/150‑1 (MQ=255)
                 tAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccag                                                                                                                                      >  1:321051/1‑150 (MQ=255)
                 aaGCAGACAACGTCGGACTGGCTGCGATAGCTGCGCCAGGCCATGTCCTTGAATCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGGAGTGCAATGTCCATTGCCAGGTCGCGCGCGCCTCGGCCCTGCTTGCTGGCTTCccag                                                                                                                                      >  1:260777/2‑150 (MQ=255)
                     aGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagccag                                                                                                                                  >  1:755721/1‑150 (MQ=255)
                         aaCGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTg                                                                                                                              <  2:764957/150‑1 (MQ=255)
                            gTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGtt                                                                                                                           <  1:837686/150‑1 (MQ=255)
                            gTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGtt                                                                                                                           >  2:324444/1‑150 (MQ=255)
                                   tGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCACGCTGGTTGTCCTgc                                                                                                                    >  2:396564/1‑150 (MQ=255)
                                       tGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGcac                                                                                                                >  2:974049/1‑150 (MQ=255)
                                          gctcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGa                                                                                                             >  2:952247/1‑150 (MQ=255)
                                           ctcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGTCTGGTTGTCCTGCGCACAGAt                                                                                                            <  1:197620/150‑1 (MQ=255)
                                            tcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATg                                                                                                           >  1:955096/1‑150 (MQ=255)
                                            tcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATg                                                                                                           >  1:315873/1‑150 (MQ=255)
                                              gctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGc                                                                                                         >  1:813234/1‑150 (MQ=255)
                                                           aTGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCggg                                                                                            >  2:98542/1‑150 (MQ=255)
                                                              tCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATg                                                                                         >  1:497404/1‑150 (MQ=255)
                                                                cTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTc                                                                                       <  2:755721/150‑1 (MQ=255)
                                                                cTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTc                                                                                       >  1:421396/1‑150 (MQ=255)
                                                                   gAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTcgcg                                                                                    <  1:98542/150‑1 (MQ=255)
                                                                     aGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCgcggc                                                                                  <  2:315873/150‑1 (MQ=255)
                                                                          aCCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGCCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAAc                                                                             <  1:1104115/150‑1 (MQ=255)
                                                                                 gcgTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCa                                                                      >  1:438347/1‑150 (MQ=255)
                                                                                     gATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGtt                                                                  >  1:975170/1‑150 (MQ=255)
                                                                                       tAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTgg                                                                <  2:438347/150‑1 (MQ=255)
                                                                                        aaTGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTggg                                                               >  1:155599/1‑150 (MQ=255)
                                                                                             gcCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCt                                                          >  1:1068683/1‑150 (MQ=255)
                                                                                             gcCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCt                                                          >  1:491210/1‑150 (MQ=255)
                                                                                                cGTCCGGCTCCGGCAATGCAATGTGCATGGCCATGTCGCGCGCGCCCCGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGATGTCCTGCGCCCAGACGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGCCGCAGGTTCGGGAGCTCgg                                                       >  2:449351/1‑150 (MQ=255)
                                                                                                     aGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTg                                                  <  2:491210/150‑1 (MQ=255)
                                                                                                        tCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGAc                                               <  2:421396/150‑1 (MQ=255)
                                                                                                           ggCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAg                                            <  2:955096/150‑1 (MQ=255)
                                                                                                                 gcaatgTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGaaa                                     >  1:255555/1‑150 (MQ=255)
                                                                                                                  cagcgagCAGGGCCTGGGCGCGCGCGACCTGGCCATGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAAc                                    >  2:885160/7‑150 (MQ=38)
                                                                                                                      gagcagggCCTGGGCGCGCGCGACCTGGCCATGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCggt                                >  1:139409/7‑150 (MQ=255)
                                                                                                                               cAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGAc                       >  1:7193/1‑150 (MQ=255)
                                                                                                                                 gggcgcgcgcgACCTGGCCATGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGACCa                     <  2:139409/147‑1 (MQ=255)
                                                                                                                                 ggTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGACCa                     <  2:7193/150‑1 (MQ=255)
                                                                                                                                  gTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGACCAg                    <  1:657803/150‑1 (MQ=255)
                                                                                                                                           gcCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAGACTGGACGAGTGACCAGCACTTCggg           >  2:527979/1‑149 (MQ=255)
                                                                                                                                           gacctggCCATGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGACCAGCACTTCGgc           <  1:521379/145‑1 (MQ=255)
                                                                                                                                             ccAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGTTGGCCTGCAGCACAGGGAGGTCGCGGCGACACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGACACCGGTGGTATTGACCAGCACTTTGgggc         >  2:424659/1‑150 (MQ=255)
                                                                                                                                                    cTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGACCAGCACTTCGGCGCCGACCTc  >  1:687978/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGACCAGCACTTCGGCGCCGACCTC  >  NC_002947/3979612‑3979909

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: