Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 722,264 (A)8→7 coding (235/798 nt) fepC → Ferric enterobactin transport ATP‑binding protein FepC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE722,2570A.92.6% 100.7 / 0.6 27coding (228/798 nt)fepCFerric enterobactin transport ATP‑binding protein FepC
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (18/7);  total (18/9)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTAGTTTCATATGGTCGTTTTCCACATCAAAAAGGATTTGGTAGATTAACTGCTGAGGATAAGAAAGAAATTGAT  >  USA300TCH1516_ALE/722119‑722392
                                                                                                                                          |                                                                                                                                       
aGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATATACCATCTACGAAAGAGATTGCaca                                                                                                                                       >  1:158584/1‑139 (MQ=255)
        cAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCaaaaaaatagc                                                                                                                               <  1:246355/141‑5 (MQ=255)
               tGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTa                                                                                                                       >  1:177389/1‑141 (MQ=255)
                           cTGCGGGAAATCTACTTTGCTAAAGGCAGTGTAACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAGAGCCATTTTACCTCAATCACCt                                                                                                           <  2:117099/141‑1 (MQ=255)
                                 gAAATCTACTGTGCGAAAGGCATTGTCACGTTTATTGGAAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTa                                                                                                     <  2:71824/141‑1 (MQ=255)
                                              cTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTaa                                                                                        >  1:132560/1‑141 (MQ=255)
                                              cTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTaa                                                                                        >  1:212863/1‑141 (MQ=255)
                                                aaaGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCGTTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACt                                                                                      <  1:122339/141‑1 (MQ=255)
                                                  aGGCATTGACACGATTAGTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGt                                                                                    <  1:167071/141‑1 (MQ=255)
                                                   ggCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGtt                                                                                   >  1:203636/1‑141 (MQ=255)
                                                   ggCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGtt                                                                                   >  1:239613/1‑141 (MQ=255)
                                                       ttGTAACGTTTATTGCGAGGTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTgggg                                                                               <  2:75586/141‑1 (MQ=255)
                                                        tGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGa                                                                              >  2:72928/1‑141 (MQ=255)
                                                            aCGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTa                                                                          >  2:215481/1‑141 (MQ=255)
                                                                ttATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTAGttt                                                                      >  1:169751/1‑141 (MQ=255)
                                                                  aTTGGCAGTTAACGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTAGGCAACTAGTTTCa                                                                    >  2:149712/1‑141 (MQ=255)
                                                                  aTTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTAGTTTCa                                                                    >  2:149706/1‑141 (MQ=255)
                                                                   ttGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTAGTTTCat                                                                   <  1:3158/141‑1 (MQ=255)
                                                                     ggCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGATTAACTGTTAGGGATTTATTTTCatat                                                                 >  2:30503/1‑141 (MQ=255)
                                                                                      aaGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTAGTTTCATATGGTCGTTTTCCACATCa                                                >  1:37553/1‑141 (MQ=255)
                                                                                           tttttAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTAGTTTCATATGGTCGTTTTCCACATCAAAAAg                                           >  2:164912/1‑141 (MQ=255)
                                                                                           tttttAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTAGTTTCATATGGTCGTTTTCCACATCAAAAAg                                           >  1:253677/1‑141 (MQ=255)
                                                                                           tttttAGATGGGTAAAATATTCATACACAAACTACGAAACAGATTGAAAAAAAAGAGTCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTAGTTTCATATGGTCGTTTTCCACATCAAAAAg                                           <  2:138457/141‑1 (MQ=255)
                                                                                            tgtgaGAAGGTGAAAATATTCAGATGCAATCTACGAAAGAGATTGAAAAAAAAGAGCCATTTTACCTAAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTGGTTTCATATGGTCGTTTTCCACATCAAAAAgg                                          <  2:248830/137‑1 (MQ=255)
                                                                                                    ggTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCGTAACTGTTGGGGAATTAGTTTCATATGGTCGTTTTCCACATCAAAAAGGATTTGGtg                                  >  2:141970/1‑140 (MQ=255)
                                                                                                       gAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATAACCTGAATAAGCATATGGTTTAACCGTCGGGGAATTATTTTCATATGGTCGTTTCCCACATCAAAAAGGATTTGGCAgct                               >  2:200967/1‑139 (MQ=255)
                                                                                                                             cGAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTAGTTTCATATGGTCGTTTTCCACATCAAAAAGGATTTGGTAGATTAACTGCTGAGGATacgaaaaa         >  2:9597/1‑139 (MQ=255)
                                                                                                                             agAAAGAGATTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTAGTTTCATATGGTCGTTTTCCACATCAAAAAGGATTTGGTAGATTAACTGCTGAGGATaagaaaga         <  2:278796/140‑1 (MQ=255)
                                                                                                                                    gaTTGCAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTAGTTTCATATGGTCGTTTTCCACATCAAAAAGGATTTGGTAGATTAACTGCTGAGGATAAGAAAGAAattgat  >  2:181188/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                       
AGTGACGTCAATCATTGGTCCTAACGGCTGCGGGAAATCTACTTTGCTAAAGGCATTGTCACGTTTATTGGCAGTTAAAGAAGGCGAAGTATTTTTAGATGGTGAAAATATTCATACACAATCTACGAAAGAGATTGCAAAAAAAATAGCCATTTTACCTCAATCACCTGAAGTAGCAGATGGCTTAACTGTTGGGGAATTAGTTTCATATGGTCGTTTTCCACATCAAAAAGGATTTGGTAGATTAACTGCTGAGGATAAGAAAGAAATTGAT  >  USA300TCH1516_ALE/722119‑722392

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