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

AGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCA  >  KT2440/3908370‑3908629
                                                                                                                                   |                                                                                                                                 
aGGTAGCTGGCGATCAGCTCGGACAGGGCGCGGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTccc                                                                                                                                <  2:380475/135‑1 (MQ=255)
       tGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCCGGTCCCCGT‑AGGa                                                                                                                        >  8:72145/1‑135 (MQ=255)
          cGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGT‑AGGAGc                                                                                                                      <  3:248642/134‑1 (MQ=255)
           gATCTGCTCGGACTGGGCGCTGTCGTCTTCGACCAGGAAAATGTTGGGCATAGGTGCTGGGGCGTCTGGGGTGAAGTAAGCCTGTGGCGAGCGCTTGGCGGGTAAACCCGGTGCTGCAGGTCCCCGT‑AGGAGCgg                                                                                                                    <  8:56440/135‑1 (MQ=255)
               aGCTCGGACAGGGCGCTGTCGTCTTAGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCCGGTAAACCCCCGGGTGCAAGGGCCCGG‑ATCAGCGGGGTc                                                                                                                >  7:59193/1‑135 (MQ=255)
                      aCAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTTAACCCGCTGCTGCAGGTCCCCGG‑AGGAGCGGGTTCGCCcgc                                                                                                          >  3:121917/1‑134 (MQ=255)
                              cTGTCGTCGTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCCGGTAAACCCGCGGCAGCAGGTCCCCGG‑AGGCGCCGGTTCCCCCCCGAAGCGGcc                                                                                                 >  4:70851/1‑135 (MQ=255)
                                       tCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGG‑AGGAGCGGGTGCACCCGCGAAGAGGCCAGCACACGc                                                                                        >  8:322929/1‑135 (MQ=255)
                                         gACCAGCAAAAGGTGCGGCATCGGGGCTGGGGCGGCGGGGTTCAAGTACGCCGGGGGCTAGCGCGGCGCGGTAAACCCCGCGGCGGCCGGCCCCCGT‑CGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGccc                                                                                      <  8:44649/135‑1 (MQ=12)
                                         ccccAGTAAAATGTGCGGCTTCGGGGGTGGGGCGTGGGGGTTCAAGTAAGCCGGGGGCGAGCGCTGCGCGGGAAACCCCGCGGCTGCAGGCCCCCGT‑AGGCGCGGGTTCACCCGCGAAGAGGCCAGCACAGGccc                                                                                      <  2:97469/133‑1 (MQ=19)
                                         caCCAGTAAAATGTTCGGCATAGGTGCTGGGGCGGCTGGGGTGAAGTACGCCTGGGGCGAGCGCTGCGGGGTAAAACCCGCTGCTGCCGGGCCCCGT‑CGGCGCGGGTTCACCCGCGACGAGGCCAGCACAgg                                                                                         <  2:220829/131‑1 (MQ=255)
                                           ccAGGAAAATGTGCGGCATCGGGGCTGGGGCGGCTGGGTTGAAGTAAGCCTGTGGGGAGCGCTGCGCGGGAAACCCCGCTGCTGCAGGGCCCCGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCAt                                                                                    <  8:370404/135‑1 (MQ=37)
                                           ccAGCAAAAGGTTCGGCAGTGGGGCTGGGGCGTCTGGGTTGACGTACGCCGGTGGCGAGCGCGGCGCGGTTAACCCCGCGGCTGCAGGCCCCCGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCAt                                                                                    <  8:192056/135‑1 (MQ=21)
                                                            aTAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCCGGTCCCCGT‑AGGAGCGGGTTCCCCCGCGAAGAGGCCAGCACAGGCCCCTGCAGAATATAGAACCCa                                                                   >  5:143375/1‑135 (MQ=255)
                                                                        gCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGT‑AGGAGCGGGGTCACCCGCGCAGACGCCAGCACAGGCCCATGCCGAAAATATAACCCCACCCCTCGACAg                                                       >  3:70478/1‑135 (MQ=37)
                                                                        gCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCCGGTAAACCCGCTGCTGCACGTACCCGG‑AGGAGCGGGGTCACCCGCCAAGAGGCCACCACAGGCCCCAGCAGCATATAGCACCCAACCAcacccaa                                                        >  6:270892/1‑128 (MQ=255)
                                                                              gggTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGT‑AGGAGCGGGTTCCCCCGCGCAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCCCGCGACAGGATCa                                                  >  1:288223/1‑134 (MQ=255)
                                                                              gggTTGAAGTAAGCCTGTGGCGAGCGCTAAGCGGGTAAACCCGCTGCTGCAGGTCCCCGT‑ATGAGCGGGTTCACCCCCGAAGAGGCCAGCACAGGCCCATCCAGAATATAGATCCCAACCACTCGACACGAACa                                                  >  1:13217/1‑134 (MQ=38)
                                                                                gTTTAATTACGCCTGTGGCTAGCGCTTCGCGGGTAAACCCGCTGCTGCCGGTCCCCGT‑AGGAGCGGTTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCAtt                                                <  3:184095/134‑1 (MQ=38)
                                                                                          gCCTGTGGCGAGCGCTTCGCGGGTTAACCCGCTGCTGCAGGGCCCCGT‑TGGCGCGGGTTCACCCGCGAAGAGGCCCGCACCGGCCCATGGAGAATATAGAACCCCACCCCTCGCCAGGATCATTTTTgacccgtt                                     >  4:11408/1‑130 (MQ=25)
                                                                                               tgGCGAGCGCTTCGCGGGTAAACCCGCTGCTCCAGGGCCCCGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  7:313932/135‑1 (MQ=255)
                                                                                               tgGCGAGCGCTGCGCGGTTACCCGCGCTGC‑GCCAGCCCCCGTGAGGCGCTGGTACCCCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  4:3471/135‑1 (MQ=255)
                                                                                               gggCTAGCGCTGCGCGGTAAACCCCGCTGCGGCTGGCCCCGGT‑CGGCGCGGTTACACCGGCTAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  4:325363/134‑1 (MQ=255)
                                                                                               gggCGAGGGCTGCGCGGGAAAACCCGGTGCTCCCGGGCCCCGT‑AGGAGCGGGTAAACCCGCTAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  2:248302/134‑1 (MQ=255)
                                                                                               gggCGAGCGCTGCGCGGGTAAACCCGCTGCTGCAGGGCCCCGA‑CGGCGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  4:177740/134‑1 (MQ=39)
                                                                                                ggCGAGCGCTTCGGGGGTAAACCCGCTGCTGCAGGGCCCCGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  4:121917/134‑1 (MQ=255)
                                                                                                ggCGAGCGCTTCGCGGGTAACCCCGCTGCTGCAGGCCCCCGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  7:388672/134‑1 (MQ=255)
                                                                                                       gcTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGG‑AGGCGCGGGTTCACCCGCGAAGAGGGCCGCACAGGCCCAAGCAGAATATAGCACCCAACCCCTAGACAGGAACATTATTAACACTTTTTGACACTCGcc                         >  8:164923/1‑134 (MQ=37)
                                                                                                       gcTGCGCGGGTAAACCCGCTGCTGCAGGGCCCCGT‑AGGAGGGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGcc                         <  7:293766/134‑1 (MQ=255)
                                                                                                                      ccGCTGCTGCGGGGCCCCGT‑ACGCGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAg         <  7:120835/135‑1 (MQ=255)
                                                                                                                        ggtgctgcAGGGCCCCGT‑AGTTGGGGTTTAACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAgcc       <  3:196126/133‑1 (MQ=255)
                                                                                                                            cggcAGGGCCCCGT‑CGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGACCTACAGGTATTAACAGCCGccc   <  2:26625/133‑1 (MQ=255)
                                                                                                                             tgcAGGTCCCCGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                >  2:329440/1‑105 (MQ=255)
                                                                                                                             tgcAGGTCCCCGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTcac                                <  1:329440/105‑1 (MQ=255)
                                                                                                                             tgcAGGTCCCCGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCa  >  7:79709/1‑135 (MQ=255)
                                                                                                                                  gTCCCCGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCAt                                                                                    >  6:167440/1‑48 (MQ=255)
                                                                                                                                  gTCCCCGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCAt                                                                                    <  5:167440/48‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                 
AGGTAGCTGGCGATCAGCTCGGACAGGGCGCTGTCGTCTTCGACCAGGAAAATGTTCGGCATAGGTGCTTGGGCGTCTGGGTTGAAGTAAGCCTGTGGCGAGCGCTTCGCGGGTAAACCCGCTGCTGCAGGTCCCCGT‑AGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCCCATGCAGAATATAGAACCCAACCACTCGACAGGATCATTCTTCACACTTTTTCACACTCGCCCTACAGGTATTAACAGCCGCCCA  >  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: ‑