Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,064,254 Δ1 bp coding (235/1218 nt) ywtF ← Putative transcriptional regulator YwtF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,064,2530C.93.8% 66.4 / ‑1.5 16coding (236/1218 nt)ywtFPutative transcriptional regulator YwtF
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (8/7);  total (8/8)
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.92e-01

TGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACGATTAGCATTTTTGCTGTCTTGAG  >  USA300TCH1516_ALE/1064108‑1064376
                                                                                                                                                 |                                                                                                                           
tgcatgGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCttt                                                                                                                          <  1:64367/149‑1 (MQ=255)
                       aGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTatcat                                                                                                  <  1:263592/149‑1 (MQ=255)
                       aGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTatcat                                                                                                  >  1:320925/1‑149 (MQ=255)
                          cAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGa                                                                                               >  1:293470/1‑149 (MQ=255)
                          cAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGa                                                                                               >  1:326468/1‑149 (MQ=255)
                             cTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATAc                                                                                            <  1:308247/149‑1 (MQ=255)
                              tttAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATAcc                                                                                           <  1:247659/149‑1 (MQ=255)
                              tttAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATAcc                                                                                           <  1:58348/149‑1 (MQ=255)
                                           cttaTAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATa                                                                              >  1:121894/4‑149 (MQ=255)
                                                                     cATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGAc                                                    >  1:123984/1‑149 (MQ=255)
                                                                          tAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTttat                                               <  1:107422/149‑1 (MQ=255)
                                                                                          ttgattAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGa                               >  1:163862/1‑149 (MQ=255)
                                                                                                                   gCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACGATTAGCATTTTTGCTGttt      >  1:12134/1‑147 (MQ=255)
                                                                                                                   gCGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACGATTAGCATTTTTGCTGTCt      <  1:294596/149‑1 (MQ=255)
                                                                                                                    cGTCTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACGATTAGCATTTTTGCTGTCtt     >  1:176704/1‑149 (MQ=255)
                                                                                                                       cTGATCTAGAGTGTTCAGCATCCTGTCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACGATTAGCATTTTTGCTGTCTTgag  <  1:214780/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                           
TGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTAGTTAAAATCATAGCGTCTGATCTAGAGTGTTCAGCATCCTGTCCTTTTTTTCTTCTTCCATCGTTATCATCGATACCTAAGAAAAGAATAGAGATAGGTTGTTCTTCGGGATTGACTTTATTATCTCTTAAGTTGGATTGACGATTAGCATTTTTGCTGTCTTGAG  >  USA300TCH1516_ALE/1064108‑1064376

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