Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,780,129 T→C 12.8% intergenic (+88/+18) 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,1290TC12.8% 112.0 / 7.4 48intergenic (+88/+18)PP_1587/dapCNhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (25/15);  new base C (3/3);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTACG  >  NC_002947/1779981‑1780274
                                                                                                                                                    |                                                                                                                                                 
ctgagggcgcagcggttggagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgt                                                                                                                                                   <  4:2868‑M2/1‑1 (MQ=255)
                 ggagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGa                                                                                                                                  >  1:45639‑M2/132‑149 (MQ=255)
                       tGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCg                                                                                                                            >  7:288110/1‑149 (MQ=255)
                        ggATATGAGTGAGGCGAACCTGGGTTCAATCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAATGCGCTGCGCGTGAAATGAGCGCTGTAACGGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATAAGCgg                                                                                                                           >  1:206759/1‑149 (MQ=38)
                            acctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGtttt                                                                                                                       <  6:179535‑M2/29‑1 (MQ=255)
                              gagtgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCCTGCCCCCGAAGGGGTTGATCAGCGGTTTTg                                                                                                                      >  7:170311/1‑148 (MQ=255)
                                  atcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGa                                                                                                                 <  2:45639‑M2/35‑1 (MQ=255)
                                           cTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGCGCGGGAAAACAGCGGTGTACCTGTGGGCGCGGGCGTGCCCGCGAAGGGGGTGATCAGCGGTTTTGCAGGAAGGCGCCaa                                                                                                        >  4:117108/1‑149 (MQ=37)
                                              tgTTCAGTCAAGTCTACCCAATCAAGGCCCTTCGCGCGCGCGCCCGCTCCCCCAGGTCCGGCCCGGCGCTTTATATCAGCGCTGTACCTGGGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATg                                                                                                      <  5:193925/147‑1 (MQ=37)
                                                     tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTccgc                                                                                              >  7:205391/1‑149 (MQ=255)
                                                               gtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATg                                                                                    >  2:207685‑M2/86‑149 (MQ=255)
                                                               ccTATCAAGGCCCTTCGCGGGCGCGCCCGCGACCACAGGGAATGCGCAGCGGGTGAAAAGAGCCCAGTACCGGGGGGAGCGGGCCTGCCCGCGAAAGGGTGGATGAGCGGGTTTGGAGGAAGGCACGGATGCGCGGCGCCGCTTCGATg                                                                                    >  8:331246/1‑149 (MQ=255)
                                                                   tCAAGGCCCTTCGCGGGCGCGCCCGCTCCCCCAGGGACGGGGCGGGGGGGGGAATGAGCGGTGTGCCGGGGGGCGCGGGGGTGCCGGCGCACGGGGGGGGCACCCGTTGTGGAGGAAGGCGCGAGTGCGCAGCGCTCCGTCGATGCAtt                                                                                >  6:299419/1‑149 (MQ=255)
                                                                    cAAGGCCCTGCGGGGGCGCGCCCGCTCCCACAGTTAGTCCGTTGCGCTTTAATTCTGCGTTTTCCTTGTGGTTGGGGGCTTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTc                                                                               <  8:170311/149‑1 (MQ=255)
                                                                     caagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGtttt                                                                                                                       <  8:267360‑M2/29‑1 (MQ=255)
                                                                     caagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGtttt                                                                                                                       >  7:267360‑M2/80‑108 (MQ=255)
                                                                      aGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGGGCTTGAAATCAGCGGTGTACCGGGGGGAGCGGGCGGGCCCGCGAAGGGGTTGATCAGCGGGTTTGCAGGAAGGCGCGAATGCGCTgcggcgg                                                                                           >  5:187422/1‑134 (MQ=37)
                                                                               gcgGGCGCGCCCGCTCCCACTGGTACTGCGCTGCTCTGTAATTCAGCGGTGTACGTGGGGGCGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTg                                                                    <  8:288110/149‑1 (MQ=255)
                                                                                   gcgcgcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGGTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCGTCGATGCATTCGGCCAGCGGTGCAAc                                                                >  7:270250/1‑149 (MQ=255)
                                                                                        cccGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTg                                                           >  3:252334/1‑149 (MQ=255)
                                                                                          atgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCg                                                                       <  1:150485‑M2/77‑1 (MQ=255)
                                                                                           tgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGcc                                                         <  1:207685‑M2/91‑1 (MQ=255)
                                                                                               ccACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgc                                                     <  4:319296/148‑1 (MQ=255)
                                                                                               ccACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgc                                                     >  3:319296/1‑148 (MQ=255)
                                                                                                 aggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCGggcgcg                                                  >  3:176914‑M2/52‑148 (MQ=255)
                                                                                                  ggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                                                 >  3:279973‑M2/51‑149 (MQ=255)
                                                                                                     gaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAAc                                                                <  6:107284‑M2/84‑1 (MQ=255)
                                                                                                     gTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAcgcg                                              >  3:13054/1‑149 (MQ=255)
                                                                                                       acctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCa                                                        <  5:194050‑M2/92‑1 (MQ=255)
                                                                                                       acctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCa                                                        >  6:194050‑M2/46‑137 (MQ=255)
                                                                                                                 agtcaagtagaccctatcaaggcccttcgcgggcgCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCGGGCGCGTGCGCATggcactggt                                  <  4:176914‑M2/114‑5 (MQ=255)
                                                                                                                     gAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGTGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGtt                               >  5:363297/1‑148 (MQ=255)
                                                                                                                     gAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCa                                       >  5:118087/1‑140 (MQ=255)
                                                                                                                     gAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCa                                       <  6:118087/140‑1 (MQ=255)
                                                                                                                       aaTCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGATGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAc                            >  5:121975/1‑149 (MQ=255)
                                                                                                                        aTCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGTCCAGCGGTGCAACCAGTGACATGCGCACGCGCACGGCACCCGGGTTCAcc                           >  3:242740/1‑149 (MQ=255)
                                                                                                                        aTCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTAAcc                           >  3:20669/1‑149 (MQ=255)
                                                                                                                          cAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgcgc                                                   <  3:112958/123‑1 (MQ=255)
                                                                                                                          cAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgcgc                                                   >  4:112958/1‑123 (MQ=255)
                                                                                                                           aGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCa                        >  1:337488/1‑149 (MQ=255)
                                                                                                                                agTACCTGTGGGAGCGGGCGCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTACCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCAc                   >  5:157399/2‑149 (MQ=255)
                                                                                                                                  tcccggtgGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACtt                 <  8:270250/144‑1 (MQ=255)
                                                                                                                                       gtgGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcgc            >  4:92410/1‑149 (MQ=255)
                                                                                                                                        tgGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcg             <  5:190229/147‑1 (MQ=255)
                                                                                                                                         gggAGCGGGCGCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCa               >  3:32526/1‑144 (MQ=255)
                                                                                                                                         gggAGCGGGCGCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCa               <  4:32526/144‑1 (MQ=255)
                                                                                                                                          ggggggggcgGGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGAc         <  7:346949‑M2/139‑1 (MQ=255)
                                                                                                                                          ggAGCGGGCGCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGa          <  2:40528/148‑1 (MQ=255)
                                                                                                                                                 gcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTACg  >  4:199739‑M2/4‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTACG  >  NC_002947/1779981‑1780274

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