Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,002,898 G→A 9.7% L70L (CTG→CTA)  bcsE → Cyclic di‑GMP binding protein BcsE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,002,8980GA 9.7% 120.9 / 3.2 41L70L (CTG→CTA) bcsECyclic di‑GMP binding protein BcsE
Reads supporting (aligned to +/- strand):  ref base G (18/19);  new base A (2/2);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGAGCGCCTGCGCGTTTTCAACCAGTGCCTCGATGGCCTGGGGCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAGCGCCCGGTACTGGAAGGCGCGCCAGCCTTCTCCGAACACTGGCA  >  KT2440/3002764‑3003029
                                                                                                                                      |                                                                                                                                   
ggcGAGCGCCTGCGCGTTTTCAACCAGTGCCTCGATGGCCTGGGGCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTg                                                                                                                                     <  5:78617/133‑1 (MQ=255)
 tCGAGCGCCTGCGCGTTTTCAACCAGTGCCTCGATGGCCTGGGGCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAACACCGCCACCCGGGACATGGAGCTGa                                                                                                                                    >  6:139236/1‑135 (MQ=255)
 tCGAGCGCCTGCGCGTTTTCAACCAGTGCCTCGATGGCCTGGGGCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGa                                                                                                                                    >  2:182393/1‑135 (MQ=255)
      cgcCTGCGCGTTTTCAACCAGTGCCTCGATGGCCTGGGGCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCgg                                                                                                                               >  2:113430/1‑135 (MQ=255)
        ccTGCGCGTTTTCAACCAGTGCCTCGATGGCCTGGGGCAGGTGTACCGCGGCTAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTAACTCGGCt                                                                                                                             >  7:105233/1‑135 (MQ=255)
        ccTGCGCGTTTTCAACCAGTGCCTCGATGGCCTGGGGCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTAACTCGGc                                                                                                                              >  7:109964/1‑134 (MQ=255)
              cgTTTTCAACCAGTGCCTCGATGGCCTGGGGCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAgcgc                                                                                                                       >  4:204747/1‑135 (MQ=255)
                                  aTGGCCTGGGGCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGTCTGAGCTCTGTCTTTGCCGTtgctgc                                                                                                     <  7:209442/133‑1 (MQ=255)
                                  aTGGCCTGGGGCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTtgctgc                                                                                                     <  3:265931/133‑1 (MQ=255)
                                  aTGGCCTGGGGCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTtgctgc                                                                                                     <  3:204747/133‑1 (MQ=255)
                                      ccTGGGGCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCcgcg                                                                                                 <  3:29855/133‑1 (MQ=255)
                                        tGGGGCAGGTGTAACGCGGCAAAGGTGGGGAGCGCTACCTGCAGCATTTCTGGATCAACCCGCTTGGCAACACCTCCACCCGGGACATGGAGCTGACTAGGCTGAGCGCGTGCTTTGCCTTTGCTTCCGCGcc                                                                                               >  8:309151/1‑133 (MQ=255)
                                          gggCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAg                                                                                           <  3:61489/135‑1 (MQ=255)
                                                  gtACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAg                                                                                   >  3:147891/1‑135 (MQ=255)
                                                   tACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAg                                                                                   >  4:136400/1‑134 (MQ=255)
                                                             aaGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACACGCTTGGCAAAACCGCAACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCggtg                                                                        <  1:56848/135‑1 (MQ=255)
                                                                ggggggTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGgtggtg                                                                     <  4:246551/133‑1 (MQ=255)
                                                                     gTGCGCTACCTGCTGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAg                                                                <  1:97521/135‑1 (MQ=255)
                                                                       gcgcTACCTGCAGCATTTGTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGtt                                                              <  2:224585/135‑1 (MQ=255)
                                                                       gcgcTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCTCCTCAGGCACCCAGCGACACCGGTGGTGACGAGtt                                                              <  5:208993/135‑1 (MQ=255)
                                                                       gcgcTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGtt                                                              >  5:62075/1‑135 (MQ=255)
                                                                       gcgcTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGtt                                                              <  1:182393/135‑1 (MQ=255)
                                                                        cgcTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGCCATGGCGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTg                                                             >  7:160368/1‑135 (MQ=255)
                                                                           tACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCtg                                                          >  3:248837/1‑135 (MQ=255)
                                                                             ccTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCCGCGACACCGGTGGTGACGAGTTGCtgtg                                                        >  3:296023/1‑135 (MQ=255)
                                                                                  agcaTTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTg                                                    <  5:111351/134‑1 (MQ=255)
                                                                                         cTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAgc                                             <  3:326039/134‑1 (MQ=255)
                                                                                         cTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAgc                                             <  8:281524/134‑1 (MQ=255)
                                                                                         cTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAgc                                             <  2:102511/134‑1 (MQ=255)
                                                                                         cTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTAACTCGGCTGAGCGCGGTCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAgc                                             <  8:109964/134‑1 (MQ=255)
                                                                                         cTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTAACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAgc                                             <  8:105233/134‑1 (MQ=255)
                                                                                          tGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAgcgc                                           >  3:198781/1‑135 (MQ=255)
                                                                                            gAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAgc                                             >  3:40113/1‑131 (MQ=255)
                                                                                            gAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAgc                                             <  4:40113/131‑1 (MQ=255)
                                                                                                    cGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAgc                                             >  1:161268/1‑123 (MQ=255)
                                                                                                    cGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAgc                                             <  2:161268/123‑1 (MQ=255)
                                                                                                    cGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAGCGCCCGGTCCTg                                  >  7:66169/1‑134 (MQ=255)
                                                                                                    cGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAGCGCCCGGTACTg                                  >  2:228607/1‑134 (MQ=255)
                                                                                                                              aTGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAGCGCCCGGTACTGGAAGGCGCGCCAGCCTTCTCCGAacac       <  2:114900/135‑1 (MQ=255)
                                                                                                                                   gCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCCGCGCCACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAGCGCCCGGTACTGGAAGGCGCGCCAGCCTTCTCCGAACACTGGCa  >  4:103309/1‑135 (MQ=255)
                                                                                                                                   gCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAGCGCCCGGTACTGGAAGGCGCGCCAGCCTTCTCCGAACACTGGCa  <  5:120565/135‑1 (MQ=255)
                                                                                                                                   gCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAGCGCCCGGTACTGGAAGGCGCGCCAGCCTTCTCCGAACACTGGCa  >  4:296235/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                   
GTCGAGCGCCTGCGCGTTTTCAACCAGTGCCTCGATGGCCTGGGGCAGGTGTACCGCGGCAAAGGTGGGGTGCGCTACCTGCAGCATTTCTGGAGCAACCCGCTTGGCAAAACCGCCACCCGGGACATGGAGCTGACTCGGCTGAGCGCGGGCTTTGCCGTTGCTGCCGCGCCTCAGGCACCCAGCGACACCGGTGGTGACGAGTTGCTGTGCCTGGCCCAGCGCCCGGTACTGGAAGGCGCGCCAGCCTTCTCCGAACACTGGCA  >  KT2440/3002764‑3003029

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