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

CCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCA  >  KT2440/2491139‑2491406
                                                                                                                                     |                                                                                                                                       
cccATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTAGGCCCCCTGCAc                                                                                                                                        >  8:104587/1‑135 (MQ=255)
     tctcATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATTCCATGATGTCTTTCAGTGCCCTGGCCCTGTTGCTCGGCCCCCTGCaacccgg                                                                                                                                   >  2:246430/1‑129 (MQ=255)
                  ttGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGGTCTTTTGACCATGCCGTGATGTGTTTCAGTGCCGTGGCCCGTGTGGTCCGCCCCCAGCCCC‑CCGGGGCCCCggcggcg                                                                                                                     >  2:291727/1‑134 (MQ=255)
                      cAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGGTGTTTTGACCAAGCCTTTGTGTGTTTGAGTGGCGTGGGCCAGTTGAGCGGCCCCCCGCCCC‑CCGGGGCCGCGCCCGCGGCCg                                                                                                                 >  1:344346/1‑135 (MQ=255)
                           aGCTAATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGGCGCGGTTATTTTCACCATGCCTTGATGTCTTGCTGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAg                                                                                                            <  8:191738/135‑1 (MQ=255)
                                              tccTCACGCCCTTCCCCCGCCGGCGGGTTTCTTTTCACCTTGCCTTGATGTCTGGCAGGGCGGGGGCCCTGTGGCGCGGCCCCCTCCACG‑CCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGccc                                                                                            <  3:280430/130‑1 (MQ=255)
                                                   aaGCCCTTGCACCGCCGTCGCGGGTGGTTTGACCCAGGCGTGGTGTGTTGCAGTGGCGGGGGCCTGTTGCTCGGCCCCCCGCCCC‑CCGGGGCCGCGCCGGCGGCCGCACACACCCCCCAGCAGGCCCGccccggc                                                                                    >  2:323266/1‑130 (MQ=255)
                                                             aCCGCCGGCGCGGTTCTTTCCACCATGCCTTGATGTCTTGCCGGGCCGGGGCCCTGTTGCTCGGCCCCCTGCACG‑CCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGAcc                                                                          <  4:352475/135‑1 (MQ=255)
                                                                   ggcggggTTATTTCCACCATCCTTGGAGGTCTTCCTGGGCGGGGGCCCGGTGGCGCGCCCCCCGCCCCG‑CGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAg                                                                    <  3:282144/130‑1 (MQ=255)
                                                                                     aTGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCACCCGCAGCACCAgg                                                  >  6:378939/1‑135 (MQ=255)
                                                                                                      aGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGc                                  >  7:502/1‑134 (MQ=255)
                                                                                                          ccGGGGCCCGGTTGCGCGCCCCCCCGCCCG‑CCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGAcc                             <  1:291727/135‑1 (MQ=255)
                                                                                                             tGGCCCTGTTGCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCCACCGCTGGCCCAGGCGCTGACCACTTTCGCCGCCCCgg                          >  4:329552/1‑135 (MQ=255)
                                                                                                                cccTGTTGCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGACCAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGCCCAGGCGCTGCCCACGTTCGCCGCCCCggccg                       >  1:311041/1‑135 (MQ=255)
                                                                                                                cccTGTTGCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGCCCACTTTCGCCGACCccgccg                       >  8:86853/1‑135 (MQ=255)
                                                                                                                cccTGTGGCGCGCCCCCCTCCCCG‑CCGGGCCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccg                       <  1:10554/135‑1 (MQ=255)
                                                                                                                  cTGTTGCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGCCCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggc                     >  7:304482/1‑135 (MQ=255)
                                                                                                                  cTGTTGCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCCGCTCGACCTGCCAGCCCAACCGCTGGCCCCGGCGCTGACCACTTTCGCCGCCCcggccggc                     >  7:70036/1‑135 (MQ=255)
                                                                                                                  cTGTTGCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCCACCGCTGGCCCAGGCGCTGCCCACTTTCGCCGACCAggccggc                     >  1:328920/1‑135 (MQ=255)
                                                                                                                   tgtgGCGCGCCCCCCTGCCCG‑CCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggcc                    <  8:54300/135‑1 (MQ=255)
                                                                                                                     ttgcggggCCCCCTGCCCC‑CCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                  <  3:19987/129‑1 (MQ=255)
                                                                                                                      ggcgcGGCCCCCGGCCCG‑CCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGc                 <  1:246430/131‑1 (MQ=255)
                                                                                                                       tgctggCCCCCCGCCCGCCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCTACCTGCCAGCCCAACCGCAGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGc                 <  7:41676/131‑1 (MQ=255)
                                                                                                                       gCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCCGGTCCCGCTCGACCTGCCAGCCCCACCGCTGGCCCAGGCGCGGACCACTTTCGCCGCCCCGGCCGGCCGGCa                >  7:156977/1‑135 (MQ=255)
                                                                                                                       gCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGCCCACTTTCGCCGACCAGGCCGGCCTGCa                >  1:227741/1‑135 (MQ=255)
                                                                                                                        cgcGGCCCCCGCCCCG‑CCGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc               <  8:37662/133‑1 (MQ=255)
                                                                                                                        cTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCCGCCCAACCGCTGGACCAGGCGCTGACCACGATCCCCGCCCACCCcggccggccc               >  3:262383/1‑133 (MQ=255)
                                                                                                                        cTCGGCCCCCGCCACG‑CTGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc               <  3:93475/135‑1 (MQ=255)
                                                                                                                             cccccTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCACCCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct          >  8:110292/1‑135 (MQ=255)
                                                                                                                             cccccTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct          >  6:159938/1‑135 (MQ=255)
                                                                                                                             cccccTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct          >  2:19791/1‑135 (MQ=255)
                                                                                                                                ccTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGCCCAGGCGCGGCCCACTTTCGCCGCCCAGGCCGGCCTGCCCCTGCGGTa       >  2:325491/1‑135 (MQ=255)
                                                                                                                                ccTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGt        >  3:173259/1‑134 (MQ=255)
                                                                                                                                  tGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCCCTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTaca     >  2:324022/1‑135 (MQ=255)
                                                                                                                                     aCG‑CAGGGGCCGCGCCGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCCACCGCTGGACCCGGCGCTGCCCACTTTCGCCGCCCCGGCCGGCCGGCACCGGCGGTAcgccc  >  1:327430/1‑134 (MQ=255)
                                                                                                                                     |                                                                                                                                       
CCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACG‑CAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCA  >  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: ‑