Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 195,286 T→A 100% intergenic (‑97/‑82) USA300HOU_RS00900 ← / → USA300HOU_RS00905 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE195,2860TA96.7% 92.4 / ‑3.2 30intergenic (‑97/‑82)USA300HOU_RS00900/USA300HOU_RS00905hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (21/8);  total (21/9)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAAAGAAG  >  USA300TCH1516_ALE/195167‑195422
                                                                                                                       |                                                                                                                                        
gCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCt                                                                                                                     <  1:604004/141‑1 (MQ=255)
gCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCt                                                                                                                     >  1:427280/1‑141 (MQ=255)
  tACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAATATCTTAATATCCTTGt                                                                                                                  <  1:180925/141‑1 (MQ=255)
  tACGACAGATTCTTTAGCCATAATGATAGCCCCCTAAATATATGTTAGTTAACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTg                                                                                                                   <  1:495708/141‑1 (MQ=255)
    cGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATACTTGtt                                                                                                                 >  1:740178/1‑140 (MQ=255)
             cTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATaa                                                                                                        >  1:399120/1‑141 (MQ=255)
                   ccATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTATTTCCCTTGTAATCCGATAAGAATta                                                                                                  >  1:444125/1‑141 (MQ=255)
                           aTAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAata                                                                                          >  1:276679/1‑141 (MQ=255)
                             aGCCCCCTATATATATGTTTATTTACTTATACCCTAACCAAAGGTTTATACTGTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTATTCCGATAAGAATTATAGTAATAtt                                                                                        >  1:681413/1‑141 (MQ=255)
                                        tatatGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCAt                                                                             >  1:800589/1‑141 (MQ=255)
                                             ttttatttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTta                                                                        <  1:61801/140‑1 (MQ=255)
                                             gTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCTATAAGCATTATAGTAATATTTTTTCAACCATTGTta                                                                        >  1:789935/1‑141 (MQ=255)
                                               ttatttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTCCAACCATTGTtata                                                                      <  1:567741/141‑1 (MQ=255)
                                                  tttaCTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATagga                                                                   >  1:693627/1‑141 (MQ=255)
                                                       ttATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTTTAGTAATATTTTTTCAACCATTGTTAAAGGAGGCCt                                                              >  1:131246/1‑141 (MQ=255)
                                                       ttATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACAATTGTGATAAGAGGACt                                                              >  1:264936/1‑141 (MQ=255)
                                                           aCCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAAGATCCGTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCttatt                                                          >  1:385044/1‑141 (MQ=255)
                                                                   aTGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAACAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACat                                                  <  1:487109/141‑1 (MQ=255)
                                                                      aTTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGCCTTATTAATGACattat                                               >  1:97617/1‑141 (MQ=255)
                                                                            ataCTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTttat                                         >  1:41372/1‑141 (MQ=255)
                                                                               ctctTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAg                                      >  1:183544/1‑141 (MQ=255)
                                                                               ctctTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAg                                      >  1:166454/1‑141 (MQ=255)
                                                                                tctTTGAAAATATATTTTACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATCGGAGGTCTTATTAGTTACATTATTTTTATTAGa                                     >  1:234457/1‑141 (MQ=255)
                                                                                        aaTATATTTTACAGAAGTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACa                             >  1:517400/1‑141 (MQ=255)
                                                                                              ttttACAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTg                       <  1:726913/141‑1 (MQ=255)
                                                                                                  aCAGAATTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATCAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTCATGACATTATTTTTATTATAATCTAACAATCTTGAttt                   <  1:486298/141‑1 (MQ=255)
                                                                                                       aTTTTATCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCAt              >  1:182293/1‑141 (MQ=255)
                                                                                                             tCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGa        >  1:63739/1‑141 (MQ=255)
                                                                                                             tCTAAATATTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGa        >  1:674443/1‑141 (MQ=255)
                                                                                                                   aatTAAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAaagaag  <  1:800857/140‑1 (MQ=255)
                                                                                                                       |                                                                                                                                        
GCTACGACAGATTCTTTATCCATAATGATAGCCCCCTATATATATGTTTATTTACTTATACCCTAACATGATTTTTATACTCTTTGAAAATATATTTTACAGAATTTTATCTAAATATTTAAAAAAATATCTTAATATCCTTGTAATCCGATAAGAATTATAGTAATATTTTTTCAACCATTGTTATAGGAGGTCTTATTAATGACATTATTTTTATTAGAAGCTAACAATCTTGATTTTGCATCAACGAAAGAAG  >  USA300TCH1516_ALE/195167‑195422

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: