Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,780,140 G→A 8.6% intergenic (+57/+49) PP_1587 → / ← dapC NhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,780,1400GA 8.6% 126.0 / 2.8 47intergenic (+57/+49)PP_1587/dapCNhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (31/11);  new base A (2/2);  total (34/13)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

CTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTT  >  NC_002947/1779991‑1780287
                                                                                                                                                     |                                                                                                                                                    
cTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTg                                                                                                                                                      <  1:2212961/150‑1 (MQ=255)
        tCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgc                                                                                                                                              <  1:2522895/150‑1 (MQ=255)
         cAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcg                                                                                                                                             <  1:3359729/150‑1 (MQ=255)
           gCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCt                                                                                                                                           <  1:4681936/150‑1 (MQ=255)
           gCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACACGTAATGCGCTGCGCATGAAATCAACGCt                                                                                                                                           >  2:4312087/1‑150 (MQ=255)
            cACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTg                                                                                                                                          <  1:4033824/150‑1 (MQ=255)
               cTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTAc                                                                                                                                       <  1:930868/150‑1 (MQ=255)
               cTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCTCCCACACGTAAGTCGCGGCGCTTGAAATAAGCGCTGTAc                                                                                                                                       >  1:3686333/1‑150 (MQ=255)
                tggtggGGGATGTGGTGGTGCGGCTAGTACTGGGTGAGCTGGTTTTGAGTGAGGCGATCCTGGGTTCAGTCATGTAGACCCTTTCAAGGCCCTTGGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGATATCAGCGCTGTAcc                                                                                                                                      <  2:203192/150‑1 (MQ=255)
                   tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                                                  <  1:1557468/150‑1 (MQ=255)
                   tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                                                  <  1:553067/150‑1 (MQ=255)
                     ggggTTGTGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGTGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                                                <  1:3883468/150‑1 (MQ=255)
                     ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCCGGTACTGCGCGGCGGGTGGAATCAGCGCGGTGACTGTgg                                                                                                                                >  1:2401824/1‑150 (MQ=255)
                     ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACCGGGACAGCGCTGCGCTTGAAATCACCGCGGTAACGGTgg                                                                                                                                >  1:2401825/1‑150 (MQ=255)
                     ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTgg                                                                                                                                <  1:2035980/150‑1 (MQ=21)
                         aTGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGc                                                                                                                            <  1:2669386/150‑1 (MQ=255)
                           gAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCgg                                                                                                                          <  1:4352409/150‑1 (MQ=255)
                           gAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAGCgg                                                                                                                          <  2:441015/150‑1 (MQ=21)
                             ggtggtGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTAACAGTGGGAGCGGGc                                                                                                                        >  1:3313696/1‑150 (MQ=255)
                             ggtggtGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCAGTACCAGTGGGAGCGGGc                                                                                                                        >  1:3313594/1‑150 (MQ=255)
                                           tACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTACGCGGCCGCGCCCGCTCACAAAGGTACAGCGCAGCGCTTGAAATCAGAGCAGTACCAGTGGGAGCGGGCGTGGCCGCGAAggg                                                                                                          >  1:2976216/1‑150 (MQ=255)
                                                ggTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCGCCCACCAGTACAGCGCGGGGCGTGAAAAAAACCCGGTAACGGTGGGGGCGGGGGGGGCCGCGGAGGGGGTgt                                                                                                     >  1:1034453/1‑149 (MQ=255)
                                                      cTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCTCCCACAGGGACTGGGCAGCGCGTGTAATAAGCCCGGTAACGGTTGGGGGGGGGCGGGCCCCCGAGGGGGTGTTCAGCg                                                                                               >  1:2357828/1‑150 (MQ=255)
                                                         gATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACACGGACAGCGCTGCGCGTGAAAAAAGCGCAGTTCCAGTGGGGGCGGGCGTGCCCGCGAACGGGTTGGTCAGCGGtt                                                                                            >  1:1154511/1‑150 (MQ=255)
                                                                 tgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAATGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGa                                                                                    >  2:3511887/1‑150 (MQ=255)
                                                                     gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTTAAATCAGCGCTTTAACTTTGGGAGCGGGCGGGCCCGCGAACGGGTTGAACACCGGTTTTGCAGGAAGgc                                                                                >  2:3937483/1‑150 (MQ=38)
                                                                                        aGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCTCCCACAGGTAATGCCCAGCCCGTGCAATAAGCGCTGTACCAGTGGGAGCGGCCGTGCCCCCGAAGGGGTTGAACACCGGATTAGAAGGAAGGCCCGAAAGCGCAacgcctca                                                             >  2:329338/1‑147 (MQ=255)
                                                                                                  aTCAAGGCCCTTCGCGGGCGCGCCCGCTCCCAAAGGTACTGCGCTGCGCGTGAAATCCCCGCTGTACCGGTGGGAGCGGGCGTGCCCGCGGAGGGGGTGATCAGCCGTGTTGCGGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCAtt                                                   >  2:4978599/1‑150 (MQ=37)
                                                                                                                 gggCGCGCCCGCTCCCACACGTACTGCCCCGCGCGTGAAATACCCCCCGTACCAGGGGGAGCGGGCGTGCCCCCGCAGGGGGTGATCAGCCGGTTTTCAGGGAGGCGCGAAAGCGCACCGCCCCCTCGATGCATTCCGCCCCCGGTGCaa                                    >  2:1403564/1‑150 (MQ=255)
                                                                                                                     gcgcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAg                                >  2:3522682/1‑150 (MQ=255)
                                                                                                                     acgcCCGCTCCCACCGGTACAGCGCTGATTTCAAGCGCAGCCCAGTACCTGTGGGAGCCGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCCGGAAGGCGCGAATGCGCACCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAg                                >  2:4473450/2‑150 (MQ=18)
                                                                                                                     acgcCCGCTCCCACACGTACAGAGCTGATTTCAATCGCAGCGCAGTAACTGTGGGAGCGGGCGTGCCCGCGAATGGGTTGAACAGCGGTTTTTCAGGAATGCCCCAATGCGCTCCGCCGCTTCGATGCAATCGGCCACCGGTGCAACCaa                                >  2:4473106/2‑149 (MQ=19)
                                                                                                                         ccGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGcc                            >  2:301536/1‑150 (MQ=255)
                                                                                                                               ccACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgcgc                      >  1:2106983/1‑150 (MQ=255)
                                                                                                                                cacaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCa                     >  2:5035270/1‑150 (MQ=255)
                                                                                                                                cacaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCa                     >  2:874701/1‑150 (MQ=255)
                                                                                                                                 acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                    >  1:1121758/1‑150 (MQ=255)
                                                                                                                                 acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                    >  1:3826808/1‑150 (MQ=255)
                                                                                                                                 acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                    >  1:2996155/1‑150 (MQ=255)
                                                                                                                                 acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                    >  1:1317749/1‑150 (MQ=255)
                                                                                                                                   aGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAcgc                  >  1:1126461/1‑150 (MQ=255)
                                                                                                                                      tACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGcc               >  2:3754009/1‑150 (MQ=255)
                                                                                                                                         tgcgctgcgcTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCgg            >  2:2409102/1‑150 (MQ=255)
                                                                                                                                         tgcgctgcgcTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCgg            >  2:4001991/1‑150 (MQ=255)
                                                                                                                                         tgcgctgcgcTTGAAATCAGCGCTGTACCTGTGGGAGCGGGAGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCgg            >  2:4520175/1‑150 (MQ=255)
                                                                                                                                           cgctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCa          >  1:3960333/1‑150 (MQ=255)
                                                                                                                                              tgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAccc       >  1:2713673/1‑150 (MQ=255)
                                                                                                                                                   ttGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGtt  >  1:3462977/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTT  >  NC_002947/1779991‑1780287

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