Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,376,913 C→G 34.1% V218L (GTA→CTA)  dnaK ← chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,376,9130CG34.1% 44.4 / 42.5 44V218L (GTA→CTA) dnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (14/15);  new base G (9/6);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

TCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGCCTTGTCCATGCCGTAAGCCAGTGCAGCGGCAGTCGGCTCGTT  >  NC_002947/5376778‑5377057
                                                                                                                                       |                                                                                                                                                
tCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTgg                                                                                                                                        <  7:110193/146‑1 (MQ=255)
        tGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTc                                                                                                                             >  4:268505/1‑149 (MQ=255)
              gCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGGCGACAAGGTAGTCGATCAGGCGCATGGCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTc                                                                                                                       <  4:213151/149‑1 (MQ=255)
               ccAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCg                                                                                                                      >  7:276355/1‑149 (MQ=255)
                             ttcAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGAc                                                                                                        >  6:48317/1‑149 (MQ=255)
                                aGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCg                                                                                                      >  4:306768/1‑148 (MQ=255)
                                aGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGa                                                                                                     <  8:42207/149‑1 (MQ=255)
                                aGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCg                                                                                                      >  5:10123/1‑148 (MQ=255)
                                         cccGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGaa                                                                                            <  6:264664/149‑1 (MQ=255)
                                            gACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGt                                                                                         <  5:48317/149‑1 (MQ=255)
                                                ctttcttGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTacca                                                                                     >  6:246366/1‑149 (MQ=255)
                                                ctttcttGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTacca                                                                                     >  7:128415/1‑149 (MQ=255)
                                                 tttcttGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTaccac                                                                                    >  2:318551/1‑149 (MQ=255)
                                                 tttcttGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTaccac                                                                                    >  7:117826/1‑149 (MQ=255)
                                                     ttGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTaccaccacc                                                                                >  3:81253/1‑149 (MQ=255)
                                                          cTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGt                                                                           <  4:81253/149‑1 (MQ=255)
                                                           tcgtcgACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGt                                                                           >  6:105787/1‑148 (MQ=255)
                                                             gtcgACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGt                                                                                         <  4:66651/132‑1 (MQ=255)
                                                             gtcgACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGt                                                                                         >  3:66651/1‑132 (MQ=255)
                                                              tcgACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATa                                                                       <  8:128415/149‑1 (MQ=255)
                                                               cgACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATaa                                                                      <  4:102786/149‑1 (MQ=255)
                                                                     ggTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGAt                                                                 >  8:178185/1‑148 (MQ=255)
                                                                         gTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACg                                                             <  1:268781/148‑1 (MQ=255)
                                                                             aTCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGtgtg                                                        <  5:105787/149‑1 (MQ=255)
                                                                              tCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTgg                                                       <  7:178185/149‑1 (MQ=255)
                                                                                  gcgcATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGc                                                   <  1:318551/149‑1 (MQ=255)
                                                                                  gcgcATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGc                                                   >  6:163271/1‑149 (MQ=255)
                                                                                    gcATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGccc                                                 >  3:289129/1‑149 (MQ=255)
                                                                                       tGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTg                                              <  8:276355/149‑1 (MQ=255)
                                                                                        gTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTgg                                             >  3:105903/1‑149 (MQ=255)
                                                                                        gTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTgg                                             >  1:226424/1‑149 (MQ=255)
                                                                                           gAAGTCTTCGCCACCCAGGAAGGTGTCGCCGGTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGcc                                           <  6:10123/148‑1 (MQ=255)
                                                                                               tCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGCCTTGTc                                      <  2:226424/149‑1 (MQ=255)
                                                                                                 ttCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGCCTTGTCCa                                    >  6:194298/1‑149 (MQ=255)
                                                                                                 ttCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGCCTTGTCCa                                    >  3:84542/1‑149 (MQ=255)
                                                                                                    gCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGCCTTGTCCATg                                  <  3:306768/148‑1 (MQ=255)
                                                                                                         ccAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGCCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGCCTTGTCCATGCCGTaa                            <  3:268505/149‑1 (MQ=255)
                                                                                                           aGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGCCTTGTCCATGCCGTAAGc                          <  5:194298/149‑1 (MQ=255)
                                                                                                             gAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGCCTTGTCCATGCCGTAAGc                          >  3:118704/1‑147 (MQ=255)
                                                                                                                ggTGTCGCCGTTGGTAGCCAGTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGCCTTGTCCATGCCGTAAGCCAGTg                     <  4:105903/149‑1 (MQ=255)
                                                                                                                ggTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGCCTTGTCCATGCCGTAAGCCAGTg                     <  8:117826/149‑1 (MQ=255)
                                                                                                                      gCCGTTGGTAGCCAGTAGTTCGAAATGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACAACCAACCAGGTCATAAAAGATGACGGTGTGGTCGACCTTGGCCTTGTCCATTCCGTAAGCCAGTGCAgagt               >  6:26188/1‑146 (MQ=255)
                                                                                                                              tAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGCCTTGTCCATGCCGTAAGCCAGTGCAGCGGCAGTCgg        >  2:355223/1‑148 (MQ=255)
                                                                                                                                    gTAGTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTAGGCCTTGTCCATGCCGTAAGCCAGTGCAGCGGCAGTCGGCTCGtt  >  4:268840/1‑148 (MQ=255)
                                                                                                                                       |                                                                                                                                                
TCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGATGACGGTGTGGTCGCCCTTGGCCTTGTCCATGCCGTAAGCCAGTGCAGCGGCAGTCGGCTCGTT  >  NC_002947/5376778‑5377057

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