Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,979,755 A→T 10.5% L307Q (CTG→CAG)  cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,979,7550AT10.5% 131.0 / 4.2 46L307Q (CTG→CAG) cobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base A (18/23);  new base T (3/2);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGACCAGCACTTCGGCGCCG  >  NC_002947/3979638‑3979904
                                                                                                                     |                                                                                                                                                     
aCGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTgg                                                                                                                       <  1:113043/150‑1 (MQ=255)
aCGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTgg                                                                                                                       <  2:664617/150‑1 (MQ=255)
 cggcagACTGGCTGCGCGCGCTGCGCCGGGTCATGTCGTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATCTCCAGGTCGCGCGCACCCAGCCCCTGCTCGCTGACTTGCCAGCCAGGCTGGt                                                                                                                      >  2:773042/6‑150 (MQ=255)
  gTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGtt                                                                                                                     <  2:4295/150‑1 (MQ=255)
      gACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTcc                                                                                                                 <  2:337397/150‑1 (MQ=255)
           gCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTgcgc                                                                                                            >  2:348736/1‑150 (MQ=255)
             tGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGcac                                                                                                          >  1:55099/1‑150 (MQ=255)
              gcgcTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGcaca                                                                                                         <  1:524475/150‑1 (MQ=255)
                               ccATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCg                                                                                        >  2:275496/1‑150 (MQ=255)
                               ccATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCg                                                                                        <  1:750235/150‑1 (MQ=255)
                                cATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCgg                                                                                       <  1:602558/150‑1 (MQ=255)
                                cATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCGTCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAACCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCgg                                                                                       <  1:485087/150‑1 (MQ=255)
                                      cTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagc                                                                                                                               >  2:191206/1‑104 (MQ=255)
                                      cTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagc                                                                                                                               <  1:191206/104‑1 (MQ=255)
                                      cTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTc                                                                                 >  2:73980/1‑150 (MQ=255)
                                       ttGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTcg                                                                                >  1:768517/1‑150 (MQ=255)
                                                 ccGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACg                                                                      >  1:166635/1‑150 (MQ=255)
                                                  cGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACgg                                                                     <  1:348736/150‑1 (MQ=255)
                                                  cGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACgg                                                                     <  1:451630/150‑1 (MQ=255)
                                                              aaTGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTggg                                                         <  1:73980/150‑1 (MQ=255)
                                                                tGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCGAGGTCGCGCGCGCCCAGGCCCTGCTCGCt                                                                                                                                          >  1:690860/1‑67 (MQ=255)
                                                                tGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCt                                                                                                                                          <  2:690860/67‑1 (MQ=255)
                                                                   gcCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCt                                                    >  1:357672/1‑150 (MQ=255)
                                                                     ccGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCg                                                  >  1:662061/1‑150 (MQ=255)
                                                                       gTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCggg                                                <  2:662061/150‑1 (MQ=255)
                                                                           aGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTg                                            <  2:55099/150‑1 (MQ=255)
                                                                                  gcaatgcaatGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGTAGGTTGGGGCGCTCGGGGCTGGACAAgg                                     <  1:582649/150‑1 (MQ=255)
                                                                                  gcaatgcaatGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAgg                                     >  2:363643/1‑150 (MQ=255)
                                                                                   caatgcaatgTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGc                                    <  2:768517/150‑1 (MQ=255)
                                                                                    aatgcaatgTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTggg                                                         <  2:540398/128‑1 (MQ=255)
                                                                                    aatgcaatgTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTggg                                                         >  1:540398/1‑128 (MQ=255)
                                                                                    aatgcaatgTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGCCg                                   >  2:830097/1‑150 (MQ=255)
                                                                                       gcaatgTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGaaa                                >  1:430901/1‑150 (MQ=255)
                                                                                        cagcgagCAGGGCCTGGGCGCGCGCGACCTGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAAc                               >  2:38388/7‑150 (MQ=255)
                                                                                             tgCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCggtg                          <  1:275496/150‑1 (MQ=255)
                                                                                             tgCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGAAGCACCGGGATGTCGGCGTGAATCATCCTCAGAGGGGGGCGCTCGGGGCTGGACAAGGCGAAACCggtg                          <  1:363643/150‑1 (MQ=255)
                                                                                               cATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCggtggt                        >  2:122982/1‑150 (MQ=255)
                                                                                                 gggCCTGGGCGCGCGCGACCTGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTAt                      >  2:407375/2‑150 (MQ=255)
                                                                                                     ctgggcgcgcgcgACCTGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAgg                                     <  2:474319/126‑1 (MQ=255)
                                                                                                     ctgggcgcgcgcgACCTGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAgg                                     >  1:474319/6‑131 (MQ=255)
                                                                                                     cAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGAc                  <  1:584860/150‑1 (MQ=255)
                                                                                                     cAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGAc                  >  1:150347/1‑150 (MQ=255)
                                                                                                         tCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTggg                                                         <  2:208557/107‑1 (MQ=255)
                                                                                                         tCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTggg                                                         >  1:208557/1‑107 (MQ=255)
                                                                                                           gcgcgcgaCCTGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGACCAGCAc            <  1:38388/150‑1 (MQ=255)
                                                                                                                 gcCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGACCAGCACTTCGgc      <  1:122982/150‑1 (MQ=255)
                                                                                                                     aGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGATCAGCACTTCGGCGCCg  >  2:424496/1‑150 (MQ=255)
                                                                                                                     |                                                                                                                                                     
ACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTTGGGGCGCTCGGGGCTGGACAAGGCGAAACCGGTGGTATTGACCAGCACTTCGGCGCCG  >  NC_002947/3979638‑3979904

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