Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,445,274	0	A	C	57.4%	21.9 / 59.0	54	T1108P (ACC→CCC)	unknown	Chromosome partition protein MukB
Reads supporting (aligned to +/- strand):  ref base A (16/7);  new base C (0/31);  total (16/38)
Fisher's exact test for biased strand distribution p-value = 1.16e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-07
Rejected as consensus: Frequency below/above cutoff threshold.

CGAAACACAGCGCGAGCAATTGCGTCAGGTTCAAGCGCAATATGCACAATACAATCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTACGTTCCATGTCGGACAAAGCCCTTGGTGCGTTGCGTACGGCGGTTGCCG  >  1/1445027‑1445514                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  cgAAACACAGCGCGAGCAATTGCGTCAGGTTCAAGCGCAATATGCACAATACAATCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCAccc                                                                                                                                                                                                                                                <  1:846350/251‑1 (MQ=255)           cgcgAGCAATTGCGTCAGGTTCAAGCGCAATATGCACAATACAATCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGaa                                                                                                                                                                                                                                        <  1:84070/249‑1 (MQ=255)           cgcgAGCAATTGCGTCAGGTTCAAGCGCAATATGCACAATACAATCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAaag                                                                                                                                                                                                                                      <  1:125685/251‑1 (MQ=255)                 cAATTGCGTCAGGTTCAAGCGCAATATGCACAATACAATCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAaagcgaag                                                                                                                                                                                                                                 >  1:859198/1‑250 (MQ=255)                                      cAATATGCACAATACAATCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCAccc                                                                                                                                                                                                                                                <  1:724601/214‑1 (MQ=255)                                      cAATATGCACAATACAATCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTg                                                                                                                                                                                                                                          <  1:576057/220‑1 (MQ=255)                                      cAATATGCACAATACAATCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCAt                                                                                                                                                                                                          >  1:761855/1‑251 (MQ=255)                                      cAATATGCACAATACAATCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCAt                                                                                                                                                                                                          >  1:555435/1‑251 (MQ=255)                                              acaatacaatCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGAAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGAAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATGGAAAAACAACTCCCCCTCATGGAAAGCGAAGCGG‑AAAATCTTCCCCGCCGCATTCGTAAAg                                                                                                                                                                                                  <  1:74361/251‑1 (MQ=255)                                                tgttcataattttccggatcgctttGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGAAAAAATCCGTCGGGATGAGTTGTATCAACAACTTC‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCCCCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAg                                                                                                                                                                                                  <  1:228652/226‑1 (MQ=255)                                                 atacaatCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACACCTCCCCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAg                                                                                                                                                                                                  <  1:560059/248‑1 (MQ=255)                                                 atacaatCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCgg                                                                                                                                                                                               >  1:110348/1‑251 (MQ=255)                                                 atacaatCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCgg                                                                                                                                                                                               >  1:506034/1‑251 (MQ=255)                                                      atCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACgg                                                                                                                                                                                          >  1:382498/1‑251 (MQ=255)                                                                 aTTCAATTGCAAAGTTCTTTCAACAATAAAACCCAAATGTGGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGGGGAAGAACGTGCAAAAATCCGGCGGGATGATTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACCCCCCCTCATTGAAAGCGAGGGGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAAc                                                                                                                                                                               <  1:180029/251‑1 (MQ=255)                                                                   taaatggaaaaGTTTTTACAACAGTAAAACCCAAATGTTGCAGAAATTAATGCAGGAAATCGGCAAATTGGTGGTACGTGCCGACGAAGGCGCGGAAGACCGTGCAAAAACCCGTCGGGATGATTTGTATCAACATCTTC‑CCACTAATCTTCAGTGCCTTTCTTAAATTGAAAAACAACTCCCCCTCATTGAAAGCGAAGCGA‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAAccc                                                                                                                                                                             <  1:18172/243‑1 (MQ=255)                                                                              tttCTTTCAACAGAAAAAACAAAATGTTGCAGGAATTAATGCAGAAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGACCGTGCAAAAATCCGTCGGGATGAGTTTTATCAACACCTTT‑CCCCTAATCGTCAGTGCCGTTCTTACATTAAAAAACAACTCCCCCTCATTGAAAGCGAAGCGA‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTg                                                                                                                                                                  <  1:691604/250‑1 (MQ=255)                                                                                   aatgcaaatttgttgttcgcgggtcAGGAATTAAGGCAGAAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAAACCGTGCAAAAATCGGGCGGGATGAGTTATAACAACAACTTT‑CCACTAATCGTCGGTGCGGTTTTTACTTTGAAAACCAACTCCCCCTCATTAAAGGCAAAGCGG‑AAATTCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAAc                                                                                                                                                             <  1:609410/227‑1 (MQ=255)                                                                                   aaggcgaatttgttgttcgccggtcAGGAATAAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGAAAGAACGTGAAAAAACCCGTCGGGATGAGTTGTATAAAAAACTTC‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACCCCCCCTCATGGAAAGCGAAGCGA‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAAc                                                                                                                                                             <  1:520667/227‑1 (MQ=255)                                                                                     tgcgaattttttgttcgccggtcAGGAATTAAGGCAGAAAATCGGCGAATGGGGTGTACGGGCCGACGAGGGCGCGAAAAAACGTGAAAAAACCCGTCGGGATGAGTTGTATCAACATCTTT‑CCACTAATCTTCAGTGCCTTTTTTACATTAAAAAACAACTCCCCCTCATAGAAGGCGAAGCGA‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAAcgg                                                                                                                                                           <  1:45187/229‑1 (MQ=255)                                                                                              aaaCCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTaaaa                                                                                                                                                  >  1:272884/1‑251 (MQ=255)                                                                                                     aTGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACCCCCCCTCATTGAAAGCGAAGCGG‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAAccc                                                                                                                                                                             <  1:663901/217‑1 (MQ=255)                                                                                                           cAGGAATTAAGGCAGGAAATCGGCGAATTGGGTGTCCGTGCCGACGAGGGGGCGAAAGAACGGGCAAAAACCCGTCGGGATGAGTTGTAACAACAACTTC‑CCACAACCCGTCAGTGCCGTTCTTACATTGAAAAACACCCCCCCCTCATTAAAGGCGAAGCGG‑AAATTCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGcgcg                                                                                                                                      <  1:406575/250‑1 (MQ=255)                                                                                                           aagaaataaaTGAAGGAAACCGGGAAATTGGGTGTACGGGCCGACGAAGGCGCGGAAGACCGTGAAAAAACCCGTGGGGATGAGTTGAACCAACAACTTC‑CCACTAACCGTCAGTGCCGTTCTTACATTAAAAAACACCCCCCCCTCATGAAAAGCGAAGCGA‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCgt                                                                                                                                     <  1:826590/243‑1 (MQ=255)                                                                                                               aaTTAATGCAGGAAATCGGGAAATTGGGTGTCCGTGCCGACGAGGGCGCGGAAGAACGTGCAAAAACCCGTCGGGATGATTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCCCCCTCATTGAAAGCGAGGCGG‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCgtggt                                                                                                                                  <  1:197817/250‑1 (MQ=255)                                                                                                                    aTGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGtgc                                                                                                                                >  1:631260/1‑247 (MQ=255)                                                                                                                      gCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGACCGTGCAAAAATCCGTCGGGATGATTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTAAAAAACAACCCCCCCTCATTGAAAGCGAAGCGA‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGt                                                                                                                           <  1:541643/250‑1 (MQ=255)                                                                                                                      gCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGAAAAAACCCGTCGGGATGAGTTGTATCAAAAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTAAAAAACAACCCCCCCTCATGAAAAGCGAAGGGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGt                                                                                                                           <  1:781744/250‑1 (MQ=255)                                                                                                                      gCAGGAAATCGGCGAATTGGGTGTACGGGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATAAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACACCCCCCCCTCATTGAAAGCGAAGCGG‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGt                                                                                                                           <  1:646441/250‑1 (MQ=255)                                                                                                                           aaaTCGGCGAATTGGGTGTCCGTGCCGCCGAAGGCGCGGAAGAACGTGAAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCGGTGCCGTTCTTACATTAAAAAACACCTCCCCCTTATTGAAAGCGAAGCGA‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTAtc                                                                                                                      <  1:412955/250‑1 (MQ=255)                                                                                                                                   aaaTTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGAAAAAATCCGTCGGGATGAGTTGTATCAACAACTTC‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCCCCCTCATGGAAAGCGAAGCGA‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTc                                                                                                             <  1:69552/250‑1 (MQ=255)                                                                                                                                    aaTTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCTTCGGGATGAGTTGTATCAACACTTTT‑CCACTAATCGTCAGTGCCTTTCTTACATTGAAAAACAACTCCCCCTCATTGAAAGCGAAGCGG‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCt                                                                                                            <  1:533662/251‑1 (MQ=255)                                                                                                                                        gggTGTACGTGCCGACGAAGGCGCGGAAGACCGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTAAAAAACAACTCCCCCTCATTGAAAGCGAAGCGA‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACaa                                                                                                                                                                                  <  1:838664/177‑1 (MQ=255)                                                                                                                                          gtgtACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAACCCGTCGGGATGAGTTGTATCAACACCTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACACCTCCCCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACaa                                                                                                                                                                                  <  1:546483/175‑1 (MQ=255)                                                                                                                                             tACGTGCCGACGAAGGCGCGGAAGACCGTGAAAAAATCCGTCGGGATGAGTTGTATCAACAACTTC‑CCACTAATCGTCGGTGCCGTTCTTACATTGAAAAACAACCCCCCCTCATTGAAAGCGAAGCGA‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGaa                                                                                                   <  1:109316/251‑1 (MQ=255)                                                                                                                                                 tGCCGACGAAGGCGCGGAAGAACGTGAAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGA‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGaaaaa                                                                                                <  1:397635/250‑1 (MQ=255)                                                                                                                                                         aaGGCGCGGAAGAACGTGAAAAAATCCGTGGGGATGAGTTGTATCAACAACTTT‑CCACTAATGGTCAGTGCCGTTCTTACATTAAAAAACACCTCCCCCTCATTGAAAGCGAAGCGA‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAAt                                                                                       <  1:663477/251‑1 (MQ=255)                                                                                                                                                                  aagaccgtgaaaaaccccGTCGGGATGAGTTGTATCAACACCTTT‑CCACTAATCGTCAGTGCCTTTCTTACATTAAAAAACAACTCCCCCTCATTGAAAGCGAAGCGA‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGaa                                                                              <  1:120886/235‑1 (MQ=255)                                                                                                                                                                      aCGTGAAAAAACCCGTCGGGATGAGTTGTATCAACACCTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTAAAAAACACCCCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACaa                                                                                                                                                                                  <  1:132017/147‑1 (MQ=255)                                                                                                                                                                        gTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCa                                                                        >  1:643798/1‑251 (MQ=255)                                                                                                                                                                                         gATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCCCCCTCATGAAAAGCGAAGCGG‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGc                                                                                                                                             <  1:253193/165‑1 (MQ=255)                                                                                                                                                                                               gttttacaaaaacTTTCCCACTAATCGTCAGTGCCGTTCTTACATTAAAAAACACCCCCCCCTAATGAAAAGCGAGGCGA‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTa                                                   <  1:91668/240‑1 (MQ=255)                                                                                                                                                                                                 gTATCAACAACTTC‑CCACAAATCGTCAGTGCCGTTCTTACATTAAAAAACAACCCCCCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAAt                                                                                       <  1:870968/211‑1 (MQ=255)                                                                                                                                                                                                   aTCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTACGTTc                                              >  1:885206/1‑250 (MQ=255)                                                                                                                                                                                                       accatcttcccattAATCGTCGGGGCCTTTTTTACATTAAAAAACAACCCCCCCTCATAAAAAGCGAAGCGA‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTTGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTACGTTCCATg                                          <  1:472357/238‑1 (MQ=255)                                                                                                                                                                                                        tgtgtaaaagagccagcGTCAGTGCCTTTCTTAAATGAAAAAACACCCCCCCCTCATTGAAAGCAAAGCGA‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTACGTTCCATGt                                         <  1:358667/235‑1 (MQ=255)                                                                                                                                                                                                        acaccttcccccTAACCGTCAGTCCCTTTCTTAAATTAAAAAACACCCCCCCCTCATAAAAAGC‑AAGGGGAAAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTACGTTCCATGt                                         <  1:88699/240‑1 (MQ=255)                                                                                                                                                                                                                        cGTCAGTGCCGTTCTTACATTGAAAAACACCTCCCCCTCATTAAAAGCGAAGCGG‑AAAATCTTCCCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGAtg                                                                                                        <  1:461875/172‑1 (MQ=255)                                                                                                                                                                                                                          tCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTACGTTCCATGTCGGACAAAGCCCTTGGtg                       >  1:478955/1‑251 (MQ=255)                                                                                                                                                                                                                             gTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTACGTTCCATGTCGGACAAAGCCCTTGGtgcgt                    >  1:893588/1‑251 (MQ=255)                                                                                                                                                                                                                               gCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTACGTTCCATGTCGGACAAAGCCCTTGGtgcgttg                  >  1:111261/1‑251 (MQ=255)                                                                                                                                                                                                                                  gTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTACGTTCCATGTCGGACAAAGCCCTTGGtgcgttgcgt               >  1:793042/1‑251 (MQ=255)                                                                                                                                                                                                                                     cTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTACGTTCCATGTCGGACAAAGCCCTTGGTGCGTTGCGTAcg            >  1:278927/1‑251 (MQ=255)                                                                                                                                                                                                                                                aaaCAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTACGTTCCATGTCGGACAAAGCCCTTGGTGCGTTGCGTACGGCGGTTTCCg  >  1:250410/1‑250 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  CGAAACACAGCGCGAGCAATTGCGTCAGGTTCAAGCGCAATATGCACAATACAATCAAGTGTTGATTCAATTGCAAAGTTCTTTCAACAGTAAAAACCAAATGTTGCAGGAATTAATGCAGGAAATCGGCGAATTGGGTGTACGTGCCGACGAAGGCGCGGAAGAACGTGCAAAAATCCGTCGGGATGAGTTGTATCAACAACTTT‑CCACTAATCGTCAGTGCCGTTCTTACATTGAAAAACAACTCACCCTCATTGAAAGCGAAGCGG‑AAAATCTTACCCGCCGCATTCGTAAAGCGGAACGGGACTACAAAACCCAGCGCGAATTGGTAACGGCGGCTAAAATGAGCTGGTGCGTGGTGCTGCGTTTATCTCGCAATTCTGATGTGGAAAAACGCTTAAATCGTCGTGAATTGGCATATTTATCCGCCGATGAGCTACGTTCCATGTCGGACAAAGCCCTTGGTGCGTTGCGTACGGCGGTTGCCG  >  1/1445027‑1445514

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