Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,496,878 T→G 9.9% intergenic (‑194/‑50) cobG ← / → cobL precorrin‑3B synthase/precorrin‑6Y C(5,15)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,496,8780TG 9.9% 141.8 / 4.0 52intergenic (‑194/‑50)cobG/cobLprecorrin‑3B synthase/precorrin‑6Y C(5,15)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (23/22);  new base G (3/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTCCGGCGGCATGAAAAGCCGGTGGTCGAGCATCCTCGGCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCGCGCCCTGCTGGGCGCTGGCCGGGTC  >  NC_002947/5496732‑5497023
                                                                                                                                                  |                                                                                                                                                 
ttGTCCGGCGGCATGAAAAGCCGGTGGTCGAGCATCCTCGGCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTg                                                                                                                                                 >  7:250889/1‑149 (MQ=255)
 tGTCCGGCGGCATGAAAAGCCGGTGGTCGAGCATCCTCGGCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTg                                                                                                                                                 <  5:5520/148‑1 (MQ=255)
                      ggtggtCGAGCATCCTCGGCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCGCCCACCGGATCACTGCAGGCGGCACATGCTGCACTGTGCCTGTAGGAGCGGGCGTGCCCGCGAAGGGGCGCGACGCgg                                                                                                                           >  2:300989/1‑149 (MQ=38)
                       gtggtCGAGCATCCTCGGCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGc                                                                                                                          <  1:300989/149‑1 (MQ=255)
                           tCGAGCATCCTCGGCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGc                                                                                                                                                >  8:146017/1‑123 (MQ=255)
                           tCGAGCATCCTCGGCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGc                                                                                                                                                <  7:146017/123‑1 (MQ=255)
                           tCGAGCATCCTCGGCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGcccc                                                                                                                       >  8:148877/1‑148 (MQ=255)
                               gCATCCTCGGCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCaaaa                                                                                                                   >  2:339595/1‑148 (MQ=255)
                                   ccTCGGCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGc                                                                                                              >  4:111722/1‑149 (MQ=255)
                                       ggCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGgcgc                                                                                                                                  >  4:310396/1‑125 (MQ=255)
                                       ggCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGgcgc                                                                                                                                  <  3:310396/125‑1 (MQ=255)
                                           aTCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGaa                                                                                                      <  1:339595/149‑1 (MQ=255)
                                           aTCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGaa                                                                                                      <  1:45991/149‑1 (MQ=255)
                                            tCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGaaa                                                                                                     <  2:30924/149‑1 (MQ=255)
                                                gTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACa                                                                                                 <  7:148877/149‑1 (MQ=255)
                                                gTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACa                                                                                                 >  7:163590/1‑149 (MQ=255)
                                                   gAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGa                                                                                              <  3:111722/149‑1 (MQ=255)
                                                           ggACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTgg                                                                                      <  5:87123/149‑1 (MQ=255)
                                                               cgcgTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGa                                                                                  <  8:163590/149‑1 (MQ=255)
                                                                       cttcttCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAgtagta                                                                          <  8:250889/149‑1 (MQ=255)
                                                                          cttcACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGtagtag                                                                         <  6:299307/147‑1 (MQ=255)
                                                                                  gCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTg                                                                                       >  1:136169/1‑125 (MQ=255)
                                                                                   cTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCATGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCg                                                                >  6:12164/1‑147 (MQ=255)
                                                                                       cccGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACg                                                          >  7:199761/1‑149 (MQ=255)
                                                                                       cccGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACg                                                          >  2:208825/1‑149 (MQ=255)
                                                                                        ccGCTCCTACAGGTACAGTGCAGCATGTGCAGCCTGCAGTGATCCTGTGGGAGCGGGCGAGCCCGTGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACgg                                                         >  4:235640/1‑149 (MQ=17)
                                                                                          gCTCCTACAGGTACAGTGCAGCATGTGCAGCCTGCAGTGATCCTGTGGGAGCGGGCGAGCCCGCGAAGGGGCGCGAAGCGGCCCCATAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCt                                                       >  2:121977/1‑149 (MQ=17)
                                                                                           cTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCATGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCtt                                                      <  5:12164/149‑1 (MQ=255)
                                                                                              ccACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGCCAGTAGAAGGCAGCGGCGAAGACGGCTTCAg                                                   >  8:94527/1‑149 (MQ=255)
                                                                                                   ggATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGcc                                              >  4:24851/1‑149 (MQ=255)
                                                                                                     aTCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTg                                            <  7:94527/149‑1 (MQ=255)
                                                                                                      tCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTgg                                           <  1:208825/149‑1 (MQ=255)
                                                                                                        aCTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTggg                                          <  8:293176/148‑1 (MQ=255)
                                                                                                        aCTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGc                                         <  5:76326/149‑1 (MQ=255)
                                                                                                            cAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAAc                                     <  3:24851/149‑1 (MQ=255)
                                                                                                                 tGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGGCACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGcc                                >  2:228531/1‑149 (MQ=255)
                                                                                                                      aTGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCAGgcg                           >  8:8974/1‑149 (MQ=255)
                                                                                                                          tgcACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCGCGcc                       >  8:194379/1‑149 (MQ=255)
                                                                                                                          tgcACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCGCGcc                       >  3:63081/1‑149 (MQ=255)
                                                                                                                             aCTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCGCGCCctg                    >  8:278787/1‑149 (MQ=255)
                                                                                                                                 tACCTGTAGGAGCGGGCTTGCCCGCGAATGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGaa                                                              <  3:298900/103‑1 (MQ=255)
                                                                                                                                 tACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACCCCCTGGATGCCAGTAGTAGGCATCGGCGaa                                                              >  4:298900/1‑103 (MQ=255)
                                                                                                                                 tACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGa                                                                                  <  1:122524/83‑1 (MQ=255)
                                                                                                                                 tACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGa                                                                                  >  2:122524/1‑83 (MQ=255)
                                                                                                                                  tccTGTGGGAGCGGGCGAGCCCGTGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCGCGCCCTGCTggg               <  3:235640/148‑1 (MQ=255)
                                                                                                                                    cTGTGGGAGCGGGCGAGCCCGCGAAGGGGCGCGAAGCGGCCCCATAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCGCGCCCTGCTGGgc              <  1:121977/148‑1 (MQ=255)
                                                                                                                                    cTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCGCGCCCTGCTGGgcg             >  5:34468/1‑149 (MQ=255)
                                                                                                                                    cTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGAAAACAGGCCCGGCGCGCCCTGCTGGgcg             >  5:72766/1‑149 (MQ=255)
                                                                                                                                     tGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCGCGCCCTGCTGGgcgc            >  5:73779/1‑149 (MQ=255)
                                                                                                                                       tAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCGCGCCCTGCTgg                <  7:245562/143‑1 (MQ=255)
                                                                                                                                       tAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCGCGCCCTGCTgg                >  8:245562/1‑143 (MQ=255)
                                                                                                                                               ggCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCGCGCCCTGCTGGGCGCTGGCCGGGTc  <  7:278787/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTGTCCGGCGGCATGAAAAGCCGGTGGTCGAGCATCCTCGGCGATCAAGTAGAAATTCAGGACCGCGTCGCCTTCTTCACGGGCTCGCCCGCTCCCACAGGATCACTGCAGGCTGCACATGCTGCACTGTACCTGTAGGAGCGGGCTTGCCCGCGAAGGGGCGCGAAGCGGCCCCAAAAATGGCGGATTGAAATACATGACACCCTGGCTGACAGTAGTAGGCATCGGCGAAGACGGCTTCAGCGGCCTGGGCAAACAGGCCCGGCGCGCCCTGCTGGGCGCTGGCCGGGTC  >  NC_002947/5496732‑5497023

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