Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 2,231,813 A→C 36.6% G291G (GGT→GGG)  emrA_2 ← Colistin resistance protein EmrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,231,8130AC36.6% 20.2 / 3.5 28G291G (GGT→GGG) emrA_2Colistin resistance protein EmrA
Reads supporting (aligned to +/- strand):  ref base A (4/13);  new base C (9/2);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 5.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGCGAATACGCCCATGCAATGGGTTGTCGGCGTGGAACGTGAGTTTTATCGGGATGCGCTGCACAACTTTGGTGAAGTTGCCGGTGGCGTTGTCCGGTGGCAGCAGGCTGAACTGCGCACCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCCGCGCGCTGGCCGGGGTGCATGCGGCCGATCTGGGTTTCCTTGAAGTTGGCCTGCACC  >  KT2440/2231693‑2231942
                                                                                                                        |                                                                                                                                 
ggggCGAATACGCCCATGCAATGGGTTGTCGGCGTGGAACGTGAGTTTTATCGGGATGCGCTGCACAACTTTGGTGAAGTTGCCGGTGGCGTTGTCCGGTGGCAGCAGGCTGAACTGCGCCCCGGCCgcggcgc                                                                                                                      >  4:129883/1‑133 (MQ=255)
    cGAATACGCCCATGCAATGGGTTGTCGGCGTGGAACGTGAGTTTTATCGGGATGCGCTGCACAACTTTGGTGAAGTTGCCGGTGGCGTTGTGCGGGGGCCGCAGGGGGGACTGCGCCCCCGGCGCGGCGCAcccg                                                                                                                 >  3:263106/1‑132 (MQ=37)
              cATGCAAGGGGTTGACGGCGTGGAACGTGAGTTTTATCGGGTTGCGCGGCCAAAATTTGGTGAAGTTGCGGGTGGCGTTTGCCGGGGGGAGCCGGGGGCAAGGGGCCCCGCACGCGGCGACCACGCGTTcccccc                                                                                                       >  8:125257/1‑135 (MQ=255)
                tGCAATGGGTTGTCGGCGTGGAACGTGAGTTTTATCGGGATGCGCTGCACAACTTTGGTGAAGTTGCCGGTGGCGTTGTACGGGTGCCGCAGGCGGAACAGCGCCCCGGACGCGGCGCACACGCTGTCCCCCCgg                                                                                                     >  2:231868/1‑135 (MQ=37)
                   aaTGGGTTGTCGGCGTGGAACGTGAGTTTTATCGGGATGCGCTGCACAACTTTGGTGAAGTTGCCGGTGGCGTTGTGCGGGGGGCGCCGGGGGGCCAGGGCCCCCGGCGCGGCGCACAGGCGGTGCCCCCCGccg                                                                                                  >  5:113884/1‑134 (MQ=19)
                   aaTGGGTTGTCGGCGTGGAACGTGAGTTTTATCGGGATGCGCTGCACAACTTTGGTGAAGTTGCCGGTGGCGTTGTACGGTGGCCGCCGGCGGGACAGCGCCCCCGCCCCGGCGCACAGGCTGTCCACCCGGccg                                                                                                  >  7:199864/1‑134 (MQ=37)
                        gTTGTCGGCGTGGAACGTGAGTTTTATCGGGATGCGCTGCACAACTTTGGTGAAGTTGCCGGTGGCGTTGTCCCGTGGCAGCAGGCTGAACTGCGCACCGGACGCGGCGCACACGCTGTCCCCCCCGCCGTCgct                                                                                             >  7:35200/1‑133 (MQ=38)
                             cggcggGTAACGTGAGTTTTTTCGGGATGCGCTGAACAACTTTGGTGATGTTGCCGGGGGCGTGGTCGGGTGGCAGCAGGCTGAACTGCGCACGGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGgtg                                                                                        <  7:214955/135‑1 (MQ=255)
                              ggCGTGGAACGTGAGTTTTATCGGGATGCGCTGCACAACTTTGGTGAAGTTGCCGGTGGCGTTGTACCGGGGCAGCCGGCGGCACTGCGCCCCCGGCGCGGGGGACCGGGTGTGCCCCCCGCCGTCCATccgggg                                                                                       >  5:282248/1‑130 (MQ=19)
                                        gggAGTTTTATGGGGAGGGGCTGCACAAGTTTGTGGATGTGGCGGGGGGCGTTGGCGGGGGGCCGCTGGCTGAACGGCGCCCCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAg                                                                             <  2:189584/133‑1 (MQ=255)
                                            gTTTTATCGGGATGCGCTGCAAAACTTTGGTGAAGTTGCGGGTGGCGTTGTCCGGGGGCAGCAGGCTGAACTCCGCACCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTg                                                                          <  2:23245/134‑1 (MQ=255)
                                                    gggATGCGCTGCACAACTTTGGTGAAGTTGCCGGTGGCGTTGTCCGGTGGCAGCAGGCTGAACTGCGCACCGGACGCGGCGAACAGGCTGTCCCCCCGGCCTTCCATCGGTGTGTCGGGGTAGCTGTcccaca                                                                   >  8:236718/1‑128 (MQ=255)
                                                          cgcTGCACAACTTTGGTGAAGTTGCCGGTGGCGTTGTCCCGTGGCAGCAGGCGGAACAGCGCCCCCGCCGCGGCGAACAGGCTGTCCACCCGGCCTTCGCTCGGTGTGTCGGGGTATCTGTCGCACAACACttac                                                           >  6:217820/1‑133 (MQ=255)
                                                            cTGCACACCTTTGGTGAAGGGGCGGGGGGCGTTGTCCGGTGGCAGCAGGCTGAACTGCGCACCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCcgc                                                         <  7:12214/135‑1 (MQ=255)
                                                             tGCACAACTTTGGTGAAGTTGCCGGTGGCGTTGGCCCGTGGCCGCCGGCGGCACAGCGCCCCCGACCCCGCGCACCCGCGGTCCACCCCGCCGGCCCTCCGTGTGTCGGGGTGGCCGGCGAACCCCAGTTCcgcc                                                        >  6:94291/1‑134 (MQ=255)
                                                                acaACTTTGGTGATGGTCCCGGGGGGGTGGGCCGGGGGCAGCAGGCTGAACGGCGCACCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCCgcgcgc                                                     <  8:199864/135‑1 (MQ=255)
                                                                       tGGTGAAGTTGCCGGTGGCGTTGTCCGGTGGCAGCAGGCTGAACTGCGCACCGGCCGCGGCGCACAGGCTGTCCCCCCGGCCGTCGCTCGGGGGGTGGGGGTTGCTGTCGCACAACACTTACCCGCGCGGGCCgg                                              >  3:210612/1‑135 (MQ=255)
                                                                                    ggTGGCGTTGTCCGGTGGCAGCAGGCTGAACTGCGCACCGGCCGCGGCGCACAGGCTGTCCACCCGGCCTTCGATGGGTGTGTCGGGGTTGCTGTCGAACAACCGTTACGCGCGCTGGCCGGGGTGGATGCGGcc                                 >  3:109174/1‑135 (MQ=255)
                                                                                          tttGTCCGGGGGCAGCTGGCTGAACTGCGCACCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCCGCGCGCTGGCCGGGGTGCATGCGGCCGATCTg                           <  3:111990/134‑1 (MQ=255)
                                                                                               ccGGGGGCAGCAGGCTTAACTGCGCACCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCCGCGCGCTGGCCGGGGTGCATGCGGCCGATCTGGGttt                      <  1:78416/135‑1 (MQ=255)
                                                                                                cGGGGGCAGCAGGCTGAACTGCGCACCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCCGCGCGCTGGCCGGGGTGCATGCGGCCGATCTGGGTTTc                     <  5:153580/135‑1 (MQ=255)
                                                                                                         cagGCTGAACTGCGCACCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCCGCGCGCTGGCCGGGGTGCATGCGGCCGATCTGGGTTTCCTTGAAGtt            <  4:263106/135‑1 (MQ=255)
                                                                                                         cagGCTGAACGGCGCACGGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCCGCGCGCTGGCCGGGGTGCATGCGGCCGATCTGGGTTTCCTTGAAGtt            <  1:59251/135‑1 (MQ=255)
                                                                                                           ggCTGAACTGCGCACCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCCGCGCGCTGGCCGGGGTGCATGCGGCCGATCTGGGTTTCCTTGAAGTTgg          <  1:48600/135‑1 (MQ=255)
                                                                                                           ggCTGAACTGCGCACCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCCGCGCGCTGGCCGGGGTGCATGCGGCCGATCTGGGTTTCCTTGAAGTTgg          <  2:242880/135‑1 (MQ=255)
                                                                                                              tGAACGGCGCACGGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCCGCGCGCTGGCCGGGGTGCATGCGGCCGATCTGGGTTTCCTTGAAGTTGGCCt       <  5:94291/135‑1 (MQ=255)
                                                                                                                accTGCGCCCGGGCCGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGGCGGGGTAGCTGTCGAACAACAGTTCCGCGCGCTGGCCGGGGTGCATGCGGCCGATCTGGGTTTCCTTGAAGTTGGCCTGc     <  5:43850/133‑1 (MQ=255)
                                                                                                                   tGCCCACCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCCGCGCGCTGGCCGGGGTGCATGCGGCCGATCTGGGTTTCCTTGAAGTTGGCCTGCAcc  <  7:18988/135‑1 (MQ=255)
                                                                                                                        |                                                                                                                                 
GGGGCGAATACGCCCATGCAATGGGTTGTCGGCGTGGAACGTGAGTTTTATCGGGATGCGCTGCACAACTTTGGTGAAGTTGCCGGTGGCGTTGTCCGGTGGCAGCAGGCTGAACTGCGCACCGGACGCGGCGAACAGGCTGTCCACCCGGCCTTCGATCGGTGTGTCGGGGTAGCTGTCGAACAACAGTTCCGCGCGCTGGCCGGGGTGCATGCGGCCGATCTGGGTTTCCTTGAAGTTGGCCTGCACC  >  KT2440/2231693‑2231942

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