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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,491,2720AC43.3% 15.0 / 7.5 31H33P (CAC→CCC) iutAFerric aerobactin receptor
Reads supporting (aligned to +/- strand):  ref base A (12/4);  new base C (2/12);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.27e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCA  >  KT2440/2491161‑2491406
                                                                                                                 |                                                                                                                                      
cAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTTCAGTGCCGGGGCCCTGTTGCTCGGCCCCCTGCCCGCCGGGGCCGCGCCGGCTGCCgc                                                                                                                >  8:313238/1‑135 (MQ=255)
 aGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGGTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCAGCCCCCAGGGGCCGCGACGGCGGCCGa                                                                                                                >  4:134763/1‑135 (MQ=255)
      gCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATCTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACCGAc                                                                                                          >  8:325247/1‑136 (MQ=255)
                                       cccgcggGCGGGGTTCTTTACCCCCGGCCTTGTTGTCTGCCGGGGCGGGGCCCGTGGTCCCCCCCCCCCGCCCCGCGGGGCCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGAcc                                                                          <  1:110850/134‑1 (MQ=255)
                                            cggggcgGTTCTTTCCACCATGCCTTGATGTATTCAAGGGCCGGGGCCCTGTGGCGCGGCCCCCTGCCCGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAg                                                                    <  7:357021/132‑1 (MQ=255)
                                                   gTTCTTTTCCCCATGCCTTGAGCTCTTGCAGGGCCGGGCCCCTGTGGCTCGGCCCCCTCCACGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAAcc                                                             <  7:325247/136‑1 (MQ=255)
                                                            aCCTTGCCTTGTTGGCTGGCCGTGCCGGGGCCCGGGTGCGCGCCCCCCGGCCCGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGAc                                                      <  4:210628/134‑1 (MQ=255)
                                                             ccATGCCTTGTTGTGTCGCTGGCCGGGGCCCCTTTTGATCGCCCCCTACCCCGACGGGCCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAg                                                   <  8:95620/136‑1 (MQ=255)
                                                                  ccTTGTTGTCTTGCCGGGCGGGGGCCGTGTGCCGCGGCCCCCTCCACGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCt                                              <  7:157538/136‑1 (MQ=255)
                                                                    ttGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGCCGGCTGCCGAACAGCCCCAGATGCAGGCCCGCCAGGTCCCGCTCGCCCTGCCCGCCCCACCGCTGGCCCAGGCGCTg                                             >  5:426433/1‑135 (MQ=255)
                                                                    ttGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCCGGTCCAGCTCGCCCTGCCAGCCCCACCGCAGGACCAGgcgcgg                                             >  5:283051/1‑133 (MQ=255)
                                                                               cAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGCCCAGGCGCTGACCACTTTCGc                                  >  2:89087/1‑135 (MQ=255)
                                                                               cAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGCCCACGTTCGc                                  >  3:422841/1‑135 (MQ=255)
                                                                               cAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGc                                  >  3:221275/1‑135 (MQ=255)
                                                                               cAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGCCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCCGCCCAACCGCTGGACCAGCCGCTGACCACTTTCGc                                  >  2:274196/1‑135 (MQ=255)
                                                                                     cggggCCGTGTTGCTCGGCCCCCTCCTCGCGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCa                            <  5:345483/132‑1 (MQ=255)
                                                                                           cccTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccg                       >  4:106969/1‑135 (MQ=255)
                                                                                           cccTGTGGCGCGCCCCCCCCCCCGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccg                       <  3:114165/135‑1 (MQ=255)
                                                                                               tGTTGCTCGGCCCCCTGCCGGCTGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggcc                    <  1:55231/135‑1 (MQ=255)
                                                                                                gTTGCGCGGCCCCCTGCCCGCGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                  <  7:443077/136‑1 (MQ=255)
                                                                                                 ttGCGCGGCCCCCTGCCCGCGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                  <  6:347465/135‑1 (MQ=255)
                                                                                                  ggCTCGGCCCCGGCCCGGCGGGGCCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGc                 <  1:288263/134‑1 (MQ=255)
                                                                                                    cgcGGCCCCCGGCAGGCGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc               <  4:180293/133‑1 (MQ=255)
                                                                                                    cTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGCCCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc               >  2:109088/1‑135 (MQ=255)
                                                                                                       ggCCCCCTGCACGCCGGGGCCGCGACGGCTGCCGAACAGGCCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctg            <  3:114628/135‑1 (MQ=255)
                                                                                                        gCCCCGCGCCCGCGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCCGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct          <  8:197282/136‑1 (MQ=255)
                                                                                                         ccccgTGCCCGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct          <  1:389795/135‑1 (MQ=255)
                                                                                                         cccccTGCCCCCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct          <  6:334989/135‑1 (MQ=255)
                                                                                                         cccccTGCACGCAGGGGCCGCGCCGGCTGCCGAACAGACCCAGATGCAGGCCCGCCCGGTCCAGCTCGACCTGCCAGCCCCCCCGCTGGCCCCGGCGCTGCCCACTTTCGCCGACCAGGCCGGCCGGCCCcggcg          >  3:17894/1‑134 (MQ=255)
                                                                                                         cccccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCGCGCCCTGCCAGCCCCCCCGCGGGCCCCGGCGCGGACCACTTTCGCCGCCCAGGCCGGCCTGCCCcggct          >  5:210628/1‑135 (MQ=255)
                                                                                                           cccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCCGGTCCAGCTCGACCTGCCAGCCCACCCGCTGGACCAGGCGCGGCCCCCTTTCGCCGCCCCGGCCGGCCGGCACCGGCTGt        >  2:389795/1‑135 (MQ=255)
                                                                                                                 aCGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCACCCGCTGGACCAGGCGCTGCCCCCTATCGCCGCACCGGCCGGCCTGCACCGGCTGtggccca  >  5:48080/1‑129 (MQ=255)
                                                                                                                 |                                                                                                                                      
CAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCA  >  KT2440/2491161‑2491406

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: