Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,344,114 G→C 15.8% A10G (GCC→GGC) ‡ PP_1168 ← TRAP dicarboxylate transporter subunit DctQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,344,1140GC15.8% 98.6 / 11.6 38A10G (GCC→GGC) ‡PP_1168TRAP dicarboxylate transporter subunit DctQ
Reads supporting (aligned to +/- strand):  ref base G (16/16);  new base C (3/3);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

AGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCG  >  NC_002947/1343971‑1344263
                                                                                                                                               |                                                                                                                                                      
aGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGACGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGCCcgccgc                                                                                                                                                  >  2:222399/1‑150 (MQ=255)
   gggCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACt                                                                                                                                               >  2:571196/1‑150 (MQ=255)
    ggCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGACGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGCCCGCCGCCCTg                                                                                                                                              <  1:222399/150‑1 (MQ=255)
        aaGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCTGCCgc                                                                                                                                          <  1:898815/150‑2 (MQ=255)
         aGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATCGTCATCGCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAg                                                                                                                                         <  2:347913/150‑1 (MQ=255)
               cccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAt                                                                                                                                   <  2:1045997/150‑1 (MQ=255)
                ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAtt                                                                                                                                  >  2:1186213/1‑150 (MQ=255)
                 cACGGGATGATCAGCGTCATGACCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGCAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTg                                                                                                                                 <  2:867218/150‑1 (MQ=255)
                           tCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCTGCCgcttt                                                                                                                                       >  1:1044289/1‑130 (MQ=39)
                           tCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCTGCCgcttt                                                                                                                                       <  2:1044289/134‑5 (MQ=39)
                                gTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGt                                                                                                                  <  1:1129509/150‑1 (MQ=255)
                                  cATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTa                                                                                                                >  1:1067494/1‑150 (MQ=255)
                                           gATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGc                                                                                                      >  1:805249/1‑150 (MQ=255)
                                              cGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGGGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTc                                                                                                   <  2:1001625/150‑1 (MQ=255)
                                                    gCCCGCGATACAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACg                                                                                             >  2:845489/1‑150 (MQ=255)
                                                        gcgATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTtc                                                                                         <  2:805249/150‑1 (MQ=255)
                                                          gATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTtcat                                                                                       <  1:1186213/150‑1 (MQ=255)
                                                               aGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGt                                                                                  <  2:256148/150‑1 (MQ=255)
                                                                         gCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTg                                                                        >  2:474146/1‑150 (MQ=255)
                                                                            gTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgcgc                                                                     >  1:888491/1‑150 (MQ=255)
                                                                                     gtCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTAttt                                                            <  1:474146/150‑1 (MQ=255)
                                                                                              gCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCCGGCCGGCGGGCCCCGCACTGGCGAGGTCGATTGGGAAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTAc                                                   <  1:571196/150‑1 (MQ=255)
                                                                                              gCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTAc                                                   >  1:731514/1‑150 (MQ=255)
                                                                                                    gcaaaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGc                                             >  1:705555/3‑150 (MQ=255)
                                                                                                     aaaaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGcc                                            <  1:845489/150‑1 (MQ=255)
                                                                                                      aaaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCg                                           >  2:471521/1‑150 (MQ=255)
                                                                                                      aaaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGC‑GCGGCCAGGTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTGGCGTACCGAAGCCg                                           >  1:811927/1‑150 (MQ=255)
                                                                                                             tttCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgtc                                    <  1:68461/150‑1 (MQ=255)
                                                                                                              ttCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgtcg                                   <  2:1067494/150‑1 (MQ=255)
                                                                                                                cATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTC‑CTGGACGTTGTCGATACGATTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTa                                 >  2:286494/1‑150 (MQ=255)
                                                                                                                     gCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcgc                                                       <  2:1155211/123‑1 (MQ=255)
                                                                                                                     gCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcgc                                                       >  1:1155211/1‑123 (MQ=255)
                                                                                                                          aCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑Ct                                                                                                             >  2:784517/1‑64 (MQ=255)
                                                                                                                          aCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑Ct                                                                                                             <  1:784517/64‑1 (MQ=255)
                                                                                                                           ccAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttct                      >  1:307581/1‑150 (MQ=255)
                                                                                                                               ctcGGCAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTc                  >  1:865430/3‑150 (MQ=255)
                                                                                                                                     aGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGc                                       <  2:377757/123‑1 (MQ=255)
                                                                                                                                     aGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGc                                       >  1:377757/1‑123 (MQ=255)
                                                                                                                                             gggcgCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACt    >  2:374240/1‑150 (MQ=255)
                                                                                                                                               gcgcCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCg  >  1:896823/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
AGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCG  >  NC_002947/1343971‑1344263

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