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

CGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAAT  >  CP000730/1439878‑1440140                                                                                                                                  |                                                                                                                                      cGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCa                                                                                                                           >  1:115282/1‑141 (MQ=255)  gCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCCTTTTCTCCTTCAACTTCGATAAACCCGACTCCTTTTT‑AGCGTTAAAccct                                                                                                                          >  2:22757/1‑139 (MQ=255)                 aaaCTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGt                                                                                                           >  1:61022/1‑141 (MQ=255)                                gccTTCTTCTGATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGCACGAATACGTCATTTTCTCCTTCAACTTCGATAACGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAAc                                                                                            <  2:231759/140‑1 (MQ=255)                                        ttttcagcaattAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGa                                                                                    <  1:127589/130‑1 (MQ=255)                                                    aaCCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCTTAATCTGAAACCTCCAAGACTAAAattcatt                                                                        >  2:40296/1‑141 (MQ=255)                                                            tGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGc                                                                >  1:153034/1‑141 (MQ=255)                                                               ttAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTa                                                             <  1:161759/141‑1 (MQ=255)                                                                        gAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCatat                                                    <  1:113258/141‑1 (MQ=255)                                                                           aaaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTAc                                                 >  2:230157/1‑141 (MQ=255)                                                                                  cGAATACGTCATTTCCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTaaaaaaaaa                                          >  1:149062/1‑140 (MQ=255)                                                                                          tCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCaa                                  <  1:36894/141‑1 (MQ=255)                                                                                                     ttCAACGTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCaa                       <  2:97610/141‑1 (MQ=255)                                                                                                      acAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAACACTCCAAGAATAAAATTCATTCAATATGCTTATAAAAAAAAAAAAAAAAATACATGCAAATTTCACTCaaa                      <  2:101169/140‑1 (MQ=255)                                                                                                             cGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCtt               <  1:188931/141‑1 (MQ=255)                                                                                                             cGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATAAATGCAAATTTCACTCAAATATtttt               >  2:225127/1‑141 (MQ=255)                                                                                                                   aGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCCat         >  1:85849/1‑141 (MQ=255)                                                                                                                     ccGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAatat       >  1:88373/1‑141 (MQ=255)                                                                                                                          tCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaat  >  2:193160/1‑141 (MQ=255)                                                                                                                                  |                                                                                                                                      CGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAAT  >  CP000730/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: ‑