Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,377,206 A→T 100% L120Q (CTG→CAG)  dnaK ← chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,377,2060AT94.7% 146.9 / ‑3.0 47L120Q (CTG→CAG) dnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (21/23);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATGCGTTTTACGTCCAGGCCAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCAGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGACATCTTCTTCGAAGCGGCGGCCAAT  >  NC_002947/5377061‑5377348
                                                                                                                                                 |                                                                                                                                                
gatgCGTTTTACGTCCAGGCCAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGa                                                                                                                                               >  2:359477/1‑149 (MQ=255)
  tgCGTTTTACGTCCAGGCCAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACt                                                                                                                                             >  4:47610/1‑149 (MQ=255)
   gCGTTTTACGTCCAGGCCAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCa                                                                                                                                                  >  6:180540/1‑143 (MQ=255)
   gCGTTTTACGTCCAGGCCAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCa                                                                                                                                                  <  5:180540/143‑1 (MQ=255)
    cGTTTTACGTCCAGGCCAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTc                                                                                                                                           <  6:242707/149‑1 (MQ=255)
      ttttACGTCCAGGCCAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCg                                                                                                                                          <  5:335425/148‑1 (MQ=255)
          aCGTCCAGGCCAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCAc                                                                                                                                      <  3:257063/148‑1 (MQ=255)
          aCGTCCAGGCCAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACt                                                                                                                                     >  4:224741/1‑149 (MQ=255)
                aGGCCAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTg                                                                                                                               >  4:224776/1‑149 (MQ=255)
                 ggCCAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTgc                                                                                                                              <  4:90256/149‑1 (MQ=255)
                    cAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTgcggc                                                                                                                           >  5:326726/1‑149 (MQ=255)
                    cAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTgcggc                                                                                                                           <  3:47610/149‑1 (MQ=255)
                     aGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCgg                                                                                                                                         <  2:161734/134‑1 (MQ=255)
                     aGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCgg                                                                                                                                         >  1:161734/1‑134 (MQ=255)
                        gATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCATCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGa                                                                                                                        <  5:269948/148‑1 (MQ=255)
                            cggccggcGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCa                                                                                                                    <  6:192122/148‑1 (MQ=255)
                                              gCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTc                                                                                                 >  1:325633/1‑149 (MQ=255)
                                                 tGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGAc                                                                                              >  2:238185/1‑149 (MQ=255)
                                                  gACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGAcc                                                                                             >  2:157777/1‑149 (MQ=255)
                                                   aCGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGAccc                                                                                            <  6:326726/149‑1 (MQ=255)
                                                     gctggctgTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAg                                                                                          <  7:47596/149‑1 (MQ=255)
                                                         gctgTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGc                                                                                        <  3:224776/147‑1 (MQ=255)
                                                               gTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCgt                                                                              >  7:284188/1‑149 (MQ=255)
                                                               gTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCgt                                                                              >  2:129906/1‑149 (MQ=255)
                                                                   aaGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCgttgtt                                                                          >  4:12647/1‑149 (MQ=255)
                                                                           cGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGa                                                                  <  1:238185/149‑1 (MQ=255)
                                                                                       gaCCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTAc                                                       <  1:157777/148‑1 (MQ=255)
                                                                                         ccGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGc                                                    <  1:359477/149‑1 (MQ=255)
                                                                                                 gTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTATTCTTCATTTTCTTCAGGACTTCGGCACTGATCTGCGGCGGCGCCATTTCCTTGCCAGCAGCCTCGCCCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTCTACGGCACAAGCtt                                            >  6:73409/1‑149 (MQ=255)
                                                                                                        gCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑CCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGt                                       <  7:27347/147‑1 (MQ=255)
                                                                                                          tCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCttt                                   >  8:207067/1‑149 (MQ=255)
                                                                                                            gCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCttt                                   <  1:129906/147‑1 (MQ=255)
                                                                                                              cGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGc                               <  3:224741/149‑1 (MQ=255)
                                                                                                                  gTAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGa                           >  6:95456/1‑149 (MQ=255)
                                                                                                                   tAGTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGAc                          <  2:325633/149‑1 (MQ=255)
                                                                                                                     gTCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGACa                         >  7:2032/1‑148 (MQ=255)
                                                                                                                      tCTTCGGCGGTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGACAtc                       <  7:207067/149‑1 (MQ=255)
                                                                                                                        ttgggcggTCTTCTTCATTTTCTTCAGTACTTCGGCACTGTTCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGTCGTCTCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGACAtc                       <  1:106570/146‑1 (MQ=255)
                                                                                                                            gcggTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGACATCttcttc                 <  5:95456/149‑1 (MQ=255)
                                                                                                                             cggTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCGATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGAt                                         <  8:276210/124‑1 (MQ=255)
                                                                                                                             cggTCTTCTTCATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGAt                                         >  7:276210/1‑124 (MQ=255)
                                                                                                                                  ttcttcATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACgg                                                     >  5:294412/1‑107 (MQ=255)
                                                                                                                                  ttcttcATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACgg                                                     <  6:294412/107‑1 (MQ=255)
                                                                                                                                  ttcttcATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATATCCTTGCCAGCAGCCTCGACCCACGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCAACAGCTTGAACTCTTTATGAAAGACATCATCATCGAAgca           >  8:292651/1‑148 (MQ=255)
                                                                                                                                  ttcttcATTTTCTTCAGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGACATCTTCTTCGAAgcg           >  2:188615/1‑149 (MQ=255)
                                                                                                                                    cttcATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGACATCTTCTTCGAAgcggc         >  8:9943/1‑149 (MQ=255)
                                                                                                                                      tcATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGACATCTTCTTCGaa              <  1:194919/142‑1 (MQ=255)
                                                                                                                                      tcATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGACATCTTCTTCGaa              >  2:194919/1‑142 (MQ=255)
                                                                                                                                      tcATTTTCTTCTGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATATTGTACGGCACCAGCTTGATGTCTTTCTGCACGACATCTTCTTCGAAGcggcgg       >  1:198038/1‑149 (MQ=255)
                                                                                                                                            tcttctGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGACATCTTCTTCGAAGCGGCGGCCAAt  <  8:284188/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GATGCGTTTTACGTCCAGGCCAGCGATGCGGCCGGCGTCTTTGGTCGCCTGACGCTGGCTGTCGTTGAAGTAGGCCGGAACGGTGATGACCGCTTCGGTGACTGGCTCGCCGAGGTAGTCTTCGGCGGTCTTCTTCATTTTCTTCAGGACTTCGGCACTGATCTGCGGCGGAGCCATTTCCTTGCCAGCAGCCTCGACCCAGGCG‑‑TCACCGTTGTTCGCCTTGACGATCTTGTACGGCACCAGCTTGATGTCTTTCTGCACGACATCTTCTTCGAAGCGGCGGCCAAT  >  NC_002947/5377061‑5377348

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: