Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,440,085 G→T 100% intergenic (‑45/+126) cspA_2 ← / ← USA300TCH1516_01338 Cold shock protein CspA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,440,0850GT97.1% 87.4 / ‑5.6 28intergenic (‑45/+126)cspA_2/USA300TCH1516_01338Cold shock protein CspA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (16/11);  total (17/11)
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.89e-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.

AAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTAATC  >  USA300TCH1516_ALE/1439950‑1440216
                                                                                                                                       |                                                                                                                                   
aaaaaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCatat                                                                                                                                <  1:543712/141‑1 (MQ=255)
      gTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACaaa                                                                                                                          <  1:756673/141‑1 (MQ=255)
         cGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACaaaaaa                                                                                                                       >  1:271066/1‑141 (MQ=255)
         cGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACaaaaaa                                                                                                                       >  1:442841/1‑141 (MQ=255)
            aTACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATAc                                                                                                                    <  1:265317/141‑1 (MQ=255)
            aTACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCAGAAACCTCCAAGACTAAAATTCATTCAATATTCTTATTCTCATCTTACAAAAAATAc                                                                                                                    >  1:855525/1‑141 (MQ=255)
                 tCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAACAAAtacataca                                                                                                               >  1:262049/1‑141 (MQ=255)
                                    cGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACAGCCAAATTTCACTCAAATATTCt                                                                                            <  1:865902/141‑1 (MQ=255)
                                              gAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACAAAAATTCATTCAATATGCTTGTTCTCCTATTACAAAAAATACATGCAACTTTACCTCAAATATTCTTTGACATATg                                                                                  >  1:282596/1‑141 (MQ=255)
                                                        tCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAatcat                                                                        <  1:384461/141‑1 (MQ=255)
                                                                aaaCCATTTAACTGTATCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTaa                                                                >  1:380977/1‑141 (MQ=255)
                                                                 aaCCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAAc                                                               >  1:833028/1‑141 (MQ=255)
                                                                  aCCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACt                                                              >  1:344376/1‑141 (MQ=255)
                                                                         aaCTGTACCTTGTTTCATAATCTGAAACCCCCCAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGtatta                                                       <  1:237832/141‑1 (MQ=255)
                                                                          aCTGTACCTTGTTTCGTAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGtattat                                                      >  1:645899/1‑141 (MQ=255)
                                                                          aCTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGtattat                                                      >  1:625799/1‑141 (MQ=255)
                                                                          aCTGTACCTTGTTTCAAAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGtattat                                                      >  1:431421/1‑141 (MQ=255)
                                                                                ccTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCa                                                <  1:608305/141‑1 (MQ=255)
                                                                                ccTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCa                                                <  1:722054/141‑1 (MQ=255)
                                                                                     tttCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGttat                                           <  1:787757/141‑1 (MQ=255)
                                                                                           aaTCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCGGCTTAACTTGTATTATAGTGCAGTTATTACCaa                                     >  1:716146/1‑141 (MQ=255)
                                                                                              cTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAAtt                                  >  1:131746/1‑141 (MQ=255)
                                                                                                             cTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAAtaataa                   <  1:397219/141‑1 (MQ=255)
                                                                                                              tAAAATTCATTCAATCTGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATaaa                  >  1:437281/1‑141 (MQ=255)
                                                                                                                     cattcaATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATAtt           >  1:501932/1‑141 (MQ=255)
                                                                                                                     cattcaATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATAtt           >  1:638439/1‑141 (MQ=255)
                                                                                                                       ttcaATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCa         >  1:517543/1‑141 (MQ=255)
                                                                                                                              tGCATATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTAATc  <  1:502298/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
AAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTAATC  >  USA300TCH1516_ALE/1439950‑1440216

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