Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,498,135 C→T 8.5% A403V (GCG→GTG) 
R2C (CGT→TGT) 		 cobL →
cbiD →		 precorrin‑6Y C(5,15)‑methyltransferase
cobalt‑precorrin‑6A synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,498,1350CT 8.5% 132.3 / 2.6 47A403V (GCG→GTG) 
R2C (CGT→TGT) 		cobL
cbiD		precorrin‑6Y C(5,15)‑methyltransferase
cobalt‑precorrin‑6A synthase
Reads supporting (aligned to +/- strand):  ref base C (22/21);  new base T (2/2);  total (24/23)
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.39e-01

ACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCACCCTGCCCAAAG  >  NC_002947/5497992‑5498276
                                                                                                                                               |                                                                                                                                              
aCTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTgaa                                                                                                                                           <  3:79154/149‑1 (MQ=255)
       cTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAAccc                                                                                                                                    >  7:67496/1‑149 (MQ=255)
          gCACAATTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCcgc                                                                                                                                  >  4:146403/1‑148 (MQ=255)
            acaTTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGaa                                                                                                                               >  4:272883/1‑149 (MQ=255)
            acaTTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGaa                                                                                                                               >  3:70214/1‑149 (MQ=255)
                          aGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCctct                                                                                                                 <  2:53236/149‑1 (MQ=255)
                                       gCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCtt                                                                                                    <  3:146403/149‑1 (MQ=255)
                                          gACCCGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ac                                                                                                 >  6:125825/1‑148 (MQ=255)
                                             ccGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCgcgc                                                                                                                     >  4:283244/1‑126 (MQ=255)
                                             ccGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCgcgc                                                                                                                     <  3:283244/126‑1 (MQ=255)
                                             ccGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑accacc                                                                                             >  6:257620/1‑149 (MQ=255)
                                                cATCCATGTCGCGCATGCCCAGCCCTTGGTTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGc                                                                                          <  8:67496/149‑1 (MQ=255)
                                                     aTGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTg                                                                                     <  6:129612/149‑1 (MQ=255)
                                                                  ccAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAAc                                                                                                                              >  1:85634/1‑96 (MQ=255)
                                                                  ccAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAAc                                                                                                                              <  2:85634/96‑1 (MQ=255)
                                                                  ccAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCctct                                                                                                                 >  7:284545/1‑109 (MQ=255)
                                                                  ccAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCctct                                                                                                                 <  8:284545/109‑1 (MQ=255)
                                                                  ccAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCgccaccgccacc                                                                        >  2:192608/1‑149 (MQ=255)
                                                                   cAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGccacca                                                                       <  3:272883/149‑1 (MQ=255)
                                                                        cTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCTTg                                                                  >  2:76/1‑149 (MQ=255)
                                                                                cATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGc                                                          >  4:170662/1‑149 (MQ=255)
                                                                                cATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGc                                                          >  7:98474/1‑149 (MQ=255)
                                                                                 aTTCGACACTTGGCGCCAGGCACTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGc                                                          >  2:245239/1‑148 (MQ=255)
                                                                                     gACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAgc                                                     <  5:257620/149‑1 (MQ=255)
                                                                                       cacTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCTTGGCCGCTGCCAGgctg                                                   <  1:76/149‑1 (MQ=255)
                                                                                          ttGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAgctgctg                                                >  6:323326/1‑149 (MQ=255)
                                                                                          ttGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAgctgctg                                                >  1:160056/1‑149 (MQ=255)
                                                                                                ccAGGCACTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTTCACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCAc                                           <  1:245239/149‑1 (MQ=255)
                                                                                                  aGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCgg                                        <  1:192608/149‑1 (MQ=255)
                                                                                                      gcTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGTGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAg                                    >  5:118390/1‑149 (MQ=255)
                                                                                                            gATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAAc                              <  3:170662/149‑1 (MQ=255)
                                                                                                            gATCACCCTTCTCGAAGTACTCAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGCCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAAc                              <  8:239474/149‑1 (MQ=255)
                                                                                                                 cccTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGc                         >  4:177452/1‑149 (MQ=255)
                                                                                                                         gAAGTACTGAAACCCGCCGATGCGTGAAGAACCCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGa                   <  5:323326/147‑1 (MQ=255)
                                                                                                                             tACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCAc             <  8:98474/149‑1 (MQ=255)
                                                                                                                               cTGAAACCCGCCGATGTGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCAcc            <  6:317382/148‑1 (MQ=255)
                                                                                                                               cTGAAACCCGCCGATGTGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCAcc            >  7:5393/1‑148 (MQ=255)
                                                                                                                               cTGAAACCCGCCGATGTGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCAcc            >  5:317382/1‑148 (MQ=255)
                                                                                                                               cTGAAACCCGCCGATGTGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCAcc            <  8:5393/148‑1 (MQ=255)
                                                                                                                                 gAAACCCGCCGATGCGTGAAGAAACCCGCGAACAACCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCACCCTg         >  7:200686/1‑149 (MQ=255)
                                                                                                                                  aaaCCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCACCCTGc        >  5:211882/1‑149 (MQ=255)
                                                                                                                                   aaCCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGTGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCACCCTGcc       <  6:118390/149‑1 (MQ=255)
                                                                                                                                       cgccgATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCACCCTGCCCa     >  8:64810/1‑147 (MQ=255)
                                                                                                                                        gccgATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGAt                           <  2:203241/124‑1 (MQ=255)
                                                                                                                                        gccgATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGAt                           >  1:203241/1‑124 (MQ=255)
                                                                                                                                        gccgATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCACCCTGCCCAAAg  <  5:125825/149‑1 (MQ=255)
                                                                                                                                        gccgATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCACCCTGCCCAAAg  <  4:70214/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
ACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTGGGTGCATTCGACACTTGGCGCCAGGCGCTGCCGATCACCCTTCTCGAAGTACTGAAACCCGCCGATGCGTGAAGAAACCCGCGAACAGCCCGCGCCTCTGCGCAGCGGCCTT‑ACCACCGGCAGCTGCGCCACCGCCACCAGCCTGGCAGCGGCCAAGCTGCTGCTCACCGGGCAGCGTAACGATGCGGTGGACATCACCCTGCCCAAAG  >  NC_002947/5497992‑5498276

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