Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,344,106 G→T 10.2% P13T (CCT→ACT)  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,1060GT10.2% 117.1 / 3.1 39P13T (CCT→ACT) PP_1168TRAP dicarboxylate transporter subunit DctQ
Reads supporting (aligned to +/- strand):  ref base G (17/18);  new base T (2/2);  total (19/20)
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.77e-01

CAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTC  >  NC_002947/1343963‑1344247
                                                                                                                                               |                                                                                                                                              
cAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGgcggcg                                                                                                                                          <  1:1096555/150‑1 (MQ=255)
        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)
                                         tCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAggg                                                                                                                                               >  2:953646/1‑104 (MQ=255)
                                         tCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAggg                                                                                                                                               <  1:953646/104‑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                                       >  1:1155211/1‑123 (MQ=255)
                                                                                                                             gCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcgc                                       <  2:1155211/123‑1 (MQ=255)
                                                                                                                                  aCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑Ct                                                                                             <  1:784517/64‑1 (MQ=255)
                                                                                                                                  aCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑Ct                                                                                             >  2:784517/1‑64 (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)
                                                                                                                                               |                                                                                                                                              
CAGCACGTAGCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTA‑CTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTC  >  NC_002947/1343963‑1344247

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: ‑