Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 2,134,732 T→G 20.7% L471V (TTG→GTG)  comEC → ComE operon protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,134,7320TG20.7% 59.0 / 2.1 47L471V (TTG→GTG) comECComE operon protein 3
Reads supporting (aligned to +/- strand):  ref base T (15/21);  new base G (7/4);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATCACAGCCTGC  >  ACB122/2134599‑2134861
                                                                                                                                      |                                                                                                                                 
gTTCCGTGGGCTGGCGCTGGGGGCGCAGCAGCGCCCGGCATGGGACCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGtt                                                                                                                                  <  3:247806/135‑1 (MQ=255)
gTTCCGTGGGCGGGCGCGGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGGGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGtt                                                                                                                                  <  3:274884/135‑1 (MQ=255)
  tCCGTGGTCTGGCGCTTGTGGCGCAGCCGCGCCCGTCATGGGTCACGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGACCCGTGGTGT‑GCCGCGGCGCGGGGTgg                                                                                                                                >  2:231541/1‑135 (MQ=255)
  tCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGtt                                                                                                                                  <  5:369203/133‑1 (MQ=255)
         gCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGGTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGGGGGTTggggggggg                                                                                                                         >  6:379270/1‑134 (MQ=255)
              cgcTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGGCCCGGGCGGGGGTGGGGGGGGGTATtg                                                                                                                     >  3:215852/1‑135 (MQ=255)
              cgcTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCGGGTATGTCTGGGTGTGCTGTTGGTTTTTCTGCCCCGTGGTGT‑GGCGCGGGGCGGGGTGGGAGGGGGAAtgc                                                                                                                    >  4:77322/1‑135 (MQ=255)
                   gtggCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGctgctg                                                                                                               <  2:308571/135‑1 (MQ=255)
                        gcagcaGCGCCCGGCAGGGGTCCCGCCGGCGCTGCCGTGGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTTGGGG‑GCGGCGGGGCGGGTTTGGAGGGGTAAAGCTGCTGGCGCt                                                                                                          <  5:165170/135‑1 (MQ=255)
                          agcagcGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGGTGGTTTTGCTGCCCCGTGGTGT‑GCCGCGGCGCGGGGTGGGCGGGGTTATGCGGGTGGCGCggt                                                                                                        >  5:160008/1‑135 (MQ=255)
                           gcagcGCCCGGCAGGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGGATGTCTGGGTGTGCTGTTGGGTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCtgtg                                                                                                       <  2:121347/135‑1 (MQ=255)
                                gcCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTGTTGCTGCCCCGGGGTGT‑GCCGGGGCGGGGGGGGGGGGGGGGTATGGGGCTGGGGCGGGGGGTgg                                                                                                  >  3:29322/1‑134 (MQ=255)
                                    ggCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCgg                                                                                             <  7:100158/136‑1 (MQ=255)
                                         gggTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTTATGCTGCTGGCGCTGTGGGTGCCCCGGGAGc                                                                                         >  5:68859/1‑135 (MQ=255)
                                               cgccgGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGGTGGGGGGGGGTATGCTGCTGGGGGTGGGGGTGCCCCGGGGGCCGGTgg                                                                                   >  5:116822/1‑134 (MQ=255)
                                                   ggCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTGTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GGCGCGGGGGGGGGTGGGGGGGGTTATGGTGGTGGGGCTGGGGGTGGCCCGGGGGGCGGTGGCggt                                                                               >  1:362527/1‑133 (MQ=255)
                                                        tGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCGGCGCGGGTTGGGAGGGGTTATGGTGCGGGCGCTGGGGGTGCCCCGGGTGCCGGTGGCGGTCGGcc                                                                          >  5:203706/1‑135 (MQ=255)
                                                                  ggCCTGGCTGCTTGTATGTCTGGGTGTGGTGTTGGTTGTGCTGGGCCGTGGTGT‑GCCGCTGCGCGGGTTGGTAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGTGAGCCGGTGCCGTTCGGCCAGGTCGAGGTc                                                               <  7:108585/136‑1 (MQ=255)
                                                                  ggCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGt                                                                <  1:134156/135‑1 (MQ=255)
                                                                   gCCTGGCTTCTGGTATGTCTGGGTGTGCTGTTGGGTTTGCTGCCCCGGGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTc                                                               <  5:133239/135‑1 (MQ=255)
                                                                        gctgcgGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCa                                                          <  4:345114/135‑1 (MQ=255)
                                                                             tgGTATGTCTGGGTGTGCTGTTGGGTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGGTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTgg                                                     <  3:98458/135‑1 (MQ=255)
                                                                             tgGTATGGCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTTGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTgg                                                     <  4:192570/135‑1 (MQ=255)
                                                                                    tCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGGGCGGGTTGGGAGGGGTAATGCTGCTGGCGCGGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGCGGTCTGGCAGCTGGATGTTgg                                              >  1:75297/1‑135 (MQ=255)
                                                                                     cTGGGGGTGCTGTTGGTTTGGCTGCCCGGTGGTGG‑GCCGCTGCGCGGGGTGGGAGGGGTAATGCTGCTGGCGCTGGGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTgg                                              <  2:353156/134‑1 (MQ=255)
                                                                                       gggTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGCGCCGGTGCCGTTAGGCCAGGTCGAGGTCTGGCAGCTGGTTGTTGGCCa                                           >  3:385204/1‑135 (MQ=255)
                                                                                        ggTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAg                                          <  3:83210/135‑1 (MQ=255)
                                                                                        ggTGTGCTGTTGGGTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAg                                          <  5:379270/135‑1 (MQ=255)
                                                                                         gtgtGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCTCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAggg                                        >  7:77682/1‑136 (MQ=255)
                                                                                              cTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTg                                    <  6:86573/135‑1 (MQ=255)
                                                                                               tGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTgg                                   >  5:394923/1‑135 (MQ=255)
                                                                                               tGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCGGCGCGGGGTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTgg                                   >  4:148680/1‑135 (MQ=255)
                                                                                                 ttGGTTTTGCTGCCCCGTGGTGT‑GCCGCGGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCGTTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGCGGTCTGGCAGCTGGGTGTTGGGCAGGGGCTGGCgg                                >  8:340621/1‑136 (MQ=255)
                                                                                                      tttGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCAGGATGTTGGCCAGGGGCTGGCGGTGct                            >  6:298514/1‑135 (MQ=255)
                                                                                                              cccGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCTCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTAc                   <  8:77682/136‑1 (MQ=255)
                                                                                                              cccGTGGTGG‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGAGCCCCGGGAGCCGGTGCCGTTCGGCCAGGGCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTAc                   <  8:348563/136‑1 (MQ=255)
                                                                                                              cccGTGGTGG‑GCCGCTGCGCGGGGGGGGAGGGGTAATGCTGCTGGCGCTGTAGGTGCCCCGGGAGCCGGGGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGt                     <  6:372086/134‑1 (MQ=255)
                                                                                                               cggtggggG‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGAGGTTGGCCAGGGGCTTGAGGTGCTCTTGCGTAc                   <  2:249226/133‑1 (MQ=255)
                                                                                                               ccGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTAc                   <  6:394923/135‑1 (MQ=255)
                                                                                                               ccGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTAc                   <  1:258406/135‑1 (MQ=255)
                                                                                                                   ggggt‑gCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTAGGCCAGGGGCTTGCGGTGCTCTTGCGTACGCgg               <  5:52504/132‑1 (MQ=255)
                                                                                                                       t‑gCCGCTGCGCGGGTTGGGACGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGGTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATca          >  7:336949/1‑136 (MQ=255)
                                                                                                                          ccGCTGCGCGGGTTGGGAGGGGTAATGCTGCGGGCGCTGGGGGTGCCCCGGGAGCCGGTGCCGGTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATcac         <  1:157350/135‑1 (MQ=255)
                                                                                                                          ccGCTGCGCGGGTGGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATcac         <  4:66251/135‑1 (MQ=255)
                                                                                                                           cGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCCGCTGGATGTTGGCCAGGGGCTGGCGGTGGTCTTGCGTACCCGGCCTcaca        >  3:170769/1‑135 (MQ=255)
                                                                                                                           cGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTTGCCAGGGGCTGGCGGTGCTAGTGCGTATGCGGCATcaca        >  6:258379/1‑135 (MQ=255)
                                                                                                                           cGCTGCGCGGGGTGGGAGGGGTAATGCGGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATcaca        <  2:172812/135‑1 (MQ=255)
                                                                                                                             cTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGGGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATCACAGc      <  3:77322/135‑1 (MQ=255)
                                                                                                                                cgcgGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATTTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATCACAGCCTg   >  6:132621/1‑135 (MQ=255)
                                                                                                                                cgcgGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATCACAGCCTGc  >  7:287100/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                 
GTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGT‑GCCGCTGCGCGGGTTGGGAGGGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATCACAGCCTGC  >  ACB122/2134599‑2134861

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