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% 80.0 / NA 18coding (34/201 nt)cspA_2Cold shock protein CspA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (11/7);  total (11/7)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAAT  >  USA300TCH1516_ALE/1439878‑1440140
                                                                                                                                 |                                                                                                                                     
cGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCa                                                                                                                           >  1:115282/1‑141 (MQ=255)
 gCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCCTTTTCTCCTTCAACTTCGATAAACCCGACTCCTTTTTAGCGTTAAAccct                                                                                                                          >  2:22757/1‑139 (MQ=255)
                aaaCTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGt                                                                                                           >  1:61022/1‑141 (MQ=255)
                               gccTTCTTCTGATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGCACGAATACGTCATTTTCTCCTTCAACTTCGATAACGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAAc                                                                                            <  2:231759/140‑1 (MQ=255)
                                                   aaCCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCTTAATCTGAAACCTCCAAGACTAAAattcatt                                                                        >  2:40296/1‑141 (MQ=255)
                                                           tGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGc                                                                >  1:153034/1‑141 (MQ=255)
                                                              ttAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTa                                                             <  1:161759/141‑1 (MQ=255)
                                                                       gAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCatat                                                    <  1:113258/141‑1 (MQ=255)
                                                                          aaaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTAc                                                 >  2:230157/1‑141 (MQ=255)
                                                                                 cGAATACGTCATTTCCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTaaaaaaaaa                                          >  1:149062/1‑140 (MQ=255)
                                                                                         tCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCaa                                  <  1:36894/141‑1 (MQ=255)
                                                                                                    ttCAACGTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCaa                       <  2:97610/141‑1 (MQ=255)
                                                                                                     acAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAACACTCCAAGAATAAAATTCATTCAATATGCTTATAAAAAAAAAAAAAAAAATACATGCAAATTTCACTCaaa                      <  2:101169/140‑1 (MQ=255)
                                                                                                            cGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCtt               <  1:188931/141‑1 (MQ=255)
                                                                                                            cGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATAAATGCAAATTTCACTCAAATATtttt               >  2:225127/1‑141 (MQ=255)
                                                                                                                  aGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCCat         >  1:85849/1‑141 (MQ=255)
                                                                                                                    ccGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAatat       >  1:88373/1‑141 (MQ=255)
                                                                                                                         tCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaat  >  2:193160/1‑141 (MQ=255)
                                                                                                                                 |                                                                                                                                     
CGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAAT  >  USA300TCH1516_ALE/1439878‑1440140

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