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
*	CP000730	1,440,007	0	C	.	100.0%	118.1 / NA	25	coding (34/201 nt)	cspA	cold shock protein CspA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (14/11);  total (14/11)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCGCGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTT  >  CP000730/1439871‑1440128                                                                                                                                         |                                                                                                                          ccGCGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCCAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGt                                                                                                                      >  1:208520/1‑141 (MQ=255)       tCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAAcc                                                                                                                <  2:217410/141‑1 (MQ=255)                              aaCAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATa                                                                                         >  1:240687/1‑141 (MQ=255)                              aaCAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATa                                                                                         <  2:151641/141‑1 (MQ=255)                                cAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAATCATTTAACTGTACCTTGTTTCATAAt                                                                                       >  2:46597/1‑141 (MQ=255)                                cAGCTGGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAAt                                                                                       <  2:243880/141‑1 (MQ=255)                                    ttGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCTTTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                                                   <  2:165705/141‑1 (MQ=255)                                      gACCTTCTTCTAATGATTTGTAACCATCTTGGTTGATTGCTGTAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGaaa                                                                                 <  2:360371/141‑1 (MQ=255)                                           tcttctAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTcc                                                                            >  2:308534/1‑141 (MQ=255)                                                tAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGAc                                                                       <  2:297314/141‑1 (MQ=255)                                                           aaCCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAattcatt                                                            >  2:212396/1‑141 (MQ=255)                                                            aCCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAattcattc                                                           >  2:280905/1‑141 (MQ=255)                                                                      ttAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTa                                                 <  2:97823/141‑1 (MQ=255)                                                                            gCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCa                                           >  2:93931/1‑141 (MQ=255)                                                                                  aaaTGTACGATTACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTAc                                     <  1:39049/141‑1 (MQ=255)                                                                                   aaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACa                                    >  1:268281/1‑141 (MQ=255)                                                                                        aCGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACaaaaaa                               >  1:258583/1‑141 (MQ=255)                                                                                         cGAATACGTCATTTGCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAAt                              <  2:225222/141‑1 (MQ=255)                                                                                               cGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATGCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATGACAAAAAATACATGc                        >  1:313151/1‑141 (MQ=255)                                                                                                          ccTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTc             <  1:79935/141‑1 (MQ=255)                                                                                                          ccTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTc             <  1:320717/141‑1 (MQ=255)                                                                                                           cTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCa            >  1:151779/1‑141 (MQ=255)                                                                                                                 cTTCGATAAAGCCGAATCCTTTTT‑ATCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATAtt      >  2:251269/1‑141 (MQ=255)                                                                                                                  ttCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTc     >  2:84210/1‑141 (MQ=255)                                                                                                                     gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCttt  >  2:147152/1‑141 (MQ=255)                                                                                                                                         |                                                                                                                          CCGCGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTT  >  CP000730/1439871‑1440128

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: ‑