Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,299,972 C→A 11.4% G60C (GGC→TGC)  recD ← multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,299,9720CA11.4% 101.4 / 2.3 35G60C (GGC→TGC) recDmultifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease
Reads supporting (aligned to +/- strand):  ref base C (14/17);  new base A (2/2);  total (16/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 = 9.36e-01

CAAAGGGCCGATAAGGGCGTCACCTTCCGGCGGCAACGACAGGGCGAAATCAGGCTCGGCCAAGGTTTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTGCCAGCGCCAGCAGGTGCTCTGCGTGCAACGGGGTGG  >  NC_002947/5299823‑5300121
                                                                                                                                                     |                                                                                                                                                     
cAAGGGGCCGATAAGGGCGTCACCTTCCGGCGGCAACGAAAGGGCGAAATCAGGCTCGGCCAAGGATTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTCGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGcc                                                                                                                                                       >  2:27920/1‑150 (MQ=255)
cAAAGGGCCGATAAGGGCGTCACCTTCCGGCGGCAACGACAGGGCGAAATCAGGCTCGGCCAAGGTTTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGcc                                                                                                                                                       >  1:419832/1‑150 (MQ=255)
        cGATAAGGGCGTCACCTTCCGGCGGCAACGACAGGGCGAAATCAGGCTCGGCCAAGGTTTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCAGGGGgcgc                                                                                                                                               >  1:135836/1‑150 (MQ=255)
            aaGGGCGTCACCTTCCGGCGGCAACGACAGGGCGAAATCAGGCTCGGCAAAGGTTTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGctc                                                                                                                                           <  1:317261/150‑1 (MQ=255)
                 cggcACCTGCCGGCGGCAACGACAGGGCGAAATCAGGCTCGGCCAAGGTTTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCa                                                                                                                                      <  2:342788/147‑1 (MQ=255)
                     aCCTTCCGGCGGCAACGACAGGGCGAAATCAGGCTCGGCCAAGGTTTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGaa                                                                                                                                  >  1:326019/1‑150 (MQ=255)
                          ccGGCGGCAACGACAGGGCGAAATCAGGCTCGGCCAAGGTTTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGaca                                                                                                                             <  1:27544/150‑1 (MQ=255)
                                  aaCGACAGGGCGAAATCAGGCTCGGCCAAGGTTTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATgcg                                                                                                                     >  1:200300/1‑150 (MQ=255)
                                          ggCGAAATCAGGCTCGGCCAAGGTTTGTTGCAGGTCGAGGCAGACCCGTGCCAGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAg                                                                                                             <  1:178240/150‑1 (MQ=255)
                                            cGAAATCAGGCTCGGCCAAGGTTTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGTAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAggg                                                                                                           <  1:507297/150‑1 (MQ=255)
                                                 tCAGGCTCGGCCAAGGTTTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGt                                                                                                      <  2:221552/150‑1 (MQ=255)
                                                     gCTCGGCCAAGGTTTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTagcc                                                                                                  >  2:341475/1‑150 (MQ=255)
                                                                 tttGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCaa                                                                                      >  1:101187/1‑150 (MQ=255)
                                                                 tttGTTGCAGGTCGAGGCAGACATGTCCGTGCCACAGTTGGTGGCTGGCCAGCTCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGAGACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCaa                                                                                      <  1:539883/150‑1 (MQ=255)
                                                                  ttgttgCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAAc                                                                                     <  2:252401/150‑1 (MQ=255)
                                                                          ggTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCAGGGGGCGCGCTCTTCCAGGAACGACACGAAt                                                                                                                        <  2:293848/107‑1 (MQ=255)
                                                                          ggTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCAGGGGGCGCGCTCTTCCAGGAACGACACGAAt                                                                                                                        >  1:293848/1‑107 (MQ=255)
                                                                             cGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCAGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGa                                                                          <  2:135836/150‑1 (MQ=255)
                                                                                    gACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGa                                                                          <  1:184698/143‑1 (MQ=255)
                                                                                    gACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGa                                                                          >  2:184698/1‑143 (MQ=255)
                                                                                                  ccAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGCAATTGCAGCAGGTCAgcgc                                                     >  1:493932/1‑150 (MQ=255)
                                                                                                    aGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGTCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCtg                                                   <  2:380402/150‑1 (MQ=255)
                                                                                                      ttGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCtgtt                                                 >  1:157740/1‑150 (MQ=255)
                                                                                                         gtggCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTgcagc                                                               >  1:337877/1‑133 (MQ=255)
                                                                                                         gtggCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTgcagc                                                               <  2:337877/133‑1 (MQ=255)
                                                                                                              tggcCAGCGCTGCCGCCAGCAATAGCAAGGGATCCCTGCCGGGGGCGAGCTCTTCCAGGAACGACACGAAGGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGGAATTGCAGCAGGTCAGCGCTGTTGTGCTGCg                                         <  2:206345/150‑1 (MQ=255)
                                                                                                                  cAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTgc                                     >  1:290665/1‑150 (MQ=255)
                                                                                                                    gcgcTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGTGCGCTCTGCCAGGAACGACACGAATGCGCGGTCTAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTgcca                                   <  1:320141/150‑1 (MQ=255)
                                                                                                                       cTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTgccagcg                                >  2:210330/1‑150 (MQ=255)
                                                                                                                       cTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTgccagcg                                >  2:360017/1‑150 (MQ=255)
                                                                                                                             ccAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCTACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTGCCAGCGCcagca                          <  2:137426/150‑1 (MQ=255)
                                                                                                                               aGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTGCCAGCGCCAGCAgg                        >  2:441618/1‑150 (MQ=255)
                                                                                                                                 cAATAGCAAGGGATCACTGCCGGGGGCGCGCTCGTCCAGGAAAGCTTTGGTCACCCTGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTGCCAGCGCCAGCAGGTg                      <  2:158490/150‑1 (MQ=255)
                                                                                                                                         aGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTGCCAGCGCCAGCAGGTGCTCTGCgg              >  1:198732/1‑149 (MQ=255)
                                                                                                                                         aGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAAAGCCAGCAGGTCAGCGCTGTTGTGCTGCGGTGCCAGCGCCAGCAGGTGCTCTGCGt              <  2:419832/150‑1 (MQ=255)
                                                                                                                                          gggATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTGCCAGCGCCAGCAGGTGCTCTGCGTg             <  1:341475/150‑1 (MQ=255)
                                                                                                                                                aCTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACCGCGAAGTGATTCCAGCAGGTCAGCGCTGTTGTGCTGCAGTCCGACCCCGAGCCGTTGCTCTGCGTGCAACgg       >  1:9986/1‑150 (MQ=255)
                                                                                                                                                  tgacGGGGGCGCGATCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCGACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTGCCAGCGCCAGCAGGTGCTCTGCGTGCAACgggg     >  2:317910/4‑150 (MQ=255)
                                                                                                                                                     cGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTGCCAGCGCCAGCAGGTGCTCTGCGTGCAACGGGGTgg  <  1:210330/150‑1 (MQ=255)
                                                                                                                                                     cGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTGCCAGCGCCAGCAGGTGCTCTGCGTGCAACGGGGTgg  <  1:524158/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAAGGGCCGATAAGGGCGTCACCTTCCGGCGGCAACGACAGGGCGAAATCAGGCTCGGCCAAGGTTTGTTGCAGGTCGAGGCAGACATGTCCGTGCCCCAGTTGGTGGCTGGCCAGCGCTGCCGCCAGCAATAGCAAGGGATCACTGCCGGGGGCGCGCTCTTCCAGGAACGACACGAATGCGCGGTCCAGGGCACGTAGCCAGCCCCGCTCAACCCAACGGTCGAGTAATTGCAGCAGGTCAGCGCTGTTGTGCTGCGGTGCCAGCGCCAGCAGGTGCTCTGCGTGCAACGGGGTGG  >  NC_002947/5299823‑5300121

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