Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,047,275	Δ1 bp	7.3%	intergenic (+146/‑28)	cspG → / → ymcE	cold shock protein homolog, cold‑inducible/cold shock gene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,047,275	0	A	.	7.3%	176.4 / 7.0	55	intergenic (+146/‑28)	cspG/ymcE	cold shock protein homolog, cold‑inducible/cold shock gene
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base . (2/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

GACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTT  >  CP009273/1047001‑1047518                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    gccattttCCACCANTTCACCGCCCCACAGCGCCATGAATTCCGCACGATGAACGAAAATCAGAAAGTTGAATGTTCTATTGAGCATGGGCAACGTGGACCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATa                                                                                                                                                                                                                                                     <  2:550215/271‑1 (MQ=255)                  aacagacATCAACAGCAATGAATCACACACCCTGAACGAAAATCAGAAAGTTTAATTTTCTATTGAGCAGGGTCCACGTGGCCCCGCGGCCGCCAACTTTGTTACGCTCTTATGTTGCCATTATTACTCCACATCTCCTTTTCCGCTGTCCATGTTGTCATTGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATa                                                                                                                                                                                                                                    <  2:505161/269‑1 (MQ=255)                          tCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGCTCACCTTTGATAAGGTCGCCTTGCAAATACGCACAGa                                                                                                                                                                                                                            >  2:33988/1‑275 (MQ=255)                                    ttaaTTCCGAACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCTAACCTTTTTACGCTCTAAGGTTGCCATTATTACTCAACATCTTCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGCCCATGcgccg                                                                                                                                                                                                                  <  1:736741/273‑1 (MQ=255)                                           cGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGAt                                                                                                                                                                                                            >  1:279351/1‑274 (MQ=255)                                                gCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACAGACCATGCGCCGCTGGATTTcac                                                                                                                                                                                                       <  1:33989/274‑1 (MQ=255)                                                                                                         gggaGCGAANGTTGTTACCCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCaa                                                                                                                                             <  2:279346/272‑1 (MQ=255)                                                                                                            ctcggaCGTTGNTAAGCTCTGAGGTTAAAATTAGTACTCACCATCTACAGTTCTGCCGTCAAGGTGATAATGGTTAGCAGTACCGCACGTCTGCATTCGGTGGGACGGAGCCAAACCCTTTGGGCGCTACGTGTCTTTTTTGGAAATCGACGATGATCACCTTTGACAACGTCGCGCTGTAAATACGCACTGACCATGCGCCGCTTGATTACACAAAATAATATCAGGCTCCCTCGTGGATCCTTCTTTATATCTGCCTTATTTTTCTGCAACGc                                                                                                                                          <  2:66978/270‑1 (MQ=255)                                                                                                            aGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGc                                                                                                                                          <  2:296936/275‑1 (MQ=255)                                                                                                                             tctAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATCCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTAATATCTGCCTTAGTTTTCTTCAACGTtgcatctttttttttct                                                                                                                         >  2:593017/1‑263 (MQ=255)                                                                                                                               tAAGGTTTCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGAt                                                                                                                       <  1:708030/275‑1 (MQ=255)                                                                                                                               tAAGGTTGCCATTATTACTCAACATCTCAATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCNCTGGATTTCACACAATAATATCAGGCTCCCt                                                                                                                                                                                   >  1:552019/1‑215 (MQ=255)                                                                                                                                 agtttGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATGTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGt                                                                                                                              <  1:593030/263‑1 (MQ=255)                                                                                                                                   gtggCCATTATTACTCAACATCTCCAGTTCCGCTGTCCATGTTGTCATGGGTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCgag                                                                                                                                                                               <  2:552008/212‑3 (MQ=255)                                                                                                                                    ttGTCATTATTTCTTAACTTCTCCATTTCCGCTGTCCCTGTTGTCATTGTTTACAGTACCGCACCTCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCt                                                                                                                  <  1:154109/275‑1 (MQ=255)                                                                                                                                                  cAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAAt                                                                                                    >  2:721960/1‑275 (MQ=255)                                                                                                                                                     cATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGAt                                                                                                                                                                                                            >  1:448648/1‑168 (MQ=255)                                                                                                                                                         tcCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTca                                                                                                                                                                                                        <  2:448641/168‑1 (MQ=255)                                                                                                                                                            aTTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATa                                                                                                                                                                                                                                                     <  2:718293/120‑1 (MQ=255)                                                                                                                                                             tttCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTcac                                                                                                                                                                                                       >  1:293269/1‑165 (MQ=255)                                                                                                                                                              ttaagatgtaaaTGTTGTAATGGTTATCAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGgaga                                                                                        <  2:572044/264‑1 (MQ=255)                                                                                                                                                               tCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGCTGTGACGGAGCGAAACCCTTTGGTCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTctc                                                                                                                                                                                                                                              >  2:509119/1‑122 (MQ=255)                                                                                                                                                                 cGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACaaaa                                                                                                                                                                                                   <  2:293264/165‑1 (MQ=255)                                                                                                                                                                   cTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGAAGATGCTCACCTTTGATAACGTCGCGCTg                                                                                                                                                                                                                                          <  1:509128/124‑1 (MQ=255)                                                                                                                                                                     gTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAGTATCAGGCTCCCTCGTGGANCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCggg                                                                                            >  1:572057/1‑264 (MQ=255)                                                                                                                                                                      tCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGa                                                                                                                        >  1:412809/1‑235 (MQ=255)                                                                                                                                                                                cATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGTGGAGATGGCTTATAACATTTCt                                                                       >  1:30853/1‑274 (MQ=255)                                                                                                                                                                                cATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTAATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGa                                                                                                                        >  1:142699/1‑225 (MQ=255)                                                                                                                                                                                         caGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTTTTGATTATTGAATATTTCGGGGAGAGGGCTTATAACATTTCTTACCTAAcc                                                              >  2:175613/1‑274 (MQ=255)                                                                                                                                                                                            tACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAatat                                                                                                                                                                                             >  2:335969/1‑144 (MQ=255)                                                                                                                                                                                                gCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAgg                                                                                                                                                                                         <  1:335975/144‑1 (MQ=255)                                                                                                                                                                                                gCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTAcc                                                      <  1:175615/275‑1 (MQ=255)                                                                                                                                                                                                 cacaTCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCa                                                                                                                                              >  2:556996/1‑186 (MQ=255)                                                                                                                                                                                                 cacaTCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGAT‑ACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCa                                                                                                                                              >  1:93125/1‑185 (MQ=255)                                                                                                                                                                                                 cacaTCGGCATTCGATGTGACGGAGCGAAACCCTGTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGAT‑ACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGACTTATTTTTCTTCa                                                                                                                                              >  1:456005/1‑185 (MQ=255)                                                                                                                                                                                                     tCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGc                                                                                                                                          <  1:557008/186‑1 (MQ=255)                                                                                                                                                                                                     tCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGAT‑ACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGc                                                                                                                                          <  2:93124/185‑1 (MQ=255)                                                                                                                                                                                                     tCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGAT‑ACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGc                                                                                                                                          <  2:455998/185‑1 (MQ=255)                                                                                                                                                                                                      cGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCg                                                                                                                                                                                 >  1:191739/1‑146 (MQ=255)                                                                                                                                                                                                          aTTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGa                                                                                                                                                                             <  2:191737/146‑1 (MQ=255)                                                                                                                                                                                                           ttagATGTGACGACGCGAAACACTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGAACATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAaca                                           <  2:30852/272‑1 (MQ=255)                                                                                                                                                                                                               aTGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATaa                                                                              <  1:422580/236‑1 (MQ=255)                                                                                                                                                                                                                          cGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGACTTATTTttc                                                                                                                                                  >  1:664070/1‑157 (MQ=255)                                                                                                                                                                                                                              aCCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCa                                                                                                                                              <  2:664056/157‑1 (MQ=255)                                                                                                                                                                                                                                          cTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATgc                                                                                                                                                                                                                      >  1:311665/1‑73 (MQ=255)                                                                                                                                                                                                                                              gtgtATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGcgccg                                                                                                                                                                                                                  <  2:311660/73‑1 (MQ=255)                                                                                                                                                                                                                                                gtATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATAc                                                                                                                                                                                                                                   >  1:553823/1‑54 (MQ=255)                                                                                                                                                                                                                                                    tttttGTAAATCGACGATGATCACCATTGATAACGTCGCTCTCCAAATACGCAc                                                                                                                                                                                                                               <  2:553811/54‑1 (MQ=255)                                                                                                                                                                                                                                                     ttttGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACAAGAAGAGGtt  <  1:721974/274‑1 (MQ=255)                                                                                                                                                                                                                                                         gTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACg                                                                                                                                           >  2:296764/1‑133 (MQ=255)                                                                                                                                                                                                                                                             aTCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAAc                                                                                                                                            <  1:636561/128‑1 (MQ=255)                                                                                                                                                                                                                                                             aTCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCtgt                                                                                                                                       <  1:296769/133‑1 (MQ=255)                                                                                                                                                                                                                                                                        tCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTgg                                                                                                                                                                              >  1:712857/1‑83 (MQ=255)                                                                                                                                                                                                                                                                            cTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGcc                                                                                                                                                                          <  2:712843/83‑1 (MQ=255)                                                                                                                                                                                                                                                                               tGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTcac                                                                                                                                                                                                       >  2:138255/1‑51 (MQ=255)                                                                                                                                                                                                                                                                               tGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCa                                                                                                                                              >  1:242871/1‑108 (MQ=255)                                                                                                                                                                                                                                                                               tGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCa                                                                                                                                              >  1:216928/1‑108 (MQ=255)                                                                                                                                                                                                                                                                                   aaCGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACaaaa                                                                                                                                                                                                   <  1:138256/51‑1 (MQ=255)                                                                                                                                                                                                                                                                                   aaCGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGc                                                                                                                                          <  2:216926/108‑1 (MQ=255)                                                                                                                                                                                                                                                                                   aaCGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGc                                                                                                                                          <  2:242868/108‑1 (MQ=255)                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    GACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAAGTGTATTTTTTGTAAATCGACGATGATCACCTTTGATAACGTCGCGCTGCAAATACGCACTGACCATGCGCCGCTGGATTTCACAAAATAATATCAGGCTCCCTCGTGGAGCCTTTTTTATATCTGCCTTATTTTTCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTT  >  CP009273/1047001‑1047518

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