Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,648,980 G→A 12.2% A29V (GCG→GTG) ‡ ygjP ← metal‑dependent hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,648,9800GA12.2% 122.5 / 6.7 49A29V (GCG→GTG) ‡ygjPmetal‑dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (3/3);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

CAGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGT  >  NC_002947/5648832‑5649119
                                                                                                                                                    |                                                                                                                                           
cAGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTgcg                                                                                                                                             <  7:278483/149‑1 (MQ=255)
        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)
                                               aGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTgcgc                                                                                                                                            <  1:207504/103‑1 (MQ=255)
                                               aGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTgcgc                                                                                                                                            >  2:207504/1‑103 (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                                                                                                         >  4:330571/1‑121 (MQ=255)
                                                                tCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCATCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACgta                                                                                                         <  3:330571/121‑1 (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                            >  3:116670/1‑144 (MQ=255)
                                                                                                                      cGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGcc                            <  4:116670/144‑1 (MQ=255)
                                                                                                                        gTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGc                     >  3:284027/1‑149 (MQ=255)
                                                                                                                                 gATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGc                                           >  3:218453/1‑118 (MQ=255)
                                                                                                                                 gATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGc                                           <  4:218453/118‑1 (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)
                                                                                                                                                    |                                                                                                                                           
CAGGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGT  >  NC_002947/5648832‑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: ‑