Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,440,007 Δ1 bp 100% coding (34/201 nt) cspA_2 ← Cold shock protein CspA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,440,0070C.100.0% 43.2 / NA 10coding (34/201 nt)cspA_2Cold shock protein CspA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (6/4);  total (6/4)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGCGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAAT  >  USA300TCH1516_ALE/1439872‑1440140
                                                                                                                                       |                                                                                                                                     
cgcgGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGtt                                                                                                                                 <  2:38472/141‑1 (MQ=255)
                            aaCAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATACTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATa                                                                                                     >  1:90140/1‑141 (MQ=255)
                               aGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTTAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATc                                                                                                  <  1:55507/141‑1 (MQ=255)
                                           ttctAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCaa                                                                                      >  1:110766/1‑141 (MQ=255)
                                                 tGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTaa                                                                                <  1:136632/141‑1 (MQ=255)
                                                                      aaTTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTAtt                                                           >  2:3098/1‑141 (MQ=255)
                                                                       aTTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTc                                                          >  1:96870/1‑141 (MQ=255)
                                                                                                                     tAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGc            >  1:49987/1‑141 (MQ=255)
                                                                                                                       aaGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCaa          >  1:113260/1‑141 (MQ=255)
                                                                                                                               tCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTCAAATTCATTCAATATGCTTATTCGAATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaat  <  1:83743/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CGCGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAAT  >  USA300TCH1516_ALE/1439872‑1440140

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