Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,919,475	G→A	P195L (CCG→CTG)	unknown ←	DNA translocase FtsK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,919,475	0	G	A	97.1%	179.9 / NA	71	P195L (CCG→CTG)	unknown	DNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (42/26);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTCAGGGTTTTCCGGTGCAAAAGCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACAAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCGGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGTTAAACTACCGCCCAGCACGCCTCCCGCAAGGTAAGAGGGCGGGTTATTCAGTGAAAGTGTAGCAA  >  1/1919247‑1919697                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                  aaTTCAGGGTTTTCCGGTGCAAAAGCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCa                                                                                                                                                                                                            <  1:467385/250‑1 (MQ=255)  aTTCAGGGTTTTCCGGTGCAAAAGCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATg                                                                                                                                                                                                          >  1:895330/1‑251 (MQ=255)  aTTCAGGGTTTTCCGGTGCAAAAGCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATg                                                                                                                                                                                                          >  1:610866/1‑251 (MQ=255)  aTTCAGGGTTTTCCGGTGCAAAAGCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCACAATCACGATCTCTTCCGGCGG‑TTCTGCCCTGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATg                                                                                                                                                                                                          >  1:332964/1‑251 (MQ=255)         gTTTTCCGGTGCAAAAGCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTcagc                                                                                                                                                                                                    >  1:350640/1‑250 (MQ=255)           tttCCGGTGCAAAAGCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTcagccag                                                                                                                                                                                                 <  1:758360/251‑1 (MQ=255)           tttCCGGTGCAAAAGCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTcagccag                                                                                                                                                                                                 <  1:320677/251‑1 (MQ=255)               cGGTGCAAAAGCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGt                                                                                                                                                                                             <  1:291139/251‑1 (MQ=255)                    cAAAAGCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGaa                                                                                                                                                                                         <  1:609097/250‑1 (MQ=255)                        aGCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCa                                                                                                                                                                                    >  1:864160/1‑251 (MQ=255)                         gCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTcagcc                                                                                                                                                                                                   >  1:672908/1‑235 (MQ=255)                         gCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCAt                                                                                                                                                                                   >  1:61000/1‑251 (MQ=255)                                       cccATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCCTGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATaa                                                                                                                                                                     >  1:265389/1‑251 (MQ=255)                                        ccATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATaaa                                                                                                                                                                    >  1:161939/1‑251 (MQ=255)                                           ttattCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATaaaaaa                                                                                                                                                                 >  1:626222/1‑251 (MQ=255)                                           ttattCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATaaaaaa                                                                                                                                                                 <  1:850604/251‑1 (MQ=255)                                                cATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGaa                                                                                                                                                                                         <  1:477000/222‑1 (MQ=255)                                                            aaaCATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAgg                                                                                                                                                                                                                        >  1:700157/1‑180 (MQ=255)                                                             aaCATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAAt                                                                                                                                              <  1:585785/251‑1 (MQ=255)                                                             aaCATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAAt                                                                                                                                              >  1:567679/1‑251 (MQ=255)                                                              aCATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATa                                                                                                                                             >  1:599715/1‑251 (MQ=255)                                                               cATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATaa                                                                                                                                            <  1:196273/251‑1 (MQ=255)                                                                  atCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAg                                                                                                                                          >  1:467200/1‑250 (MQ=255)                                                                  atCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGc                                                                                                                                         <  1:495762/251‑1 (MQ=255)                                                                  atCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TACTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGc                                                                                                                                         <  1:488230/251‑1 (MQ=255)                                                                   tCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGc                                                                                                                                         >  1:339230/1‑250 (MQ=255)                                                                    ccTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGccc                                                                                                                                       <  1:165857/251‑1 (MQ=255)                                                                                   aGCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAAt                                                                                                                         >  1:215187/1‑250 (MQ=255)                                                                                    gCCGGAAATATTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATca                                                                                                                       >  1:164166/1‑251 (MQ=255)                                                                                       atgtgtaaaagagccagccGCAGGATTGACTAATCCCGTAAGCTGCTCCCCCGCTTCCTGCCCCGGATTTTCCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTCCCTCTCTTTCAGTTTCGTCGTTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCACCCAGCGGTAGAATTTCATTAATAGACAAAAAAAAAA‑CGCGCCCGAACAGAAAAAAAACCCCACCGCCACAGAAATCACAg                                                                                                                    <  1:132433/234‑1 (MQ=255)                                                                                         cagtatTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACaa                                                                                                                             >  1:626557/4‑240 (MQ=255)                                                                                           atatTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATaaaaaa                                                                                                                                                                 <  1:215778/203‑1 (MQ=255)                                                                                           atatTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAg                                                                                                                >  1:413064/1‑251 (MQ=255)                                                                                            tatTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTcagccag                                                                                                                                                                                                 >  1:160196/1‑170 (MQ=255)                                                                                            tatTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGt                                                                                                               >  1:46537/1‑251 (MQ=255)                                                                                            tatTAATGCGATCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGt                                                                                                               >  1:495189/1‑251 (MQ=255)                                                                                                 aTGCGACCCGCGGAATTGACTATTTCCGTAAGCTCCCCCCCCGTTTCCGGCGCCGGATTTTCCCCAACCCCGATCTC‑TCCGGGGGTTTCTGCCTTGTGTGTTCCCTCTTTTTCATTTTCGTGGTTTTTTT‑CCGGTTTTT‑AAGGGGTTTTTT‑GAATGGCCAGCCAGGGGTAAAATTTCATAAATAGACAAAAAAAAAC‑CGCCCCAGAAAAAAAAAAAAACCCCACCGCCACAAAAATCACAGCCAGTAAGGa                                                                                                          <  1:385443/251‑1 (MQ=255)                                                                                                    cGATCCGCAGGATTGTCTATTCCCGCAGGCTGCTCCCCCTCTCCCTGCGCGGGATTTTCCCCATTCCCGTTCTCTTCGGGCGT‑TTTTGCAGTGTGCTTTTCCTCTTTTTCATTTTCGTCGTTTTTTT‑CCAGTTCTT‑CGGGCGGATTTT‑GCATGGTTAGCCAGGGGTAAATTTTAATAAAAAGACGAAAAAAAAA‑CGCGCCCAAAAAAAAAATAAGGCCCACCGCCACAAAAATCACAGCCAGTAAGGATa                                                                                                        <  1:212200/250‑1 (MQ=255)                                                                                                      aTCCGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGa                                                                                                          >  1:183907/1‑246 (MQ=255)                                                                                                        ccGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGc                                                                                                                                  <  1:253379/220‑1 (MQ=255)                                                                                                        ccGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGcc                                                                                                    >  1:551943/1‑250 (MQ=255)                                                                                                        ccGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCg                                                                                                   >  1:691903/1‑251 (MQ=255)                                                                                                        ccGCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCg                                                                                                   >  1:345686/1‑251 (MQ=255)                                                                                                          gCAGGATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGaa                                                                                                 >  1:229813/1‑251 (MQ=255)                                                                                                             ggATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAAttt                                                                                              >  1:797776/1‑251 (MQ=255)                                                                                                             ggATTGACTAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAAttt                                                                                              >  1:393817/1‑251 (MQ=255)                                                                                                                  gactaaTTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTa                                                                                         >  1:448656/5‑251 (MQ=255)                                                                                                                         tCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGa                                                                                                          >  1:144353/1‑227 (MQ=255)                                                                                                                            ctAAGCTCCTCCCCCTCTCCCTGCGCCGGATTTTCCCCAATCCCGTTCTCTCCCGGCTT‑TTCTGCCTTGTGCTTTTCCTCTTTTTCATTTTCTTCTTTTTTTC‑CCAGTTCTT‑CGGGCGGTTTTT‑GCATGGTCACCCGGCGGTAGAATTTAATAAATAGACAAAAAAAAAC‑CCCGCCCAAACAAAAAAAAAACCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGta                                                                               <  1:610637/250‑1 (MQ=255)                                                                                                                              aaGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGtata                                                                             >  1:362862/1‑251 (MQ=255)                                                                                                                                 ctgccccccccctcccTGCGCCGGTTTTTTCCCAATCACGATCTC‑TCCGGGGGTTTCTCCCATGTGCGTTTCCTCTCTTTCTGTTTCGTCGTTTTTTC‑CCAGGTCTT‑CAGGGGGTTTTT‑GCATGGCCAGCCAGCGGTAAAATTCAATAAAAAGAAAAAAAAAAAA‑CGCCCCAAAAAAAAAAAAAAACCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAAt                                                                          <  1:379995/237‑1 (MQ=255)                                                                                                                                       cccccctccctgcgcgggttttttCCCAATCACTATCTCTCCCGGCGT‑TTCTCCCTTGTGCTTTCCCTCTCTTTCATTTTCGTCGTTTTTTC‑CCAGTTCTT‑CGGGGGGTTTTT‑GCAGGGCCACCCAGGGGTAAAATTTAATAAAAAGACAAAAAAAAAC‑CGCGCCCAAAAAAAAAAAAAACCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAggg                                                                    <  1:756873/231‑1 (MQ=255)                                                                                                                                        cccGCTCCCGGCGCGGGAGTTTCCCCAACCACGATCTCTCCGGGCGG‑TTCTGCAATGTGCGTTTCCTCTCTTTCATTTTCTTGGTTTTTTT‑CCGGGTCTG‑CAGGTGGATTTT‑GCAGGGACAGCCGGCGGAAAAATTTCATAAATAGA‑GAAAAAAAAACCGCGCCCGAAAAAAAAAAAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGt                                                                   <  1:298373/250‑1 (MQ=255)                                                                                                                                        aCCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGt                                                                   >  1:139215/1‑251 (MQ=255)                                                                                                                                            cTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAg                                                                                                                                                                             <  1:190897/142‑1 (MQ=255)                                                                                                                                                     ccGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCG‑CACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGTTAAACTACCGccc                                                      <  1:20115/250‑1 (MQ=255)                                                                                                                                                              tCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGCAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGAATAAATTTAGGGTTAAACTACCGCCCAGCACGcc                                              >  1:849843/1‑250 (MQ=255)                                                                                                                                                                ccAATCACGTTCTTTTCGGGCGT‑TTCTGCCTTGTGCTTTTCCTCTTTTTCAGTTTCGTCTGTTTTTT‑CCAGGTCTA‑GAGGCGGTTTTG‑GCATGGCCAGCCAGCGGAAAAATTTCATAAAAAGAAAAAAAAAAAC‑CGCGCCAAAAAAAAAAAAAACCCCCCCCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGTTAAACTACCGCCCAGCACGCCTcc                                           <  1:22963/251‑1 (MQ=255)                                                                                                                                                                        gATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGTTAAACTACCGCCCAGCACGCCTCCCGCAAGGt                                   >  1:855241/1‑251 (MQ=255)                                                                                                                                                                           ctctTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGcc                                                                                                                  >  1:372414/1‑169 (MQ=255)                                                                                                                                                                             ctTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATaaaa                                                                                                                                                                   <  1:726301/119‑1 (MQ=255)                                                                                                                                                                             ctTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGTTAAACTACCGCCCAGCACGCCTCCCGCAAGGTAagag                              >  1:757398/1‑251 (MQ=255)                                                                                                                                                                               cccGGCGT‑TTCTGCCTTGTGCTTTCCCTCCTTTTCTTTTTCGT‑GGTTTTTT‑CCCGGTTTTCCAGGGGGATTTT‑GCAGGGCCAGCCAGGGGTAAAATTTCATAAAAAAAAAAAAAAAAAC‑CGCGCCAAAAAAAAAAATAACCCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGTTAAACTACCGCCCAGCACGCCTCCCGCAAGGTAAGAggg                            <  1:366745/250‑1 (MQ=255)                                                                                                                                                                                      g‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGTTAAACTACCGCCCAGCACGCCTCCCGCAAGGTAAGAGGGCGGGTTa                     >  1:484874/1‑251 (MQ=255)                                                                                                                                                                                          ctgccttgtgtttttCCTCTCTTTCTGTTTCGT‑GGTTTTTTCCCAGGTCTT‑CAGG‑GGTTTTTGGCAGGGCCAGCCAGGGGAAAATTTTCATAAAAAGACAAAAAAAAAG‑CGCCCCCAAAAAAAAAAAAAACCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGTTAAACTACCGCCCAGCACGCCTCCCGCAAGGTAAGAGGGCGGGTTATTc                  <  1:575736/239‑1 (MQ=255)                                                                                                                                                                                            gCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATa                                                                                                        >  1:83532/1‑163 (MQ=255)                                                                                                                                                                                            gCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACt                                                                                     >  1:502799/1‑182 (MQ=255)                                                                                                                                                                                            gCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGTTAAACTACCGCCCAGCACGCCTCCCGCAAGGTAAGAGGGCGGGTTATTCAg                >  1:844063/1‑251 (MQ=255)                                                                                                                                                                                                 gtgCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGAc                                                                                                                                                                           <  1:661257/92‑1 (MQ=255)                                                                                                                                                                                                          ctctctttaggTTTCGTTGTTTTTTT‑CCAGGTCTG‑CAGGCGGTTTTT‑GCAGGGCCAGCCGGCGGTAAAATTTCATAAAAAGACAAAAAAAAAC‑CGCGCCAAAAAAAAAAAAAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGTTAAACTACCGCCCAGCACGCCTCCCGCAAGGTAAGAGGGCGGGTTATTCAGTGAAAGTGTagcaa  <  1:332027/241‑1 (MQ=255)                                                                                                                                                                                                              cTTTCAGTTTCGTCGGTTTTTT‑CCAGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGTT‑AACTACCGCCCAGCACGCCTCCCGCAAGGTAAGAGGGCGGGTTATTCAGTGAAAGTGTagca   >  1:540682/1‑245 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                  AATTCAGGGTTTTCCGGTGCAAAAGCATTGACCTCGTTCCCATTATTCATCGGCACATAAAACATATCCTGCGCTTTTGATAAGCCGGAAATATTAATGCGATCCGCAGGATTGACAAATTCCGTAAGCTGCTCCACCGCTTCCTGCGCCGGAGTTTTCCCAATCACGATCTCTTCCGGCGG‑TTCTGCCATGTGCGTTTCCTCACTTTCAGTTTCGTCGGTTTTTT‑CCGGGTCTT‑CAGGCGGATTTT‑GCATGGTCAGCCAGCGGTAGAATTTCATTAATAGACGAATAAAAAA‑CGCGCCCGAACAGAAAATAAAGCCCACCGCCACAAAAATCACAGCCAGTAAGGATACGCCGAATTTGCCTAAACTCGGGTATAAATTCAGGGTTAAACTACCGCCCAGCACGCCTCCCGCAAGGTAAGAGGGCGGGTTATTCAGTGAAAGTGTAGCAA  >  1/1919247‑1919697

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