Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,245,923 C→T 6.3% E29K (GAA→AAA)  rpoE ← putative DNA‑directed RNA polymerase subunit delta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,245,9230CT 6.3% 176.9 / 2.3 63E29K (GAA→AAA) rpoEputative DNA‑directed RNA polymerase subunit delta
Reads supporting (aligned to +/- strand):  ref base C (19/40);  new base T (2/2);  total (21/42)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

CTCCAACATTTAAAAAACGACCATCTGTGTTTAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCATTGCCCAATACTTCCTTTCGAAGGATATAATCTAAATTTATACGTCT  >  USA300TCH1516_ALE/2245791‑2246053
                                                                                                                                    |                                                                                                                                  
ctcCAACATTAAAAAAACGACCATCTGTGTTTAAATCCGTGTAACATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATCATCACATATAAGTTCATTGTTTCg                                                                                                                                   >  2:62532/1‑134 (MQ=255)
 tcCAACATTTAAAAAACGACCATCTGTGTTTAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTctc                                                                                                                                  >  1:286337/1‑132 (MQ=255)
 tcCAACATTTAAAAAACGACCATCTGTGTTTAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGc                                                                                                                                  <  1:534226/134‑1 (MQ=255)
  ccAACATTTAAAAAACGACCATCTGACGTTAAATCCGTGTAAAATGGTACAACACGATTGTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTTGcc                                                                                                                                 <  1:31543/134‑1 (MQ=255)
          tAAAAAACTACCATCTGTGTTTAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCAt                                                                                                                         <  2:552088/134‑1 (MQ=255)
          tAAAAAACGACCATCTGTGTTTAAATCCGTGTAATATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCAt                                                                                                                         <  2:63177/134‑1 (MQ=255)
          tAAAAAACGACCATCTGTGTTTAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCAt                                                                                                                         <  1:588830/134‑1 (MQ=255)
          tAAAAAACGACCATCTGTGTTTAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCAt                                                                                                                         <  2:120856/134‑1 (MQ=255)
          tAAAAAACGACCATCTGTGTTTAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCAt                                                                                                                         <  1:495848/134‑1 (MQ=255)
           aaaaatCGATCTTCTGTGTTTAACTCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCAtt                                                                                                                        <  1:248695/134‑1 (MQ=255)
            aaaaaCGACCATCTGTGTTTAAATCCGTGTAAAATTGTACAACACGATATTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTc                                                                                                                       >  2:493479/1‑134 (MQ=255)
                 cGACCATCTGTGTTTAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCaataa                                                                                                                  >  1:117563/1‑134 (MQ=255)
                  caCCATCTGTGTTTGCATCCGTGTAAAATTTTACACCACGATTTCCAATTTCTTCGTACTCATAATCACCTAACTCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCaataat                                                                                                                 <  2:346767/133‑1 (MQ=255)
                   aCCAGCTGTGTTTAGGTCCGTGTAAACTTGTACAACACGATTTTCAATGTCTTCGTACTCATAATCACCTCACGTTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATg                                                                                                                <  2:204625/134‑1 (MQ=255)
                    ccATCTGTGTTTAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATACTCACCTAACGCTCTAAATTCATCGATGATATAATAAAAGTACAATGTTTCGCCTTTATCATTCAATAATGt                                                                                                               >  1:489366/1‑134 (MQ=255)
                    ccATCTGTGTTTAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTATAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGt                                                                                                               >  1:507618/1‑134 (MQ=255)
                      ttttgtgTTTAAATCCATGTTAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGtat                                                                                                             <  1:242518/131‑1 (MQ=255)
                       tCTGTGTTTAAATCCGTGTACAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGtata                                                                                                            <  2:146806/134‑1 (MQ=255)
                          gtgtTTAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCATTATCATACAATAATGTATAAGc                                                                                                         <  2:531117/134‑1 (MQ=255)
                          gtgtTAAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGc                                                                                                         <  1:192550/134‑1 (MQ=255)
                                      gtgtAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATTATCACCTAACGCTCTACAATTGTTGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATaaa                                                                                             <  2:5970/134‑1 (MQ=255)
                                        gtAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATg                                                                                           <  1:478112/134‑1 (MQ=255)
                                            aagtGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGAttt                                                                                       <  2:536002/131‑1 (MQ=255)
                                            aaTTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGAttt                                                                                       <  1:361932/134‑1 (MQ=255)
                                            aaTTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGAttt                                                                                       <  1:464521/134‑1 (MQ=255)
                                             aTTGTACAACACGATTTTCAATTTCTTCGTACTCATAATTATCTAACGCTTTGAATTCATCGATGATATCATATAAGTGCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGAtttt                                                                                      <  1:98626/134‑1 (MQ=255)
                                             aTTGTACAACACGATTTTCAATTTCTTCGTACTCATAATTATCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGAtttt                                                                                      <  1:471085/134‑1 (MQ=255)
                                             aTTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAGATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGAtttt                                                                                      <  2:582730/134‑1 (MQ=255)
                                             aTTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGAtttt                                                                                      <  1:466046/134‑1 (MQ=255)
                                                 tACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTtcat                                                                                  <  2:353881/134‑1 (MQ=255)
                                                  acaacaCGATTTTCAATTTCTTCGTACTCATAATCACCTAACGTTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTTGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTtcatc                                                                                 <  1:167173/134‑1 (MQ=255)
                                                    aacaCGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTtcatcat                                                                               >  2:278711/1‑133 (MQ=255)
                                                     acacGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAAc                                                                              <  1:167391/134‑1 (MQ=255)
                                                     acacGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAAc                                                                              <  1:313919/134‑1 (MQ=255)
                                                        ctaTTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCAt                                                                           <  1:429221/132‑1 (MQ=255)
                                                        cGATTTTCAATTTCTTCGTACTCATAATCACCTATCGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCAt                                                                           <  1:359568/134‑1 (MQ=255)
                                                        cGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCAt                                                                           <  1:173342/134‑1 (MQ=255)
                                                           ttttCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTTGCCTTTATCATTCAATAATGTATAAGCCATATCCATACATGATTTTTTATCAACTAtttt                                                                        >  1:12194/1‑133 (MQ=255)
                                                           ttttCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTTGCCTTTATCATTCAATAATGTATAAGCCATATCACTAACTCATTCTTCATCTACCAtttg                                                                        >  1:127697/1‑134 (MQ=255)
                                                             ttCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCCACCAtttgtt                                                                      >  2:315708/1‑134 (MQ=255)
                                                               cAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCAtttgtttt                                                                    >  2:446026/1‑134 (MQ=255)
                                                                aaTTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCAtttgttttg                                                                   <  2:454148/134‑1 (MQ=255)
                                                                     cttcGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATGTATGATTTTTCATCAACCCTTTCTTTTGTATaa                                                              >  1:35027/1‑134 (MQ=255)
                                                                     cttcGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTTTTTTGTATaa                                                              >  1:579177/1‑134 (MQ=255)
                                                                     cttcGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATaa                                                              >  1:294258/1‑134 (MQ=255)
                                                                     cttcGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATAACATATAAGTTCATTATTTCGCCTTTATCATTCAATAATGTTTATGCGGCACCAATCATTGCTTTTTCCACCACGATTTGTGTTGTATaa                                                              >  1:136745/1‑134 (MQ=255)
                                                                     cttcGTACTCAAAATCAGTTGGTGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATaa                                                              >  2:370564/1‑134 (MQ=255)
                                                                             tCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAAt                                                      <  1:525964/134‑1 (MQ=255)
                                                                                  aTCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAACTGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCa                                                 >  1:516986/1‑134 (MQ=255)
                                                                                  aTCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATAGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCa                                                 >  1:318777/1‑134 (MQ=255)
                                                                                   tCACCTAACGCTCTAAATTCATCGATGATATCATATATGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTCGAATTTTCAt                                                <  1:309950/134‑1 (MQ=255)
                                                                                   tCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCAt                                                <  1:523788/134‑1 (MQ=255)
                                                                                         aaCGCTCTAAATTCATCGATGATATCATACAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTCTAATCTTGAATTTTCATTGCCCa                                          >  2:197781/1‑134 (MQ=255)
                                                                                                tAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCATTGCCCAATACTTc                                   <  1:551551/134‑1 (MQ=255)
                                                                                                tAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCATTGCCCAATACTTc                                   <  1:131992/134‑1 (MQ=255)
                                                                                                    ttCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATTTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCATTGCCCAATACTTCCttt                               <  1:148833/134‑1 (MQ=255)
                                                                                                      cttCGATGATATCATATATGTTCCTTGTTTCGCCTTTATCATTCAATAATGTATTAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCATTGCCCAATACTTCCTTTCg                             <  2:316993/132‑1 (MQ=255)
                                                                                                       atctttgatATCATATAATTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTCTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCATTGCCCAATACTTCCTTTCGa                            <  2:305104/129‑1 (MQ=255)
                                                                                                       aTCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCATTGCCCAATACTTCCTTTCGa                            >  2:12794/1‑134 (MQ=255)
                                                                                                                       aaGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCATTGCCCAATACTTCCTTTCGAAGCATATATTCTaact            >  1:270305/1‑132 (MQ=255)
                                                                                                                              ttGTTTCGCTTTTGTCATTCTATATTGTATAGTCCATATCAATAAATGTTTTTTCATCAACCATTTGTTTTGTATAATCTCGATTTTTCATTGCCCATTACTTCCTTTCGAAGGATATAATCTAATTTTATACg     <  2:297653/134‑1 (MQ=255)
                                                                                                                              ttGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCATTGCCCAATACTTCCTTTCGAAGGATATAATCTAAATTTATACg     <  1:446026/134‑1 (MQ=255)
                                                                                                                              ttGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTCTTGTATAATCTTGAATTTTCATTGCCCAATACTTCCTTTCGAAGGATATAATCTAAATTTATACg     <  2:363599/134‑1 (MQ=255)
                                                                                                                                 tttCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACAATTTGTTTTGTATAATCTTGAATTTTCATTGCCCAATACTTCCTTTCGAAGGATATAATCTAAATTTATACGTCt  >  2:262044/1‑134 (MQ=255)
                                                                                                                                    |                                                                                                                                  
CTCCAACATTTAAAAAACGACCATCTGTGTTTAAATCCGTGTAAAATTGTACAACACGATTTTCAATTTCTTCGTACTCATAATCACCTAACGCTCTAAATTCATCGATGATATCATATAAGTTCATTGTTTCGCCTTTATCATTCAATAATGTATAAGCCATATCAATAAATGATTTTTCATCAACCATTTGTTTTGTATAATCTTGAATTTTCATTGCCCAATACTTCCTTTCGAAGGATATAATCTAAATTTATACGTCT  >  USA300TCH1516_ALE/2245791‑2246053

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: