Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,300,642 T→C 10.8% D1044G (GAC→GGC)  recB ← Chi activated ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,300,6420TC10.8% 90.6 / 3.5 38D1044G (GAC→GGC) recBChi activated ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base T (16/16);  new base C (2/2);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

GGTCGGCTGCGCAGCTGGGCGCGCCAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTCGCTGAAG  >  NC_002947/5300494‑5300788
                                                                                                                                                     |                                                                                                                                                   
ggTCGGCTGCGCAGCTGGGCGCGCCAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGt                                                                                                                                                     <  3:271351/149‑1 (MQ=255)
    ggCTGCGCAGCTGGGCGCGCCAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCt                                                                                                                                                 >  2:39050/1‑149 (MQ=255)
             gCTGGGCGCGCCAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGCGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCa                                                                                                                                         <  1:39050/149‑1 (MQ=255)
                gggCGCGCCAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCgg                                                                                                                                     <  8:179840/149‑1 (MQ=255)
                 ggCGCGCCAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCTCGGCAGCgg                                                                                                                                     >  2:228868/1‑148 (MQ=255)
                 ggCGCGCCAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGACCGGCAGCgg                                                                                                                                     >  6:330015/1‑148 (MQ=255)
                    gcgcCAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAAc                                                                                                                                <  6:273573/149‑1 (MQ=255)
                     cgcCAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACg                                                                                                                               <  4:72147/149‑1 (MQ=255)
                       ccAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGcc                                                                                                                             <  4:36079/149‑1 (MQ=255)
                        cAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGcc                                                                                                                             >  3:136745/1‑148 (MQ=255)
                        cAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGcc                                                                                                                             <  4:136745/148‑1 (MQ=255)
                        cAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGCCGCTGCCGGGCAGCgacca                                                                                                                                  >  2:229615/1‑140 (MQ=255)
                        cAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACCCAGCTGGCCCAGCGTTACGCTCTGGCCGCTGCCGGGCAGCgacca                                                                                                                                  <  1:229615/144‑5 (MQ=255)
                           aCCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGGCCAGCGTTACGCTCTGGTCGCTGCTCGGCAGCGGCAACGCCTCGc                                                                                                                         <  1:228868/149‑1 (MQ=255)
                                    ggtgTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCTATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAAc                                                                                                                >  5:50484/1‑149 (MQ=255)
                                              gCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGCCGCTGCCGGGCAGCgaccaga                                                                                                                                >  5:294887/1‑118 (MQ=255)
                                              gCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGCCGCTGCCGGGCAGCgaccaga                                                                                                                                <  6:294887/124‑7 (MQ=255)
                                                 accaGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCTATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTg                                                                                                   <  6:50484/149‑1 (MQ=255)
                                                      gCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCt                                                                                                >  6:34274/1‑147 (MQ=255)
                                                                            aCCCGGTGGCTGGTGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCa                                                                        <  1:130745/149‑1 (MQ=255)
                                                                                    gctggcGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCgcggc                                                                >  1:117279/1‑149 (MQ=255)
                                                                                    gctggcGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCgcgg                                                                 <  1:280905/148‑1 (MQ=255)
                                                                                                    tCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGc                                                 >  5:252873/1‑148 (MQ=255)
                                                                                                      cATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCAGCGATTGCAGCGCTGGCCaa                                              >  8:35005/1‑149 (MQ=255)
                                                                                                          tCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGACGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACgg                                           >  2:164303/1‑148 (MQ=255)
                                                                                                           cAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGACGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGt                                          <  1:164303/148‑1 (MQ=255)
                                                                                                           cAATCTGGTAATGACCCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTa                                         >  3:101680/1‑149 (MQ=255)
                                                                                                                tGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGctc                                    >  2:25836/1‑149 (MQ=255)
                                                                                                                        aCGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGtt                            >  2:221258/1‑149 (MQ=255)
                                                                                                                         cGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGAGAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGtt                            >  1:228760/1‑149 (MQ=255)
                                                                                                                          gCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGTTGa                          <  2:228760/149‑1 (MQ=255)
                                                                                                                              cTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGTTGAGGGt                      >  5:93675/1‑149 (MQ=255)
                                                                                                                              cTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCAGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGTTGAGGGt                      <  7:35005/149‑1 (MQ=255)
                                                                                                                                       gCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTg              <  1:221258/149‑1 (MQ=255)
                                                                                                                                            aCGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTc         >  1:203115/1‑149 (MQ=255)
                                                                                                                                                tctGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGtt                            <  6:20143/126‑1 (MQ=255)
                                                                                                                                                tctGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGtt                            >  5:20143/1‑126 (MQ=255)
                                                                                                                                                    gTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTCGCTGaa   <  5:34274/147‑1 (MQ=255)
                                                                                                                                                    gTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTCGCTGAAg  <  2:117279/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGTCGGCTGCGCAGCTGGGCGCGCCAGACCGGGGTGGGTGTGGCGTGCCACCAGGCGATCGAGTTGTTCTGCGTCGACCCGGTGGCTGGCGAACCAGAACTCCATTTCAATCTGGTAATGACGCAGCTGGCCCAGCGTTACGCTCTGGTCGCTGCCCGGCAGCGGCAACGCCTCGCTCAGCAAACGCTGCATCCAGTGGCTCAGGGTGGGGATCCAGCCGGCCCAGTCGCGGCGATTGCAGCGCTGGCCAACGGTACGCTCTATCAGTTGAGGGTTAGCGCTGACCTCGCTGAAG  >  NC_002947/5300494‑5300788

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