Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,440,007	Δ1 bp	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	.	96.9%	141.7 / ‑1.1	32	coding (34/201 nt)	cspA	cold shock protein CspA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (16/15);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

CGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATA  >  CP000730/1439874‑1440134                                                                                                                                      |                                                                                                                                cgGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTaa                                                                                                                         >  1:14552/1‑141 (MQ=255)   gTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAAc                                                                                                                       <  1:292470/141‑1 (MQ=255)    tCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAAcc                                                                                                                      <  1:244661/141‑1 (MQ=255)          tCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTaa                                                                                                                <  1:371192/141‑1 (MQ=255)             actactTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTg                                                                                                             >  2:274321/1‑141 (MQ=255)               tactTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTa                                                                                                           >  2:334758/1‑141 (MQ=255)                     aaaCTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGt                                                                                                     >  1:305974/1‑141 (MQ=255)                        cTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTc                                                                                                  <  1:350840/141‑1 (MQ=255)                             cAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATATAGCCGAATACTTTTT‑AGCTTTAAACCATTTAACTGTACGTTGTTTCATAAt                                                                                             >  2:259109/1‑141 (MQ=39)                                  tGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGaa                                                                                        >  1:262413/1‑141 (MQ=255)                                                    ttGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAatt                                                                      >  2:360273/1‑141 (MQ=255)                                                      gTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAattc                                                                     >  1:191072/1‑141 (MQ=255)                                                       tAACCATCGTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTCTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAattcat                                                                   >  2:145213/1‑141 (MQ=255)                                                               ttGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATg                                                           <  2:245400/141‑1 (MQ=255)                                                                tGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGc                                                          >  1:330438/1‑141 (MQ=255)                                                                 ggTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCt                                                         <  1:227476/141‑1 (MQ=255)                                                                 ggTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCt                                                         >  1:275554/1‑141 (MQ=255)                                                                    tAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGcgaat                                                      <  1:234155/141‑5 (MQ=255)                                                                        tGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑TGCTTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGc                                                  >  2:264350/1‑141 (MQ=255)                                                                                 aTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACaa                                         <  1:289262/141‑1 (MQ=255)                                                                                       gAATACGTCATTTTCTCCTTCAACTTCGATAAGGCCGAATCCTTTTT‑AGCGTTAAACCGTTTATCTGTACCTTGTTTCATAATCTGAAACCTTAAAGATTAATATTCATTCAATATGCTTTTTCGCATATTACAAAAAATa                                   <  2:157633/141‑1 (MQ=255)                                                                                         aTACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACa                                 <  1:74463/141‑1 (MQ=255)                                                                                           aCGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATg                               <  1:50076/141‑1 (MQ=255)                                                                                                        cTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCa                  <  1:369871/141‑1 (MQ=255)                                                                                                               ttCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAGATATTc           <  1:289220/141‑1 (MQ=255)                                                                                                                 cGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCtt         <  2:37042/141‑1 (MQ=255)                                                                                                                 cGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCtt         <  2:349401/141‑1 (MQ=255)                                                                                                                  gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTTTTAATCTGAAACCTCCAAGACTAAAATTCATTCAATATACTTATTCGCATATGACAAAAAATACATGTAAATTTCACTCAAATAttatgt        >  2:152910/1‑137 (MQ=255)                                                                                                                  gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCttt        >  2:106692/1‑141 (MQ=255)                                                                                                                   aTAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCtttt       >  2:175405/1‑140 (MQ=255)                                                                                                                     aaaGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCa     >  1:9627/1‑141 (MQ=255)                                                                                                                        gCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAata  >  2:148173/1‑141 (MQ=255)                                                                                                                                      |                                                                                                                                CGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATA  >  CP000730/1439874‑1440134

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