Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,081,788 G→A 7.8% P50P (CCC→CCT)  dpkA ← delta 1‑piperideine‑2‑carboxylate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,081,7880GA 7.8% 143.6 / 2.7 51P50P (CCC→CCT) dpkAdelta 1‑piperideine‑2‑carboxylate reductase
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base A (2/2);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CGCCACCAACAGCTCACGGGCCGCCGCCAGTGCCGGCTGGGCAAAACCGCCCGCAGCATCGACGCGCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCTTTGCAGCTCGGTAAAAGGCACGCGCA  >  NC_002947/4081641‑4081936
                                                                                                                                                   |                                                                                                                                                    
cGCCACCAACAGCTCACGGGCCGCCGCCAGTGCCGGCTGGGCAAAACCGCCCGCAGCATCGACGCGCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCg                                                                                                                                                      <  1:317301/148‑1 (MQ=255)
   caccaACAGCTCACGGGCCGCCGCCAGTGCCGGCTGGGCAAAACCGCCCGCAGCATCGACGCGCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCa                                                                                                                                                  >  2:44670/1‑149 (MQ=255)
                         gccAGTGCCGGCTGGGCAAAACCGCCCGCAGCATCGACGCGCACATAGCCGGCGGCCACGTAGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCAAGCGGAACACCCCATGGCTATg                                                                                                                            >  3:250669/1‑149 (MQ=255)
                                     tGGGCAAAACCGCCCGCAGCATCGACGCGCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGcc                                                                                                                      <  1:221730/143‑1 (MQ=255)
                                     tGGGCAAAACCGCCCGCAGCATCGACGCGCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGcc                                                                                                                      >  2:221730/1‑143 (MQ=255)
                                     tGGGCAAAACCGCCCGCAGCATCGACGCGCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATcgcg                                                                                                                >  3:63649/1‑149 (MQ=255)
                                         cAAAACCGCCCGCAGCATCGACGCGCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCgcgctgg                                                                                                            <  5:239177/149‑1 (MQ=255)
                                               cgcccgcAGCATCGACGCGCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCTAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCgcgctgggcgctg                                                                                                      <  2:285845/149‑1 (MQ=255)
                                                             aCGCGCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGGCCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGc                                                                                                                                >  1:84000/1‑109 (MQ=255)
                                                             aCGCGCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGGCCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGc                                                                                                                                <  2:84000/109‑1 (MQ=255)
                                                             aCGCGCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGc                                                                                        <  2:317289/149‑1 (MQ=255)
                                                              cgcgCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCgcgc                                                                                                               <  8:202057/125‑1 (MQ=255)
                                                              cgcgCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCgcgc                                                                                                               >  7:202057/1‑125 (MQ=255)
                                                               gcgcACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGc                                                                                        >  6:206555/1‑147 (MQ=255)
                                                               gcgcACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGc                                                                                        >  1:263958/1‑147 (MQ=255)
                                                                 gcACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGc                                                                                    >  5:270286/1‑149 (MQ=255)
                                                                  cacaTAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGcc                                                                                       >  3:149618/1‑145 (MQ=255)
                                                                  cacaTAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGcc                                                                                       <  4:149618/145‑1 (MQ=255)
                                                                     aTAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCAccc                                                                                >  5:111254/1‑149 (MQ=255)
                                                                       aGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCAGGCGCAGTTGTGGGCCAGCACCCgg                                                                              >  7:261979/1‑149 (MQ=255)
                                                                        gCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCgg                                                                              <  5:206555/148‑1 (MQ=255)
                                                                            gcggcCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGc                                                                            <  3:189455/146‑1 (MQ=255)
                                                                            gcggcCACGTCGCTGACCTCTGGCGTGGCCTGCCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGCGCATGCGCAACACCCCATCGCTATGGGCGCCATCGCGCAGGGCGCTCGCGCAGATGTGCGCCAGCACCCGGGc                                                                            >  4:189455/1‑146 (MQ=255)
                                                                                 cACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCt                                                                    >  2:124277/1‑149 (MQ=255)
                                                                                 cACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGAGGGCCAGCACCCGGGCCACGGCCt                                                                    >  3:217548/1‑149 (MQ=255)
                                                                                           aCCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCAc                                                            <  6:270286/147‑1 (MQ=255)
                                                                                                  gcgtggcCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCgcgctg                                                                                                             >  4:116377/1‑91 (MQ=255)
                                                                                                  gcgtggcCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCgcgctg                                                                                                             <  3:116377/91‑1 (MQ=255)
                                                                                                  gcgtggcCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGgcgc                                                   >  3:291627/1‑149 (MQ=255)
                                                                                                   cgtggcCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGAGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCt                                                  <  4:217548/149‑1 (MQ=255)
                                                                                                       gcctggccATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGa                                               >  3:188738/1‑148 (MQ=255)
                                                                                                             ccATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGc                                        <  4:291627/149‑1 (MQ=255)
                                                                                                                 cGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCAc                                                                         <  8:110918/112‑1 (MQ=255)
                                                                                                                 cGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCAc                                                                         >  7:110918/1‑112 (MQ=255)
                                                                                                                          cGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCtt                           >  1:171255/1‑149 (MQ=255)
                                                                                                                          cGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCtt                           >  6:46381/1‑149 (MQ=255)
                                                                                                                           gCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCttt                          <  2:171255/149‑1 (MQ=255)
                                                                                                                            cTGGCCAAGGTCGAGACATAACCGGGGATGCGGAACTGCCCATGGCTATGGGCGCCTTCGCGGTGGGCGCTGGCGCAGTTGTTGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCTGTGGAAAATGGCCTGCAACAGGCTTTg                         <  4:250669/149‑1 (MQ=255)
                                                                                                                             tGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGTCCTGCAACa                                <  6:278160/141‑1 (MQ=255)
                                                                                                                             tGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGCCCTGCAACa                                >  5:278160/1‑141 (MQ=255)
                                                                                                                               gCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCTTTGc                        <  6:111254/147‑1 (MQ=255)
                                                                                                                                  aaGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCACCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCTTTGCAGCTc                   >  2:89870/1‑149 (MQ=255)
                                                                                                                                      tCGAGACATAACCAGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGTCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATg                                          >  5:309568/1‑122 (MQ=255)
                                                                                                                                      tCGAGACATAACCAGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGTCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATg                                          <  6:314978/122‑1 (MQ=255)
                                                                                                                                      tCGAGACATAACCAGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGTCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATg                                          <  6:309568/122‑1 (MQ=255)
                                                                                                                                      tCGAGACATAACCAGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGTCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATg                                          >  5:314978/1‑122 (MQ=255)
                                                                                                                                         agaCATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCTTTGCAGCTCGGTaaaa            <  4:188738/149‑1 (MQ=255)
                                                                                                                                          gaCATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTCTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAAAGGCCTGc                                    >  1:222058/1‑124 (MQ=255)
                                                                                                                                            cATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCTTTGCAGCTCGGTAAAAGGc         <  2:263958/149‑1 (MQ=255)
                                                                                                                                              tAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCTTTGCAGCTCGGTAAAAGGCAc       >  3:95279/1‑149 (MQ=255)
                                                                                                                                                  cGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCTTTGCAGCTCGGTaaaa            >  3:227337/1‑140 (MQ=255)
                                                                                                                                                  cGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCTTTGCAGCTCGGTaaaa            <  4:227337/140‑1 (MQ=255)
                                                                                                                                                   gggCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCTTTGCAGCTCGGTAAAAGGCACGCGca  <  1:137581/149‑1 (MQ=255)
                                                                                                                                                   gggCATGCGGAACACCCCATGGCTATGGGCGCCACCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCTTTGCAGCTCGGTAAAAGGCACGCGca  <  1:89870/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCCACCAACAGCTCACGGGCCGCCGCCAGTGCCGGCTGGGCAAAACCGCCCGCAGCATCGACGCGCACATAGCCGGCGGCCACGTCGCTGACCTGTGGCGTGGCCTGGCCATCGACCCAGCCGCTGGCCAAGGTCGAGACATAACCGGGCATGCGGAACACCCCATGGCTATGGGCGCCATCGCGCTGGGCGCTGGCGCAGTTGTGGGCCAGCACCCGGGCCACGGCCTCGCTGCACCCATGGCGCTGGAAAATGGCCTGCAACAGGCTTTGCAGCTCGGTAAAAGGCACGCGCA  >  NC_002947/4081641‑4081936

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