Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 396,540 G→A A208T (GCA→ACA)  efeM → putative iron uptake system component EfeM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE396,5400GA96.9% 104.6 / NA 32A208T (GCA→ACA) efeMputative iron uptake system component EfeM
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/14);  minor base C (1/0);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATG  >  USA300TCH1516_ALE/396406‑396672
                                                                                                                                      |                                                                                                                                    
tAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTa                                                                                                                         >  1:507459/1‑148 (MQ=255)
   atatCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGa                                                                                                                      <  1:565155/148‑1 (MQ=255)
            cAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaa                                                                                                             <  1:562962/148‑1 (MQ=255)
             aaaaTTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaaa                                                                                                            <  1:454668/148‑1 (MQ=255)
                           aaGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaa                                                                                              >  1:180393/1‑148 (MQ=255)
                             ggTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaa                                                                                            >  1:255733/1‑148 (MQ=255)
                                        cccaTTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATaaaaaa                                                                                 <  1:509862/145‑1 (MQ=255)
                                         cTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAAt                                                                                <  1:335059/148‑1 (MQ=255)
                                           aTTAAATGAAGTTGCAACTTCTAAAATGACAGGGGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACAAAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTa                                                                              <  1:606263/148‑1 (MQ=255)
                                               aaTGAAGGTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtg                                                                          <  1:417234/148‑1 (MQ=255)
                                                aTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtga                                                                         >  1:323001/1‑148 (MQ=255)
                                                                      cacaGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGCAGGCACACATAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAGAAAATTACGGGCTGATATCCAAATGAACTTAGAt                                                   >  1:163269/1‑148 (MQ=255)
                                                                          ggTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAg                                               <  1:591582/148‑1 (MQ=255)
                                                                           gTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGGTGATATCCAAATGAACTTCGATAAAGt                                              >  1:38302/1‑148 (MQ=255)
                                                                             gaagaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCCAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCGAAATGAACTTCGATAAAGTGa                                            >  1:439464/1‑148 (MQ=255)
                                                                               agaagaaATTTATTCACATACAGAATTATATGATTTTAAAGCGAACATTGAAGGCACACAAAAAATTTATGACGTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAAt                                          <  1:143811/148‑1 (MQ=255)
                                                                                  agaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATGTCCAAATGAACTTCGATAAAGTGAATCaa                                       >  1:522425/1‑148 (MQ=255)
                                                                                  agaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCaa                                       <  1:322016/148‑1 (MQ=255)
                                                                                    aaaTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAGGTGTTGCTTTACAAATGAACTTCGATAAAGTGAATCAatt                                     >  1:240200/1‑148 (MQ=255)
                                                                                             cacaTACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATaa                            <  1:214300/148‑1 (MQ=255)
                                                                                                 tACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAatat                        <  1:505715/148‑1 (MQ=255)
                                                                                                  aCAGATTTATATGATTTTAAAACGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAAtata                       <  1:218522/148‑1 (MQ=255)
                                                                                                     ggatcaaatGATTTTAAAGGCAAGAGCGAAGGCCCTCAAACACTTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGACGATATCCAAATGAACTTCGATGAAGTGAATCAATTATTGGATAAAAATAAAg                    >  1:468292/8‑148 (MQ=255)
                                                                                                            atGATTTTAAAGCGAACGTTGAAGGCACGCAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTAGAAATAGTAAGACAATTATTGGATAAATATAAAGATaacaa             >  1:391739/1‑148 (MQ=255)
                                                                                                            atGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATaacaa             <  1:90540/148‑1 (MQ=255)
                                                                                                              gATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAACAAGATGCAAAATTAAGTGATGATATCCAAATGAACTTCGGTAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcg           >  1:34097/1‑148 (MQ=255)
                                                                                                               aTTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAGGATAACAAcgg          >  1:127253/1‑148 (MQ=255)
                                                                                                               aTTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcgg          >  1:42187/1‑148 (MQ=255)
                                                                                                               aTTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAAcgg          >  1:247227/1‑148 (MQ=255)
                                                                                                                      aGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTAt   >  1:346816/1‑148 (MQ=255)
                                                                                                                      aGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTAt   >  1:220273/1‑148 (MQ=255)
                                                                                                                       gCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATAAATTATTGGATAAATATAAAGATAACAACGGCGGTTATg  >  1:576293/1‑148 (MQ=255)
                                                                                                                                      |                                                                                                                                    
TAGATATCATACCAAAATTAATGTTACAAGGTTCTGTTGACCTATTAAATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTGAATCAATTATTGGATAAATATAAAGATAACAACGGCGGTTATG  >  USA300TCH1516_ALE/396406‑396672

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