Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 591,178 C→T 42.4% A465V (GCA→GTA)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE591,1780CT42.4% 11.6 / 25.0 28A465V (GCA→GTA) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (8/8);  new base T (5/7);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TAGACTTGGTATTCAAGCATTTGAACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGTACTTGGTAACTATTACCTTACTTTAGAAAGAAAAGAT  >  USA300TCH1516_ALE/591041‑591308
                                                                                                                                         |                                                                                                                                  
tAGACTCGGTATTCAAGCATTTACACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCCCTTGTAACAACATCCTACAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGcaca                                                                                                                                 <  2:264275/141‑1 (MQ=255)
        gTATTCAAGCATTTGAACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAg                                                                                                                         <  2:338593/141‑1 (MQ=255)
           ttCAAGCATTTGAACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCaa                                                                                                                      <  1:101106/141‑1 (MQ=255)
            tCAAGCATTTGAACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCACGAGAGGTACAACCTGACGCCAg                                                                                                                     >  2:169421/1‑141 (MQ=255)
                gCATTTGAACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATg                                                                                                                 <  2:281656/141‑1 (MQ=255)
                        aCCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCGCAAGCTGAAGCAAGAATGTTGATGtt                                                                                                         >  2:175235/1‑141 (MQ=255)
                               ttAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTAAAGCAAGAATGTTGATGTTAgccgca                                                                                                  >  2:290172/1‑141 (MQ=255)
                                 aGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGcaca                                                                                                <  2:319994/141‑1 (MQ=255)
                                    tGAAGGTCGTGCGATTCGTTTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACaaaa                                                                                             <  2:29333/141‑1 (MQ=255)
                                    tGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACaaaa                                                                                             <  2:194706/141‑1 (MQ=255)
                                       aGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACAt                                                                                          >  1:85865/1‑141 (MQ=255)
                                       aGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACAt                                                                                          <  1:331519/141‑1 (MQ=255)
                                        ggTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATc                                                                                         <  1:124280/141‑1 (MQ=255)
                                            gTGCGATTCGTCTACATCCACTTGTAACAACCGATTATAAAGCTGACTTTGACGGTGCCCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCACGAATGTTGATGTTAGCAGCACAAAACATCTTGa                                                                                     >  1:74323/1‑141 (MQ=255)
                                            gTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGa                                                                                     >  1:68762/1‑141 (MQ=255)
                                            gTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGa                                                                                     >  1:150780/1‑141 (MQ=255)
                                                gATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAAccc                                                                                 <  1:238172/141‑1 (MQ=255)
                                                  ttCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACACAACATCTTGAACCCTa                                                                               >  1:101271/1‑141 (MQ=255)
                                                                  tGTAAAAACAGCGTATAACGCTTACTTTGACTGTGACCAAATGGCTGTTCACGTTCCTTTATCAATAGGAGTACACGCTGAAGCAAGAATGTTGATGTTAGCAGCATAACACATCTTGATCCCTAAAGATGGTAAACCTGt                                                               <  2:310781/141‑1 (MQ=255)
                                                                   gTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTa                                                              <  1:173662/141‑1 (MQ=255)
                                                                         acaGCTTATAACGCTCACTTTGACGGTGACCAAATGGCTGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTGAACCTGTAGTTaca                                                        <  2:270112/141‑1 (MQ=255)
                                                                           aGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACAcc                                                      <  1:212223/141‑1 (MQ=255)
                                                                                ataACGCTGAGTTTGACGATGACCAAATGGCGGGTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCACGAATGTTCATGTTAGCAGCACAAACCAGCTTGAACCCTAAATATGGCAAACCTGTAGTTACACCATcac                                                 >  2:132155/1‑141 (MQ=255)
                                                                                      cTGACTTTGAGGGTGAAAAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAACTATGTTGATGTTAGTAGCACAAAACAGCTTGAACCCTAAAGATGGTAAACCTGTAGTAACACCATCACAAGata                                           <  2:184867/141‑1 (MQ=255)
                                                                                          cTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATGTTCATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGt                                       <  1:182878/141‑1 (MQ=255)
                                                                                                                gTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGTACTTGGTAACTATTACCTTACt                 >  2:20588/1‑141 (MQ=255)
                                                                                                                 ttCACGTTCCTTTATCAAAAGAGGTACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGTACTTGGTAACTATTACCTTACtt                >  1:181500/1‑141 (MQ=255)
                                                                                                                         ccTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTGGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGTACTTGGTAACTATTACCTTACTTTagaaaga        >  2:50873/1‑141 (MQ=255)
                                                                                                                               tCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGTACTTGGTAACTATTACCTTACTTTAGAAAGAAAAGAt  >  2:229100/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                  
TAGACTTGGTATTCAAGCATTTGAACCAACTTTAGTTGAAGGTCGTGCGATTCGTCTACATCCACTTGTAACAACAGCTTATAACGCTGACTTTGACGGTGACCAAATGGCGGTTCACGTTCCTTTATCAAAAGAGGCACAAGCTGAAGCAAGAATGTTGATGTTAGCAGCACAAAACATCTTGAACCCTAAAGATGGTAAACCTGTAGTTACACCATCACAAGATATGGTACTTGGTAACTATTACCTTACTTTAGAAAGAAAAGAT  >  USA300TCH1516_ALE/591041‑591308

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: