Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 2,491,272 A→C 31.5% 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,2720AC31.5% 18.3 / 4.9 33H33P (CAC→CCC) iutAFerric aerobactin receptor
Reads supporting (aligned to +/- strand):  ref base A (18/4);  new base C (3/8);  total (21/12)
Fisher's exact test for biased strand distribution p-value = 5.47e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCA  >  KT2440/2491139‑2491406
                                                                                                                                     |                                                                                                                                       
cccATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTAGGCCCCCTGCAc                                                                                                                                        >  8:104587/1‑135 (MQ=255)
     tctcATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATTCCATGATGTCTTTCAGTGCCCTGGCCCTGTTGCTCGGCCCCCTGCaacccgg                                                                                                                                   >  2:246430/1‑129 (MQ=255)
                  ttGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGGTCTTTTGACCATGCCGTGATGTGTTTCAGTGCCGTGGCCCGTGTGGTCCGCCCCCAGCCCCCCGGGGCCCCggcggcg                                                                                                                     >  2:291727/1‑134 (MQ=255)
                      cAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGGTGTTTTGACCAAGCCTTTGTGTGTTTGAGTGGCGTGGGCCAGTTGAGCGGCCCCCCGCCCCCCGGGGCCGCGCCCGCGGCCg                                                                                                                 >  1:344346/1‑135 (MQ=255)
                           aGCTAATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGGCGCGGTTATTTTCACCATGCCTTGATGTCTTGCTGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAg                                                                                                            <  8:191738/135‑1 (MQ=255)
                                              tccTCACGCCCTTCCCCCGCCGGCGGGTTTCTTTTCACCTTGCCTTGATGTCTGGCAGGGCGGGGGCCCTGTGGCGCGGCCCCCTCCACGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGccc                                                                                            <  3:280430/130‑1 (MQ=255)
                                                   aaGCCCTTGCACCGCCGTCGCGGGTGGTTTGACCCAGGCGTGGTGTGTTGCAGTGGCGGGGGCCTGTTGCTCGGCCCCCCGCCCCCCGGGGCCGCGCCGGCGGCCGCACACACCCCCCAGCAGGCCCGccccggc                                                                                    >  2:323266/1‑130 (MQ=255)
                                                             aCCGCCGGCGCGGTTCTTTCCACCATGCCTTGATGTCTTGCCGGGCCGGGGCCCTGTTGCTCGGCCCCCTGCACGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGAcc                                                                          <  4:352475/135‑1 (MQ=255)
                                                                   ggcggggTTATTTCCACCATCCTTGGAGGTCTTCCTGGGCGGGGGCCCGGTGGCGCGCCCCCCGCCCCGCGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAg                                                                    <  3:282144/130‑1 (MQ=255)
                                                                                     aTGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCACCCGCAGCACCAgg                                                  >  6:378939/1‑135 (MQ=255)
                                                                                                      aGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGc                                  >  7:502/1‑134 (MQ=255)
                                                                                                          ccGGGGCCCGGTTGCGCGCCCCCCCGCCCGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGAcc                             <  1:291727/135‑1 (MQ=255)
                                                                                                             tGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCCACCGCTGGCCCAGGCGCTGACCACTTTCGCCGCCCCgg                          >  4:329552/1‑135 (MQ=255)
                                                                                                                cccTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGACCAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGCCCAGGCGCTGCCCACGTTCGCCGCCCCggccg                       >  1:311041/1‑135 (MQ=255)
                                                                                                                cccTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGCCCACTTTCGCCGACCccgccg                       >  8:86853/1‑135 (MQ=255)
                                                                                                                cccTGTGGCGCGCCCCCCTCCCCGCCGGGCCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccg                       <  1:10554/135‑1 (MQ=255)
                                                                                                                  cTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGCCCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggc                     >  7:304482/1‑135 (MQ=255)
                                                                                                                  cTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCCGCTCGACCTGCCAGCCCAACCGCTGGCCCCGGCGCTGACCACTTTCGCCGCCCcggccggc                     >  7:70036/1‑135 (MQ=255)
                                                                                                                  cTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCCACCGCTGGCCCAGGCGCTGCCCACTTTCGCCGACCAggccggc                     >  1:328920/1‑135 (MQ=255)
                                                                                                                   tgtgGCGCGCCCCCCTGCCCGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggcc                    <  8:54300/135‑1 (MQ=255)
                                                                                                                     ttgcggggCCCCCTGCCCCCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                  <  3:19987/129‑1 (MQ=255)
                                                                                                                      ggcgcGGCCCCCGGCCCGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGc                 <  1:246430/131‑1 (MQ=255)
                                                                                                                       tgctggCCCCCCGCCCGCCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCTACCTGCCAGCCCAACCGCAGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGc                 <  7:41676/131‑1 (MQ=255)
                                                                                                                       gCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCCGGTCCCGCTCGACCTGCCAGCCCCACCGCTGGCCCAGGCGCGGACCACTTTCGCCGCCCCGGCCGGCCGGCa                >  7:156977/1‑135 (MQ=255)
                                                                                                                       gCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGCCCACTTTCGCCGACCAGGCCGGCCTGCa                >  1:227741/1‑135 (MQ=255)
                                                                                                                        cgcGGCCCCCGCCCCGCCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc               <  8:37662/133‑1 (MQ=255)
                                                                                                                        cTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCCGCCCAACCGCTGGACCAGGCGCTGACCACGATCCCCGCCCACCCcggccggccc               >  3:262383/1‑133 (MQ=255)
                                                                                                                        cTCGGCCCCCGCCACGCTGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc               <  3:93475/135‑1 (MQ=255)
                                                                                                                             cccccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCACCCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct          >  8:110292/1‑135 (MQ=255)
                                                                                                                             cccccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct          >  6:159938/1‑135 (MQ=255)
                                                                                                                             cccccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct          >  2:19791/1‑135 (MQ=255)
                                                                                                                                ccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGCCCAGGCGCGGCCCACTTTCGCCGCCCAGGCCGGCCTGCCCCTGCGGTa       >  2:325491/1‑135 (MQ=255)
                                                                                                                                ccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGt        >  3:173259/1‑134 (MQ=255)
                                                                                                                                  tGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCCCTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTaca     >  2:324022/1‑135 (MQ=255)
                                                                                                                                     aCGCAGGGGCCGCGCCGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCCACCGCTGGACCCGGCGCTGCCCACTTTCGCCGCCCCGGCCGGCCGGCACCGGCGGTAcgccc  >  1:327430/1‑134 (MQ=255)
                                                                                                                                     |                                                                                                                                       
CCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCA  >  KT2440/2491139‑2491406

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