Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,864,440 G→T 24.9% R94L (CGC→CTC)  cobM → precorrin‑4 C(11)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,864,4400GT24.9% 37.0 / 3.1 16R94L (CGC→CTC) cobMprecorrin‑4 C(11)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (9/3);  new base T (2/2);  total (11/5)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

TCAACAGCGCCGAACTGCACCTGGAACAGATCATCGCCGCCATGCGCAGCGCACATGAAAAAGGTCAGGATGTGGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATG  >  NC_002947/3864311‑3864571
                                                                                                                                 |                                                                                                                                   
tCAACAGCGCCGAACTGCACCTGGAACAGATCATCGCCGCCATGCGCAGCGCACATGAAAAAGGTCAGGATGTGGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCTCCACCt                                                                                                                               >  1:284337/1‑136 (MQ=255)
tCAACAGCGCCGAACTGCACCTGGAACAGATCATCGCCGCCATGCGCAGCGCACATGAAAAAGGTCAGGATGTGGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCTCCACCt                                                                                                                               >  1:287046/1‑136 (MQ=255)
    cAGCGCCGAACTGCACCTGGAACAGATCATCGCCGCCATGCGCAGCGCACATGAAAAAGGTCAGGATGTGGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCaa                                                                                                                           <  1:469258/136‑1 (MQ=255)
           gAACTGCACCTGGAACAGATCATCGCCGCCATGCGCAGCGCACATGAAAAAGGTCAGGATGTGGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTgg                                                                                                                    <  1:84575/136‑1 (MQ=255)
               tGCACCTGGAACAGATCATCGCCGCCATGCGCAGCGCACATGAAAAAGGTCAGGATGTGGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCAt                                                                                                                >  1:427788/1‑136 (MQ=255)
                    cTGGAACAGATCATCGCCGCCATGCGCAGCGCACATGAAAAAGGTCAGGATGTGGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGAtt                                                                                                           >  2:497075/1‑136 (MQ=255)
                                                      aTGAAAAAGGTCAGGATGTGGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCagc                                                                         >  2:104074/1‑136 (MQ=255)
                                                              tgTCAGGATGTGGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGTGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGcgccc                                                                 >  2:344108/2‑135 (MQ=255)
                                                                     aTGTGGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCt                                                          >  1:447778/1‑136 (MQ=255)
                                                                       gtgGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCg                                                        >  1:439285/1‑136 (MQ=255)
                                                                       gtgGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCg                                                        >  1:49938/1‑136 (MQ=255)
                                                                                      aGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCTCCACCTGCAAGCGCTGGGCATCTATTACCAGATCATCCCCGGCGTTACCGCCACTGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGt                                         <  2:284337/136‑1 (MQ=255)
                                                                                      aGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCTCCACCTGCAAGCGCTGGGCATCTATTACCAGATCATCCCCGGCGTTACCGCCACTGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGt                                         <  2:287046/136‑1 (MQ=255)
                                                                                           cgACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCt                                    >  2:198754/1‑136 (MQ=255)
                                                                                                 cAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTg                              <  1:104074/136‑1 (MQ=255)
                                                                                                                             aTTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATg  >  1:306849/1‑136 (MQ=255)
                                                                                                                                 |                                                                                                                                   
TCAACAGCGCCGAACTGCACCTGGAACAGATCATCGCCGCCATGCGCAGCGCACATGAAAAAGGTCAGGATGTGGCCCGGGTGCACAGCGGCGACCCCAGCCTGTACGGCGCCATCGGTGAGCAGATTCGCCACCTGCAAGCGCTGGGCATCGATTACCAGATCATCCCCGGCGTTACCGCCACGGCAGCGAGCGCCGCCCTGCTCGGCTGCGAGCTGACGTTGCCTCAAGTGGCGCAGACCGTGATCCTCACCCGCTATG  >  NC_002947/3864311‑3864571

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