Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,979,691 T→G 7.7% M314L (ATG→CTG)  cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,979,6910TG 7.7% 158.3 / 3.9 52M314L (ATG→CTG) cobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base T (26/22);  new base G (2/2);  total (28/24)
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.50e-01

GGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGT  >  NC_002947/3979551‑3979838
                                                                                                                                             |                                                                                                                                                   
ggCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAgg                                                                                                                                              <  1:231841/148‑1 (MQ=255)
 gCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTc                                                                                                                                            <  5:139064/149‑1 (MQ=255)
   gACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCAGGGcc                                                                                                                                                 >  7:293107/1‑142 (MQ=255)
   gACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCAGGGcc                                                                                                                                                 <  8:293107/142‑1 (MQ=255)
     cAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTcgcgc                                                                                                                                        >  2:35678/1‑149 (MQ=255)
            cATGCGCTCGGGGTTGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGccc                                                                                                                                  <  7:345859/148‑1 (MQ=255)
            cATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGcc                                                                                                                                   >  7:314958/1‑147 (MQ=255)
             aTGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAg                                                                                                                                >  5:219115/1‑149 (MQ=255)
             aTGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAg                                                                                                                                >  7:758/1‑149 (MQ=255)
                cgcTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGcc                                                                                                                             <  6:37205/149‑1 (MQ=255)
                   tCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTg                                                                                                                          >  3:211194/1‑149 (MQ=255)
                                    gCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagc                                                                                                         <  6:28106/149‑1 (MQ=255)
                                    gCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagc                                                                                                         >  1:167445/1‑149 (MQ=255)
                                      aGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGAGCCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagc                                                                                                         >  1:281441/1‑148 (MQ=255)
                                      aGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagcca                                                                                                       <  1:35678/149‑1 (MQ=255)
                                            cGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTgg                                                                                                 <  6:208043/149‑1 (MQ=255)
                                              tCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGtt                                                                                               <  6:219115/149‑1 (MQ=255)
                                                ggACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGt                                                                                             <  4:211194/149‑1 (MQ=255)
                                                 gACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCAGGGCCTGGGCGCGCGCGACCTGGCCATGcaca                                                                                                                      >  7:275642/1‑120 (MQ=38)
                                                 gACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCAGGGCCTGGGCGCGCGCGACCTGGCCATGcaca                                                                                                                      <  8:275642/123‑4 (MQ=38)
                                                    tGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTcc                                                                                                             <  5:41407/129‑1 (MQ=255)
                                                    tGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTcc                                                                                                             >  6:41407/1‑129 (MQ=255)
                                                                     ccAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTgca                                                                        <  2:167445/149‑1 (MQ=255)
                                                                         gCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCAc                                                                    >  6:315460/1‑149 (MQ=255)
                                                                                  ttGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTc                                                           >  4:162634/1‑149 (MQ=255)
                                                                                          gA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCgg                                                                 <  4:195504/135‑1 (MQ=255)
                                                                                          gA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCgg                                                                 >  3:195504/1‑135 (MQ=255)
                                                                                           ggacGGGCGCATTATCACGCGCCCGGTCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGa                                                    >  1:261877/4‑148 (MQ=255)
                                                                                              cGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACg                                                >  4:38446/1‑149 (MQ=255)
                                                                                                   gcgTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGgcgc                                           <  3:38446/149‑1 (MQ=255)
                                                                                                     gTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAg                                         >  7:352832/1‑149 (MQ=255)
                                                                                                         tAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTg                                     >  5:76007/1‑149 (MQ=255)
                                                                                                           aTGCGCCCGTCCAGCTCCGGCAATGCAATGTGCTTGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTggg                                   >  6:88828/1‑149 (MQ=255)
                                                                                                            tGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTgggg                                  >  6:263282/1‑149 (MQ=255)
                                                                                                             gcgcCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGgc                                 >  4:89149/1‑149 (MQ=255)
                                                                                                               gcCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCgcgg                                                       <  8:248624/125‑1 (MQ=255)
                                                                                                               gcCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCgcgg                                                       >  7:248624/1‑125 (MQ=255)
                                                                                                                cccGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCt                              >  4:330996/1‑149 (MQ=255)
                                                                                                                  cGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCg                            >  5:73708/1‑149 (MQ=255)
                                                                                                                    tCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCggg                          <  5:263282/149‑1 (MQ=255)
                                                                                                                     ccAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCgggg                         >  6:303670/1‑149 (MQ=255)
                                                                                                                      cAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGgcgc                                           >  7:324753/1‑130 (MQ=255)
                                                                                                                      cAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGgcgc                                           <  8:324753/130‑1 (MQ=255)
                                                                                                                       aGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAACCGGGGGCAGGTTGGGGCGCTCGGGGcg                       >  4:54247/1‑148 (MQ=255)
                                                                                                                            cGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGa                    <  8:758/147‑1 (MQ=255)
                                                                                                                            cGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACa                  <  5:315460/149‑1 (MQ=255)
                                                                                                                               caatgcaatgTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAgg               <  8:352832/149‑1 (MQ=255)
                                                                                                                                      atgTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGTCAAGGCGaaa          >  3:121326/1‑147 (MQ=255)
                                                                                                                                      atgTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGaaa          <  4:121326/147‑1 (MQ=255)
                                                                                                                                         tgCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCgg      <  2:261877/148‑1 (MQ=255)
                                                                                                                                           cATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCggtg    <  3:54247/148‑1 (MQ=255)
                                                                                                                                            aTGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCggtggt  >  7:330835/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGA‑CCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGT  >  NC_002947/3979551‑3979838

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: ‑