Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,821,793 T→A 8.7% T30S (ACC→TCC)  ptxD ← phosphonate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,821,7930TA 8.7% 148.9 / 2.6 54T30S (ACC→TCC) ptxDphosphonate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (19/30);  new base A (3/2);  total (22/32)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCACGCCTACCGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCAGAG  >  NC_002947/3821647‑3821939
                                                                                                                                                   |                                                                                                                                                   
gTCCACGCCTACCGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGgtggtg                                                                                                                                                    <  1:206482/148‑1 (MQ=255)
   cACGCCTACCGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCa                                                                                                                                                <  6:276835/149‑1 (MQ=255)
        cTACCGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAg                                                                                                                                          <  1:262900/149‑1 (MQ=255)
        cTACCGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAg                                                                                                                                          <  3:81353/149‑1 (MQ=255)
           ccGAAACGCTCGACACCTCTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGTATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTc                                                                                                                                       <  1:54851/149‑1 (MQ=255)
           ccGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGt                                                                                                                                        <  5:93317/148‑1 (MQ=255)
           ccGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTc                                                                                                                                       <  2:11271/149‑1 (MQ=255)
           ccGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCAAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGACACCACCCAgc                                                                                                                                         <  2:164025/147‑2 (MQ=255)
            cGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCCACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACACAGGTCa                                                                                                                                      >  6:243519/1‑149 (MQ=255)
            cGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCa                                                                                                                                                >  4:171701/1‑140 (MQ=255)
            cGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCa                                                                                                                                                <  3:171701/140‑1 (MQ=255)
            cGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCCCCCAGGTCa                                                                                                                                      >  5:218176/1‑149 (MQ=255)
              aaaCGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCAc                                                                                                                                     >  4:31199/1‑148 (MQ=255)
                      gACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGAcgc                                                                                                                             <  7:34143/148‑1 (MQ=255)
                        caccacTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGAt                                                                                                                          <  5:328955/149‑1 (MQ=255)
                         accacTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATc                                                                                                                         >  7:129233/1‑149 (MQ=255)
                              ttCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGCGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTg                                                                                                                      >  1:231219/1‑148 (MQ=255)
                              ttCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGACACCACCCAg                                                                                                                                          <  5:58647/127‑1 (MQ=255)
                              ttCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGACACCACCCAg                                                                                                                                          >  6:58647/1‑127 (MQ=255)
                                ccAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACTACCCGTGCGCCCCGGGTAACGCATCGCGCTGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCAGCCAGGTCACCTCGACGCGATCTTGCAg                                                                                                                   <  2:238720/148‑1 (MQ=255)
                                ccAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCa                                                                                                                    <  3:58893/147‑1 (MQ=255)
                                 cAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGt                                                                                                                                        <  7:233510/126‑1 (MQ=255)
                                 cAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGt                                                                                                                                        >  8:233510/1‑126 (MQ=255)
                                  aGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCg                                                                                                                <  1:1670/149‑1 (MQ=255)
                                   gTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGt                                                                                                               <  8:215901/149‑1 (MQ=255)
                                   gTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGt                                                                                                               <  7:193190/149‑1 (MQ=255)
                                                  cGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAgtcgtc                                                                                                 <  1:67608/149‑1 (MQ=255)
                                                     gCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGa                                                                                              >  6:127101/1‑148 (MQ=255)
                                                     gCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGAc                                                                                             <  2:231219/149‑1 (MQ=255)
                                                           tGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGt                                                                                       <  3:31199/149‑1 (MQ=255)
                                                               aTCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTtat                                                                                   <  5:23934/149‑1 (MQ=255)
                                                                tCCATGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTtata                                                                                  <  4:92994/149‑1 (MQ=255)
                                                                  cAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACa                                                                                <  4:40143/149‑1 (MQ=255)
                                                                     gcgcAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGa                                                                             <  5:126081/149‑1 (MQ=255)
                                                                            gCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACCGGACGATTc                                                                       >  6:77737/1‑149 (MQ=255)
                                                                                gCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGc                                                                     <  5:117088/147‑1 (MQ=255)
                                                                                    ccAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGcttctt                                                                >  5:291979/1‑148 (MQ=255)
                                                                                        caacCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTg                                                           >  4:256464/1‑149 (MQ=255)
                                                                                           cccGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAg                                                        <  3:256464/149‑1 (MQ=255)
                                                                                              gTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCa                                                                                                                                                >  2:262213/1‑59 (MQ=255)
                                                                                                     cccGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAAt                                              >  3:268876/1‑149 (MQ=255)
                                                                                                              gCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATgg                                     <  5:127101/149‑1 (MQ=255)
                                                                                                              gCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATgg                                     >  3:222114/1‑149 (MQ=255)
                                                                                                              gCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCACCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATgg                                     >  7:238860/1‑149 (MQ=255)
                                                                                                                   gcgcAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGa                                >  8:34886/1‑149 (MQ=255)
                                                                                                                      cAGGCGGGCCAGGGCATCGGGCTGGCTGGTGTCCCCCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGAcgc                             >  3:131271/1‑149 (MQ=255)
                                                                                                                       aGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGAcgcg                            <  4:268876/149‑1 (MQ=255)
                                                                                                                        ggCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCACCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGa                           <  8:238860/149‑1 (MQ=255)
                                                                                                                         gCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAg                          >  6:200739/1‑149 (MQ=255)
                                                                                                                            ggCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTa                       <  5:200739/149‑1 (MQ=255)
                                                                                                                              ccAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATc                     <  4:222114/149‑1 (MQ=255)
                                                                                                                                  ggCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCgcgc                 >  3:292027/1‑149 (MQ=255)
                                                                                                                                          gCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCg         <  6:291979/149‑1 (MQ=255)
                                                                                                                                           cTGGGTGGACACCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCgg        >  8:4911/1‑149 (MQ=255)
                                                                                                                                            tGGGTGGACACCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGc       >  7:89965/1‑149 (MQ=255)
                                                                                                                                                  gtgtCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCagag  <  2:323536/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTCCACGCCTACCGAAACGCTCGACACCACTTCCAGTTGCGGGGCAAGGTCGAGCAGGCTGGCATCCAGGCGCAGGCTGGCCCCCAACAACCCGTGCGCCCCCGGTAACGCATCGCGCAGGCGGGCCAGGGCATCGGGCTGGGTGGTGTCCACCCAGGTCACCTCGACGCGATCTTGCAGGCGTGCCATCAAGTCGTCGGACAGGCGTTTATACAGGACGATTCGCTTCTTCATTGCAGGTTTCCTAATGGTTCATGGCCAGACGCGAGCTATCGCGCTTTGCCCGGCCAGAG  >  NC_002947/3821647‑3821939

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: