Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 2,491,272 A→C 30.2% H33P (CAC→CCC)  iutA → Ferric aerobactin receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,491,2720AC30.2% 43.7 / 7.0 40H33P (CAC→CCC) iutAFerric aerobactin receptor
Reads supporting (aligned to +/- strand):  ref base A (22/5);  new base C (3/10);  total (25/15)
Fisher's exact test for biased strand distribution p-value = 1.04e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCAC  >  ACB122/2491139‑2491407
                                                                                                                                      |                                                                                                                                       
cccATTCTCATTTACGCCTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGGTGGTTTCCCCCAGCCGTGATGTGTGTGAGTGCCGGGGCCC‑AGTTTCGGCGCCCCCGGCAc                                                                                                                                        >  6:97414/1‑135 (MQ=255)
     tctcATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGGTCGTTTGACCAATCCGTGATGTGTTGCAGTGCCGTGGCCC‑GGTTGCTCGGCCCCCAGCACCCAggg                                                                                                                                  >  7:237807/1‑136 (MQ=255)
             aCGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGGTGTTTTTAACAATCCTTGTTTTGTTGGAGTGGCGTGGGCC‑GTGTGGTCGGCCCCCCGCCCCCCCGCGCCGCGAc                                                                                                                          >  7:281978/1‑136 (MQ=255)
              cGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGGTGTTTTAACCATGCCGTGATGTATTGCAGTGCCGTGGGCC‑TGTTTCTCGGCCCCCCGCCCGCAGGGGCCGCGCCg                                                                                                                         >  4:89941/1‑136 (MQ=255)
                      cAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTTACCATGCCGTGATGTCTTGGAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCGGCACGCAGGGGCCCCGACGGCTGCCgc                                                                                                                 >  7:263840/1‑135 (MQ=255)
                      cAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGCCGCGGGTCTTTTGAACCAGCCGTGATGTATTGCACTGGCGTGGCCC‑TGTTGCACCGCCCCCAGCACGCCGGGGCCGCGCCGGCTGCCg                                                                                                                  >  5:296209/1‑135 (MQ=255)
                                          cgcatccTCAACCCCTTGCACCCCCGCGGGGGTTCTTACCACCATCCCTTGTGGTATTCCCGGGCGGGGGCCG‑TGTTCCCCGGCCCCCTCCCCCCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGccc                                                                                             <  8:4343/130‑1 (MQ=255)
                                                      cccTTCCCCCGCCGGCGCGTTTCTTTACCCCATCCCTGTATGGCTGGCCGGGCCGGGCCCC‑TGTGGCGCGGCCCCCGCCCCGCGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGc                                                                                 <  7:50361/136‑1 (MQ=255)
                                                       ccTTCCCCCGCCGTCGGGGTTGTTACCACCATGCATTGTTGGCTTGACGGGCGGGGGCCG‑TGTTGCGCGCCCCCGTCCCCGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCt                                                                                <  8:281978/136‑1 (MQ=255)
                                                             aCCGCGGGCGGGGTTCTTTTCACCCTGCCTTGCTGTC‑TGCCTGGCCGGGGCCCGTGTGGCTCGGCCCCCGGCCCGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCt                                                                          <  7:24754/136‑1 (MQ=255)
                                                                  cggcggggTTCTTACCACCCTGCCTTGAGGGCTGGCGGGGCGGTGGCCC‑TGTGGCGGGGCCCCCTCCCGGCCGGGGCCGCGACGGCTGCCGAACAGACACAGATTCAGGCCCGCCAGGTCCAGCTCGACCTGCCAg                                                                     <  4:227092/130‑1 (MQ=255)
                                                                    tCGCGGTTCGTTTCACCCTGCCGTGATGTCTTGCAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCCCGCCGGGGCCGCGACGGCGGGCGCACAGCCCCCGATGCCGGCCCGCCCGGTCCAGCTCGCCCTGCCAGcc                                                                   >  4:385864/1‑136 (MQ=255)
                                                                               ttCCCCAGGCCTTGATGTCTCGCCGGGCCGGGGCCC‑TGTTGCTCGGCCCCCACCACGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTgg                                                         <  4:273735/135‑1 (MQ=255)
                                                                               ttCACCATGCCTTGATGTCTGGCAGGGCCGGGGCCC‑TGTTGCGCGCCCCCCGCCACGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTgg                                                         <  8:326693/135‑1 (MQ=255)
                                                                                       gCCTTGATGTCTTGCAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCCGCCCAACCGCGGGCCCCGgcg                                                 >  5:294801/1‑135 (MQ=255)
                                                                                       gCCTTGATGTCTTGCAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGgcg                                                 >  1:88378/1‑135 (MQ=255)
                                                                                         cTTGATGTCTTGCAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTg                                              >  4:268509/1‑136 (MQ=255)
                                                                                          ttGATGTCTTGCAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGCACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCACCCGCTGGCCCAGGCGCTg                                              >  1:275206/1‑135 (MQ=255)
                                                                                          ttGATGTCTTGCAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCCGCTCGACCTGCCAGCCCAACCGCGGGACCAGGCGCTg                                              >  1:331350/1‑135 (MQ=255)
                                                                                                   tGCAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGCCCACTTTc                                     >  1:209947/1‑135 (MQ=255)
                                                                                                     cAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGCCCTGCCCGCCCAACCGCGGGCCCCCGCGCTGACCCCATTCGc                                   >  2:154386/1‑135 (MQ=255)
                                                                                                     cAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGCCCACTTTCGcc                                  >  8:91564/1‑136 (MQ=255)
                                                                                                     cAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGc                                   >  5:163235/1‑135 (MQ=255)
                                                                                                       gTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGCCCAGGCGCTGACCACTTTCGCCGa                                >  8:145878/1‑136 (MQ=255)
                                                                                                          ccgggcccc‑TGTTGCTCGGCCCCCTCCCCGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCa                             <  8:345150/130‑1 (MQ=255)
                                                                                                                ccc‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCCGCCCACCCGCAGGACCAGGCGCGGACCACTTTCGCCGCCCCcgccgg                       >  7:343332/1‑136 (MQ=255)
                                                                                                                ccc‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGCCCACTTTCGCCGCCCAggccg                        >  2:278646/1‑135 (MQ=255)
                                                                                                                ccc‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGCCCACTTTCGCCGACCAggccgg                       >  3:115146/1‑136 (MQ=255)
                                                                                                                  c‑TGTTGCGCGGCCCCCTCCACGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggcc                     <  4:26783/136‑1 (MQ=255)
                                                                                                                       tGCGCGGCCCCCGGCCCGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                   <  6:141108/134‑1 (MQ=255)
                                                                                                                        gcgcGGCCCCCGCCCCGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc                <  3:148535/133‑1 (MQ=255)
                                                                                                                        gCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGCTGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCAACCCCCTGGCCCCGGCCCTGACCACTTTCGCCGACCACGCCGGCCGGCAc                >  8:200826/1‑136 (MQ=255)
                                                                                                                        gCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGCCCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc                >  4:380351/1‑136 (MQ=255)
                                                                                                                        gCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGCCCACATTCGCCGCCCAGGCCGGCCTGCa                 >  1:251847/1‑135 (MQ=255)
                                                                                                                        gCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGCCCAGGCCGGCCTGccc                >  7:345150/1‑134 (MQ=255)
                                                                                                                         cgcgcccccccGCCGGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc                <  5:97414/130‑1 (MQ=255)
                                                                                                                         cgcGGCCCCCGCCACCCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc                <  2:255416/133‑1 (MQ=255)
                                                                                                                             ccccccGGCACGCTGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct           <  8:237807/135‑1 (MQ=255)
                                                                                                                                   tGCCCGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTaca      <  5:287878/135‑1 (MQ=255)
                                                                                                                                   tGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCGGGCCCCGGCGCTGACCCCGTTCGCCGCCCCGGCCGGCCGCCCCCTGCTGGacc      >  1:151254/1‑134 (MQ=255)
                                                                                                                                      aCGCAGGGGCCGCGACGGCTGCCGAACAGCCCCAGATGCAGGCCCGCCCGGTCCAGCTAGACCAGCCAGCCCAACCGCCGGCCCAGGCGCGGACCCCGTTCGCCGACCCGGCCGGCCTGCACCAGCTGTAcaccac  >  8:392810/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                       
CCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCC‑TGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCAC  >  ACB122/2491139‑2491407

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