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

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

CCGCGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGA  >  CP000730/1439871‑1440138
                                                                                                                                        |                                                                                                                                   
ccGCGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGt                                                                                                                                <  2:184864/141‑1 (MQ=255)
                         aCTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGttt                                                                                                       >  2:5225/1‑141 (MQ=255)
                                                  tGATGTGTACGCAGCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTAGTTTCGCAGGCTGAAACCTCCAAGACTaa                                                                              <  2:82327/141‑1 (MQ=255)
                                                       tGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAattc                                                                         <  1:22168/141‑1 (MQ=255)
                                                                                   aTGTCCGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACaa                                             <  1:134429/141‑1 (MQ=255)
                                                                                   aTGTACGAATACGTCATTTTCTCCTTCGACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAACCCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACaa                                             <  1:320560/141‑1 (MQ=255)
                                                                                                tCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCaa                                <  1:329519/141‑1 (MQ=255)
                                                                                                  ctttgattcgtCAACTTCGATACAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAATCCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAAt                              <  2:171855/131‑1 (MQ=255)
                                                                                                         ccTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTc                       <  1:327284/141‑1 (MQ=255)
                                                                                                            tCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCaaa                    <  2:190127/141‑1 (MQ=255)
                                                                                                             cAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACACAAAATACATGCAAATTTCACTCAAat                   >  1:235098/1‑141 (MQ=255)
                                                                                                                    gATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCttt            <  1:19626/141‑1 (MQ=255)
                                                                                                                            cGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATCTTCTTTGCAATATg    >  2:50033/1‑141 (MQ=255)
                                                                                                                            cGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTCCATAATCTGAAACCTCCAACACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATg    <  2:248115/141‑1 (MQ=255)
                                                                                                                              aaTCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACCAAAATATTCTTTGCAATATGGa  >  2:152884/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                   
CCGCGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGA  >  CP000730/1439871‑1440138

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: