Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,648,982 C→G 14.0% Q28H (CAG→CAC) ‡ ygjP ← metal‑dependent hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,648,9820CG14.0% 118.1 / 11.4 49Q28H (CAG→CAC) ‡ygjPmetal‑dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base C (22/20);  new base G (3/4);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

GCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGT  >  NC_002947/5648840‑5649119
                                                                                                                                              |                                                                                                                                         
gcgcACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAg                                                                                                                                      <  7:236837/148‑1 (MQ=255)
     cTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCACGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCAGGCGCAGGTGCgg                                                                                                                                <  3:19846/149‑1 (MQ=255)
      ttCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCggg                                                                                                                               >  2:196062/1‑149 (MQ=255)
        cGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCggg                                                                                                                               <  7:90971/147‑1 (MQ=255)
        cGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCAGGCGCAGGTGCggc                                                                                                                               <  8:31258/147‑2 (MQ=255)
        cGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCAGGCGCAGGTGCggc                                                                                                                               <  4:102663/147‑2 (MQ=255)
                   ttGGCGCTCATCGTGGGCATAGCCATACAGCGCCTTGTCGCTTTGAACGTCATGGCGGTCAGGGAACCGGGGCTGCAGGTATTCACCCTGGGGGGCGCTGTCGATCATCTCCCCAACCTGCGCGTGCTGGTGGGGGGGGTAGGCTTg                                                                                                                    <  6:45296/147‑1 (MQ=255)
                        gcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTaa                                                                                                             >  7:284691/1‑149 (MQ=255)
                        gcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTaa                                                                                                             >  5:6804/1‑149 (MQ=255)
                               cTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACgtagta                                                                                                      >  5:230278/1‑149 (MQ=255)
                                tGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTAc                                                                                                     >  4:334165/1‑149 (MQ=255)
                                   gCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGt                                                                                                  <  6:6804/149‑1 (MQ=255)
                                     aTAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCgg                                                                                                                                <  6:204093/117‑1 (MQ=255)
                                     aTAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCgg                                                                                                                                >  5:204093/1‑117 (MQ=255)
                                       aGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCtt                                                                                                                     >  3:51383/1‑126 (MQ=255)
                                       aGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCtt                                                                                                                     <  4:51383/126‑1 (MQ=255)
                                        gCCATACAGCGCCTTGGCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATgg                                                                                             <  3:204321/149‑1 (MQ=255)
                                         ccATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGACAGGGTAAAGGCGCTGCAGGTATTCACCCACGCTGTCGCTGTCCATCAACTGCCGCACCTGCGACTGCAGGTCCCGGGGGTAGTCTTGCAAGTAACTTAGTACCGTCATGgc                                                                                            >  5:173982/1‑149 (MQ=255)
                                                gcgcCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCggggg                                                                                                                             <  6:99663/109‑1 (MQ=255)
                                                gcgcCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCAGCGCCTGCAGGTGCggggg                                                                                                                             >  5:99663/1‑109 (MQ=255)
                                                  gcCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCggg                                                                                                                               >  6:286510/1‑105 (MQ=255)
                                                  gcCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCggg                                                                                                                               <  5:286510/105‑1 (MQ=255)
                                                   ccTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTTACGTAGTACGGTCATGGCGCTGTATTCg                                                                                  <  1:196062/149‑1 (MQ=255)
                                                   ccTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCg                                                                                  >  4:249227/1‑149 (MQ=255)
                                                        tCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCATCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACgta                                                                                                         <  3:330571/121‑1 (MQ=255)
                                                        tCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCATCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACgta                                                                                                         >  4:330571/1‑121 (MQ=255)
                                                          gCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTGATGGCGCTGTATTCGGCGTAAc                                                                           <  6:230278/149‑1 (MQ=255)
                                                              tGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGaa                                                                        <  3:249227/148‑1 (MQ=255)
                                                              tGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAAt                                                                       >  6:304772/1‑149 (MQ=255)
                                                                 aCGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGcc                                                                    >  3:48357/1‑149 (MQ=255)
                                                                      aTGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGAtt                                                                <  8:284691/148‑1 (MQ=255)
                                                                          cGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGc                                                           <  4:208794/149‑1 (MQ=255)
                                                                                                 gTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGgag                                     <  4:159709/148‑1 (MQ=255)
                                                                                                       acacAGCCGGTCGCAGTCGATCATCAGCCGCACCTGCGCCTGCAGGTGCGGGGGGATGGCTAGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGAAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTg                              >  8:153714/4‑149 (MQ=255)
                                                                                                          cAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCTGGTGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCg                           <  3:334165/149‑1 (MQ=255)
                                                                                                              cGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGcc                            <  4:116670/144‑1 (MQ=255)
                                                                                                              cGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGcc                            >  3:116670/1‑144 (MQ=255)
                                                                                                                gTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGc                     >  3:284027/1‑149 (MQ=255)
                                                                                                                         gATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGc                                           <  4:218453/118‑1 (MQ=255)
                                                                                                                         gATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGc                                           >  3:218453/1‑118 (MQ=255)
                                                                                                                         gATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAAtt                               >  6:119395/1‑130 (MQ=255)
                                                                                                                         gATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAAtt                               <  5:119395/130‑1 (MQ=255)
                                                                                                                              tgtGCCGGACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGATGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTc       <  6:173982/147‑1 (MQ=255)
                                                                                                                                 gCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcgcg    >  2:287644/1‑149 (MQ=255)
                                                                                                                                 cccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcgcg    >  5:76181/2‑149 (MQ=255)
                                                                                                                                 cccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcgcg    >  2:250286/2‑149 (MQ=255)
                                                                                                                                  ccGCACCTGCGCCTCCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGAAAGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCACAGGCCACCAGAAAATTGACAGGAACCCGCAAGAGTCAGTcgcc    >  2:325718/1‑147 (MQ=255)
                                                                                                                                  ccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCgg   >  7:167873/1‑149 (MQ=255)
                                                                                                                                   cGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGt  >  2:324310/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                         
GCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGT  >  NC_002947/5648840‑5649119

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