Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,491,259 A→C 33.3% A34A (GCA→GCC)  iutA → Ferric aerobactin receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,491,2590AC33.3% 58.5 / 7.8 57A34A (GCA→GCC) iutAFerric aerobactin receptor
Reads supporting (aligned to +/- strand):  ref base A (22/7);  new base C (2/12);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 2.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCAC  >  ACB143/2491125‑2491390
                                                                                                                                       |                                                                                                                                   
aTTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCCCCATGCCTTGATTTCTTGGAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAg                                                                                                                                    >  8:374618/1‑136 (MQ=255)
  tctcATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGGTATTTTAACCATGCCTTGATTTATTGCAGTTCCGTGGCCCTGTTGCTCGGCCCCCGGCAC‑GCAggg                                                                                                                                  >  3:287766/1‑136 (MQ=255)
          aCGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTGTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCAGCAC‑GCAGGGGCcgcgcc                                                                                                                          >  5:76814/1‑134 (MQ=255)
           cGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTGACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCGGCCC‑CCAGGGGCCGCGACg                                                                                                                         >  5:307621/1‑136 (MQ=255)
           cGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCAATCCTTGATGTCTTGCAGTGCCCTGGCCCAGTTGCTCGGCCCCCAGCAC‑GCAGGGGCCGCGCCg                                                                                                                         >  5:147868/1‑136 (MQ=255)
             aCTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGGTCTTTTTACCAAGCCATGATGTGTTGCAGTGCCGTGGGCCAGTTTCTCGGCCCCCAGCAC‑GCCCGGGCCGCGACCGc                                                                                                                       >  1:353334/1‑136 (MQ=255)
                   cAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGGTGTTTTCACCAAGCCGTGGTGTGTTTCCGTTCCCTGGCCCAGGTTGTCGGCCCCCGGCAC‑CCCGGGGCCCCGCCGGcaggcgc                                                                                                                 >  1:147614/1‑130 (MQ=255)
                               tCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGAt                                                                                                     >  3:23810/1‑136 (MQ=255)
                                  tGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGGGGTTCTTTTCACCATGCATTGTTGTGTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAG‑GCTGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCa                                                                                                  <  7:235275/136‑1 (MQ=255)
                                  tGCCTCTCATGCACACGCCTTTGCACCGCCGGCGGGGTTTTTTTCACCTTGCCTTTTTGTTTTGCAGTGCGGTGGCCTTGTGGCTCGGCCCCCTCCAC‑GCTGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCa                                                                                                  <  6:151623/136‑1 (MQ=255)
                                        tGATGCTCACCCCGTTGCACCGCGGGCGCGGTTTTTTTCCCCTTCCGTTGTTGTCTTGCCGGGCGGGGGCCCTGTGGCGCGGCCCCCTCCAC‑GCGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGcc                                                                                              <  4:154930/134‑1 (MQ=255)
                                           tGCGCAAGCCCTTGCACCGCCGGCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCGCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCa                                                                                         <  1:235578/136‑1 (MQ=255)
                                                   cccTTGCACCGCCGTCGCGGGTGTTTTCACCATGCCGTGAAGTCTTGCACTGCCGTGGCCCCGTTGCACGGCCCCCAGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGCCCCCGATGCAGGCCCGCCACGTCCAGc                                                                                 >  6:199973/1‑136 (MQ=255)
                                                     cTTGCACCGCCGTCGCGGTTCTTTTGACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGCCGGCTGCCGAACAGACCCACAAGCAGGCCCCCCAGGTCCACCTc                                                                               >  8:102729/1‑136 (MQ=255)
                                                          aCCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCt                                                                          >  7:9492/1‑136 (MQ=255)
                                                               ggTCGGGGTTTTTTGCACCTTCCGTTGTGGTCTTGCTGGGCGGGGGCCCTGTTGCTCGGCCCCCTCCAG‑GCTGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAg                                                                     <  4:357725/135‑1 (MQ=255)
                                                               cgggggggTTCTTTGCACCATGCGTTTAGGTGTTCCTGTGCCGGGGCGATGTTGCGCGGCCCCCTCCAC‑CCCGGGGCCGCGACGGCTGCCGAACAGACGCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAg                                                                     <  2:301221/130‑1 (MQ=255)
                                                                      gTTCTTTTCACCTGGCCTTGATGTCTTGCCGTGCGGGGGCCCGGTGGCTCGGCCCCCTGCAC‑GCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAAcc                                                              <  2:221133/136‑1 (MQ=255)
                                                                          ttttCACCATGCGTTGATGTCTTGCAGGCCGGTGGCCCTGTTGCTCGGCCCCCTCCAC‑GCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTg                                                          <  8:44838/136‑1 (MQ=255)
                                                                            ttCACCATCCCTTGATGTCTTGCCGGGCGGTGGCCCTGTGGCGCGCCCCCCTGCCC‑GCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTgg                                                         <  1:7408/135‑1 (MQ=255)
                                                                              cACCAGGCCTTGATGTCTTGCAGGGCGGGGGCCCTGTTGCGCGCCCCCCTCCAC‑GCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGAc                                                       <  1:34620/135‑1 (MQ=255)
                                                                                  aTGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGgc                                                  >  5:49599/1‑136 (MQ=255)
                                                                                   tcccTTGCTGTCTGGCCGGGCGGGGGCCGTGTGGCGCGGCCCCCGCCCG‑GCGGGGCCCGCGACGGCTGCCGAACAGACCCAGAAGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGgcg                                                 <  7:136575/134‑1 (MQ=255)
                                                                                      cTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCt                                               >  2:56968/1‑135 (MQ=255)
                                                                                      cTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCt                                               >  2:331484/1‑135 (MQ=255)
                                                                                       ttGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGCACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGCCCAGGCGCTGa                                             >  3:234575/1‑136 (MQ=255)
                                                                                                 ccggTGCCGGGGCCGTGTTGCTCGCCCCCCTCCCC‑GCGGGGCCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGc                                   <  3:147779/133‑1 (MQ=255)
                                                                                                  cAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGCCCACTTTCGcc                                  >  2:206487/1‑136 (MQ=255)
                                                                                                   aGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCACCCGCTGGACCAGGCGCTGACCACTTTCGCCg                                 >  7:384574/1‑136 (MQ=255)
                                                                                                      gcgggggCCCTGTGGCGCGCCCCCCTGCAC‑GCTGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGAcc                              <  6:307621/131‑1 (MQ=255)
                                                                                                      gCCGGGGCCCTGTTGCTCGGACCCCTGCCC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGAcc                              <  4:303940/136‑1 (MQ=255)
                                                                                                       ccGGGGCCCTGGGCCGCGGCCCCCTGCAC‑CCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCa                             <  8:123553/136‑1 (MQ=255)
                                                                                                          tGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggc                          >  2:192522/1‑136 (MQ=255)
                                                                                                          tGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggc                          >  4:281588/1‑136 (MQ=255)
                                                                                                          tGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCGGGACCAGGCGCTGCCCACTTTCGCCGACCAggc                          >  1:350363/1‑136 (MQ=255)
                                                                                                            gCCCTGTTGCGCGCCCCCCTGCCC‑GCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCTCTGACCACTTTCGCCGACCAggccg                        <  7:228467/136‑1 (MQ=255)
                                                                                                            gCCCTGTGGCTCGGCCCCCTGCAC‑GCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccg                        <  6:76814/136‑1 (MQ=255)
                                                                                                               gggTTGCTCGGCCCCCTGCCC‑GCTGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggcc                     <  8:231681/134‑1 (MQ=255)
                                                                                                               cTGTTGCTCGGCCCCCTGCAC‑GCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggcc                     <  4:23810/136‑1 (MQ=255)
                                                                                                               cTGTTGCCGGGCCCCGGCCAC‑GCTGGGCCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggcc                     <  2:323065/136‑1 (MQ=255)
                                                                                                               cTGTGGCGCGGCCCCCGCCCC‑GCGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggcc                     <  1:88777/136‑1 (MQ=255)
                                                                                                                 gtggcccGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                   <  2:273710/130‑1 (MQ=255)
                                                                                                                 gtggcccGGCCCCCGCCCC‑GCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                   <  7:348857/130‑1 (MQ=255)
                                                                                                                 gTTGCGCGCCCCCCTGCAC‑GCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                   <  2:248712/136‑1 (MQ=255)
                                                                                                                  ttGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                   <  1:331484/135‑1 (MQ=255)
                                                                                                                  ttGCTCGGCCCCCTCCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                   <  8:301571/135‑1 (MQ=255)
                                                                                                                  tggcgcGGCCCCCTGCCC‑GCGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                   <  4:347275/130‑1 (MQ=255)
                                                                                                                    gcccGGCCCCCGGCCC‑GCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc                <  8:2863/133‑1 (MQ=255)
                                                                                                                    gCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCCCCCGCTGGCCCAGGCGCTGCCCACTTTCGCCGCCCAGGCCGGCCTGCAc                >  8:285302/1‑136 (MQ=255)
                                                                                                                     cgcGGCCCCCGCCCC‑GCTGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc                <  5:382499/133‑1 (MQ=255)
                                                                                                                         gCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCACCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct           <  5:372670/136‑1 (MQ=255)
                                                                                                                         ccccccGGCAC‑CCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct           <  8:12835/135‑1 (MQ=255)
                                                                                                                          cccccTGCCCGGCGGGGCCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct           <  2:353334/136‑1 (MQ=255)
                                                                                                                             ccTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTa        >  1:305478/1‑135 (MQ=255)
                                                                                                                               tGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCACGACCTGCCAGCCCAACCGCTGGACCAGGCCCTGCCCACTATCGCCGACCAGGCCGGCCTGCACCTGCTGTacac     >  1:88817/1‑136 (MQ=255)
                                                                                                                                gCAC‑GCTGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTAcacc    <  2:165396/136‑1 (MQ=255)
                                                                                                                                  aC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGCCCACTTTCGCCGACCAGGCCGGCCTGCACCTGCGGTAcaccac  >  3:217893/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                   
ATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAC‑GCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCAC  >  ACB143/2491125‑2491390

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