Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,977,355 G→C 8.4% H1092Q (CAC→CAG)  cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,977,3550GC 8.4% 144.4 / 2.3 47H1092Q (CAC→CAG) cobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base G (23/20);  new base C (2/2);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCTTGAGGGTACGGATGCGCGGTCGGTCGGCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGG  >  NC_002947/3977216‑3977490
                                                                                                                                           |                                                                                                                                       
ccTTGAGGGTACGGATGCGCGGTCGGTCGGCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTg                                                                                                                                 <  1:155486/148‑1 (MQ=255)
 cTTGAGGGTACGGATGCGCGGTCGGTCGGCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCACGACCTGCTCgatt                                                                                                                                <  1:270033/148‑4 (MQ=255)
  ttGAGGGTACGGATGCGCGGTCGGTCGGCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGtt                                                                                                                              >  8:255144/1‑149 (MQ=255)
          aCGGATGCGCGGTCGGTCGGCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGt                                                                                                                      >  4:268313/1‑149 (MQ=255)
             gATGCGCGGTCGGTCGGCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCt                                                                                                                   >  5:119237/1‑149 (MQ=255)
             gATGCGCGGTCGGTCGGCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGATCGTGGTTGTCCTGGTTCt                                                                                                                   >  4:39441/1‑149 (MQ=255)
               tGCGCGGTCGGTCGGCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGTGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTgc                                                                                                                 <  8:292133/149‑1 (MQ=255)
                  gcgGTCGGTCGGCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTg                                                                                                                        <  8:274386/139‑1 (MQ=255)
                  gcgGTCGGTCGGCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTg                                                                                                                        >  7:274386/1‑139 (MQ=255)
                        ggtcggCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGc                                                                                                         <  2:243135/148‑1 (MQ=255)
                             gCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCAc                                                                                                   >  7:289731/1‑149 (MQ=255)
                                  gCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCtt                                                                                                <  7:65728/147‑1 (MQ=255)
                                   cTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGc                                                                                             <  2:1936/149‑1 (MQ=255)
                                      tgGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGc                                                                                                     >  2:184224/1‑138 (MQ=255)
                                      tgGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGc                                                                                                     <  1:184224/138‑1 (MQ=255)
                                        gTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGGCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGgc                                                                                        >  2:241801/1‑149 (MQ=255)
                                         tCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTgc                                                                                                                 >  5:287618/1‑123 (MQ=255)
                                         tCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTgc                                                                                                                 <  6:287618/123‑1 (MQ=255)
                                                 ttAGCTGGCCACGAGCGCCCCGGACAGGGTTTCCGAGGCCGCCATCATGCCGCCCTGGAACTTGTAGTAATCGTTGGAATCGAGCATGTCGTGCTCGTGGTTGTGCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCg                                                                                <  3:39441/147‑1 (MQ=255)
                                                      tGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTgg                                                                          >  3:127269/1‑149 (MQ=255)
                                                          caccacCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTg                                                                                               >  1:287902/1‑124 (MQ=255)
                                                          caccacCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTg                                                                                               <  2:287902/124‑1 (MQ=255)
                                                           accaccGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCg                                                                     >  1:123003/1‑149 (MQ=255)
                                                             caccGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGAACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCggg                                                                   >  3:215933/1‑149 (MQ=255)
                                                                                 ccGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCt                                                <  2:16378/148‑1 (MQ=255)
                                                                                    aGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGCCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGGGCGGGCCGGGGTGCCGTCGTCGCTgccg                                            <  3:268313/149‑1 (MQ=255)
                                                                                         gccAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTa                                       <  6:119237/149‑1 (MQ=255)
                                                                                          ccAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAg                                      >  5:112996/1‑149 (MQ=255)
                                                                                           cAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAgg                                     <  2:209656/149‑1 (MQ=255)
                                                                                           cAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAgg                                     <  7:255144/149‑1 (MQ=255)
                                                                                                 gccgccCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAgc                               >  1:6807/1‑149 (MQ=255)
                                                                                                   cgccCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAgccg                             >  8:212424/1‑149 (MQ=255)
                                                                                                     cccTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAgccgcc                           <  3:94229/149‑1 (MQ=255)
                                                                                                               ggTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAgg                 >  4:5740/1‑149 (MQ=255)
                                                                                                                gtagtaATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGt                >  7:151739/1‑149 (MQ=255)
                                                                                                                     aaTCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGAcc           >  8:273747/1‑149 (MQ=255)
                                                                                                                     aaTCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGAc            <  2:294586/148‑1 (MQ=255)
                                                                                                                     aaTCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGAc            >  1:294586/1‑148 (MQ=255)
                                                                                                                            ggAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCgtcgtc                                                   <  6:86991/102‑1 (MQ=255)
                                                                                                                            ggAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCgtcgtc                                                   >  5:86991/1‑102 (MQ=255)
                                                                                                                              ccacGAGCACGACCTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAgg  >  3:64103/4‑149 (MQ=255)
                                                                                                                              ccacGAGCACGACCTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCCTAGCCATAGCCGCCGTGGt                      >  8:174170/4‑129 (MQ=255)
                                                                                                                              ccacGAGCACGACCTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGGCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGt                      <  7:174170/126‑1 (MQ=255)
                                                                                                                                  gAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCcgtc                                                      >  2:67811/1‑93 (MQ=255)
                                                                                                                                  gAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCcgtc                                                      <  1:67811/93‑1 (MQ=255)
                                                                                                                                   aGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAgg                 <  3:95252/129‑1 (MQ=255)
                                                                                                                                   aGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAgg                 >  4:95252/1‑129 (MQ=255)
                                                                                                                                           |                                                                                                                                       
CCTTGAGGGTACGGATGCGCGGTCGGTCGGCCTGGCTGTGGTCGCCATGGTAGCTGGCCACCACCGCCCCGGACAGGGTTTCCGAGGCCGCCAGCATGCCGCCCTGGAACTGGTAGTAATCGTTGGAATCGAGCAGGTCGTGCTCGTGGTTGTCCTGGTTCTGCAGCACCGCCTGCACCTTGGCCAGGCGCTGGGCGAACTGGGTGCGGGCCGGGGTGCCGTCGTCGCTGCCGCCGTAGGCATAGCCGCCGTGGTTGAGGTAGACCTCTGCCAGG  >  NC_002947/3977216‑3977490

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