Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,908,501 T→G 41.7% intergenic (‑70/‑176) rstA_2 ← / → acrE_2 Transcriptional regulatory protein RstA/Multidrug export protein AcrE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,908,5010TG41.7% 20.1 / 2.4 37intergenic (‑70/‑176)rstA_2/acrE_2Transcriptional regulatory protein RstA/Multidrug export protein AcrE
Reads supporting (aligned to +/- strand):  ref base T (13/6);  new base G (2/10);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 9.15e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCA  >  KT2440/3908370‑3908629
                                                                                                                                   |                                                                                                                                 
aGGTAGCTGGCGATCAGCTCGGACAGGGCGCGGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTccc                                                                                                                                <  2:380475/135‑1 (MQ=255)
       tGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCCGGTCCCCGTAGGa                                                                                                                        >  8:72145/1‑135 (MQ=255)
          cGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGc                                                                                                                      <  3:248642/134‑1 (MQ=255)
           gATCTGCTCGGACTGGGCGCTGTCGTCTTCGACCAGGAAAATGTTGGGCATAGGTGCTGGGGCGTCTGGGGTGAAGTAAGCCTGTGGCGAGCGCTTGGCGGGTAAACCCGGTGCTGCAGGTCCCCGTAGGAGCgg                                                                                                                    <  8:56440/135‑1 (MQ=255)
               aGCTCGGACAGGGCGCTGTCGTCTTAGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCCGGTAAACCCCCGGGTGCAAGGGCCCGGATCAGCGGGGTc                                                                                                                >  7:59193/1‑135 (MQ=255)
                      aCAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTTAACCCGCTGCTGCAGGTCCCCGGAGGAGCGGGTTCGCCcgc                                                                                                          >  3:121917/1‑134 (MQ=255)
                              cTGTCGTCGTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCCGGTAAACCCGCGGCAGCAGGTCCCCGGAGGCGCCGGTTCCCCCCCGAAGCGGcc                                                                                                 >  4:70851/1‑135 (MQ=255)
                                       tCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGGAGGAGCGGGTGCACCCGCGAAGAGGCCAGCACACGc                                                                                        >  8:322929/1‑135 (MQ=255)
                                         gACCAGCAAAAGGTGCGGCATCGGGGCTGGGGCGGCGGGGTTCAAGTACGCCGGGGGCTAGCGCGGCGCGGTAAACCCCGCGGCGGCCGGCCCCCGTCGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGccc                                                                                      <  8:44649/135‑1 (MQ=12)
                                         ccccAGTAAAATGTGCGGCTTCGGGGGTGGGGCGTGGGGGTTCAAGTAAGCCGGGGGCGAGCGCTGCGCGGGAAACCCCGCGGCTGCAGGCCCCCGTAGGCGCGGGTTCACCCGCGAAGAGGCCAGCACAGGccc                                                                                      <  2:97469/133‑1 (MQ=19)
                                         caCCAGTAAAATGTTCGGCATAGGTGCTGGGGCGGCTGGGGTGAAGTACGCCTGGGGCGAGCGCTGCGGGGTAAAACCCGCTGCTGCCGGGCCCCGTCGGCGCGGGTTCACCCGCGACGAGGCCAGCACAgg                                                                                         <  2:220829/131‑1 (MQ=255)
                                           ccAGGAAAATGTGCGGCATCGGGGCTGGGGCGGCTGGGTTGAAGTAAGCCTGTGGGGAGCGCTGCGCGGGAAACCCCGCTGCTGCAGGGCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCAt                                                                                    <  8:370404/135‑1 (MQ=37)
                                           ccAGCAAAAGGTTCGGCAGTGGGGCTGGGGCGTCTGGGTTGACGTACGCCGGTGGCGAGCGCGGCGCGGTTAACCCCGCGGCTGCAGGCCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCAt                                                                                    <  8:192056/135‑1 (MQ=21)
                                                            aTAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCCGGTCCCCGTAGGAGCGGGTTCCCCCGCGAAGAGGCCAGCACAGGCCCCTGCAGAATATAGAACCCa                                                                   >  5:143375/1‑135 (MQ=255)
                                                                        gCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGGTCACCCGCGCAGACGCCAGCACAGGCCCATGCCGAAAATATAACCCCACCCCTCGACAg                                                       >  3:70478/1‑135 (MQ=37)
                                                                        gCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCCGGTAAACCCGCTGCTGCACGTACCCGGAGGAGCGGGGTCACCCGCCAAGAGGCCACCACAGGCCCCAGCAGCATATAGCACCCAACCAcacccaa                                                        >  6:270892/1‑128 (MQ=255)
                                                                              gggTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCCCCCGCGCAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCCCGCGACAGGATCa                                                  >  1:288223/1‑134 (MQ=255)
                                                                              gggTTGAAGTAAGCCTGTGGCGAGCGCTAAGCGGGTAAACCCGCTGCTGCAGGTCCCCGTATGAGCGGGTTCACCCCCGAAGAGGCCAGCACAGGCCCATCCAGAATATAGATCCCAACCACTCGACACGAACa                                                  >  1:13217/1‑134 (MQ=38)
                                                                                gTTTAATTACGCCTGTGGCTAGCGCTTCGCGGGTAAACCCGCTGCTGCCGGTCCCCGTAGGAGCGGTTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAtt                                                <  3:184095/134‑1 (MQ=38)
                                                                                          gCCTGTGGCGAGCGCTTCGCGGGTTAACCCGCTGCTGCAGGGCCCCGTTGGCGCGGGTTCACCCGCGAAGAGGCCCGCACCGGCCCATGGAGAATATAGAACCCCACCCCTCGCCAGGATCATTTTTgacccgtt                                     >  4:11408/1‑130 (MQ=25)
                                                                                               tgGCGAGCGCTTCGCGGGTAAACCCGCTGCTCCAGGGCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  7:313932/135‑1 (MQ=255)
                                                                                               tgGCGAGCGCTGCGCGGTTACCCGCGCTGCGCCAGCCCCCGTGAGGCGCTGGTACCCCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  4:3471/135‑1 (MQ=255)
                                                                                               gggCTAGCGCTGCGCGGTAAACCCCGCTGCGGCTGGCCCCGGTCGGCGCGGTTACACCGGCTAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  4:325363/134‑1 (MQ=255)
                                                                                               gggCGAGGGCTGCGCGGGAAAACCCGGTGCTCCCGGGCCCCGTAGGAGCGGGTAAACCCGCTAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  2:248302/134‑1 (MQ=255)
                                                                                               gggCGAGCGCTGCGCGGGTAAACCCGCTGCTGCAGGGCCCCGACGGCGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  4:177740/134‑1 (MQ=39)
                                                                                                ggCGAGCGCTTCGGGGGTAAACCCGCTGCTGCAGGGCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  4:121917/134‑1 (MQ=255)
                                                                                                ggCGAGCGCTTCGCGGGTAACCCCGCTGCTGCAGGCCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  7:388672/134‑1 (MQ=255)
                                                                                                       gcTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGGAGGCGCGGGTTCACCCGCGAAGAGGGCCGCACAGGCCCAAGCAGAATATAGCACCCAACCCCTAGACAGGAACATTATTAACACTTTTTGACACTCGcc                         >  8:164923/1‑134 (MQ=37)
                                                                                                       gcTGCGCGGGTAAACCCGCTGCTGCAGGGCCCCGTAGGAGGGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGcc                         <  7:293766/134‑1 (MQ=255)
                                                                                                                      ccGCTGCTGCGGGGCCCCGTACGCGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAg         <  7:120835/135‑1 (MQ=255)
                                                                                                                        ggtgctgcAGGGCCCCGTAGTTGGGGTTTAACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAgcc       <  3:196126/133‑1 (MQ=255)
                                                                                                                            cggcAGGGCCCCGTCGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGACCTACAGGTATTAACAGCCGccc   <  2:26625/133‑1 (MQ=255)
                                                                                                                             tgcAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                >  2:329440/1‑105 (MQ=255)
                                                                                                                             tgcAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  1:329440/105‑1 (MQ=255)
                                                                                                                             tgcAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCa  >  7:79709/1‑135 (MQ=255)
                                                                                                                                  gTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCAt                                                                                    >  6:167440/1‑48 (MQ=255)
                                                                                                                                  gTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCAt                                                                                    <  5:167440/48‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                 
AGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCA  >  KT2440/3908370‑3908629

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