Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,891,904 (T)6→5 coding (4/402 nt) rot ← HTH‑type transcriptional regulator rot

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,891,8990T.97.5% 176.5 / ‑1.7 40coding (9/402 nt)rotHTH‑type transcriptional regulator rot
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (19/20);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

TTTTGCCATAAAGTTAGTAAAATTAAAATTTCTTCTCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATCCCAACAATCCCAAAACTTGTATGTGCTAACTTATGCATAAAACTACAAGTGTAAATAAA  >  USA300TCH1516_ALE/1891762‑1892027
                                                                                                                                         |                                                                                                                                
ttttGCCATAAAGTTAGTAAAATTAAAATTTCTTCTCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACtttt                                                                                                                               >  2:156476/1‑141 (MQ=255)
    gCCATAAAGTTAGTAAAATTAAAATTTCTTCTCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATCCCAAATACAGTGTCGTTATTTACTTTTTTCa                                                                                                                           >  1:350218/1‑141 (MQ=255)
       aTAAAGTTAGTAAAATTAAAATTTCTTCTCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCt                                                                                                                       >  1:98349/1‑141 (MQ=255)
       aTAAAGTTAGTAAAATTAAAATTTCTTCTCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTNTCATGCt                                                                                                                       >  2:310426/1‑141 (MQ=255)
                     attaaaattTCTTCTCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAAt                                                                                                         >  2:373073/1‑141 (MQ=255)
                      ttaaaattTCTTCTCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCACATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAAtt                                                                                                        <  2:413501/141‑1 (MQ=255)
                         aaattTCTTCTCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATtaat                                                                                                     <  1:183166/141‑1 (MQ=255)
                         aaattTCTTCTCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATtaat                                                                                                     <  2:379354/141‑1 (MQ=255)
                         aaattTCTTCTCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATtaat                                                                                                     <  1:389068/141‑1 (MQ=255)
                              tcttctCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATca                                                                                                <  2:52514/141‑1 (MQ=255)
                              tcttctCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATca                                                                                                <  2:231872/141‑1 (MQ=255)
                               cttctCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATcac                                                                                               <  2:98385/141‑1 (MQ=255)
                                         cATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGct                                                                                     >  1:27183/1‑141 (MQ=255)
                                         cATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGct                                                                                     <  1:296539/141‑1 (MQ=255)
                                                 aTTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACa                                                                             <  1:158513/141‑1 (MQ=255)
                                                 aTTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACa                                                                             <  2:375961/141‑1 (MQ=255)
                                                 aTTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACa                                                                             <  2:218181/141‑1 (MQ=255)
                                                    cGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGt                                                                          >  1:7072/1‑141 (MQ=255)
                                                    cGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGt                                                                          >  2:126426/1‑141 (MQ=255)
                                                    cGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGt                                                                          >  2:155049/1‑141 (MQ=255)
                                                           aaTCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTaaca                                                                   >  1:291796/1‑141 (MQ=255)
                                                           aaTCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGATTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTaaca                                                                   >  1:218705/1‑141 (MQ=255)
                                                             tctcGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTaacaaa                                                                 >  2:356923/1‑141 (MQ=255)
                                                                cGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACAt                                                              >  2:215783/1‑141 (MQ=255)
                                                                cGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACAt                                                              >  1:117284/1‑141 (MQ=255)
                                                                cGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACAt                                                              >  2:8021/1‑141 (MQ=255)
                                                                cGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACAt                                                              >  1:22937/1‑141 (MQ=255)
                                                                cGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACAt                                                              >  2:28927/1‑141 (MQ=255)
                                                                 gCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATc                                                             >  2:378217/1‑141 (MQ=255)
                                                                              ttAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATCCCAACAATCCCaa                                                <  2:27183/141‑1 (MQ=255)
                                                                                  tGTCACCCAAACGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATCCCAACAATCCCAAAACt                                            <  2:159261/141‑1 (MQ=255)
                                                                                     cACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACCGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATCCAAGTATTAACAAACACCCCAACAATCCCAAAACTTGt                                         >  2:327597/1‑141 (MQ=255)
                                                                                     cACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATCCCAACAATCTCAAAACTTGt                                         >  2:307123/1‑141 (MQ=255)
                                                                                                  ttCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTGTTTCATGCTAAACATCTCCCAATTAATTATCACATATTTAGCTACTATACCAGTATTAACAAACATCCCAACAATCCCACAACTTGTATGTGCTAACTTa                            <  1:402965/141‑1 (MQ=255)
                                                                                                        ttGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATCCCAACAATCCCAAAACTTGTATGTGCTAACTTATTCATa                      <  1:72212/141‑1 (MQ=255)
                                                                                                        ttGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATCCCAACAATCCCAAAACTTGTATGTGCTAACTTATGCATa                      <  2:231073/141‑1 (MQ=255)
                                                                                                          gCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATCCCAACAATCCCAAAACTTGTATGTGCTAACTTATGCATaaa                    >  2:155900/1‑141 (MQ=255)
                                                                                                                    aaaTACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATCCCAACAATCCCAAAACTTGTATGTGCTAACTTATGCATAAAACTACAAgtg          <  1:28927/141‑1 (MQ=255)
                                                                                                                    aaaTACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATCCCAACAATCCCAAAACTTGTATGTGCTAACTTATGCATAAAACTACAAgtg          <  1:8021/141‑1 (MQ=255)
                                                                                                                    aaaTACAGTGTCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATCCCAACAATCCCAAAACTTGTATGTGCTAACTTATGCATAAAACTACAAgtg          <  2:22937/141‑1 (MQ=255)
                                                                                                                            tgtCGTTATTTACTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATCCCAACAATCCCAAAACTTGTATGTGCTAACTTATGCATAAAACTACAAGTGtaaataaa  <  1:373073/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                
TTTTGCCATAAAGTTAGTAAAATTAAAATTTCTTCTCTAGACATTTTGTATTCGCTTTCAATCTCGCTGAAAATTGAGTTAATGTCACCCAAAAGTGTTTCTAATTGCAAAATTCCAAATACAGTGTCGTTATTTACTTTTTTCATGCTAAACATCTCCCAATTAATAATCACATATTTAGCTACTATACAAGTATTAACAAACATCCCAACAATCCCAAAACTTGTATGTGCTAACTTATGCATAAAACTACAAGTGTAAATAAA  >  USA300TCH1516_ALE/1891762‑1892027

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