Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,571	A→G	E323E (GAA→GAG)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,571	0	A	G	98.7%	247.8 / ‑4.7	79	E323E (GAA→GAG)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (34/44);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAA‑AAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCC  >  NC_000913/4620325‑4620817                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       cgTGCTGCAGAAACTGGTTGTTGACACACAAGGTCACGTGGTTTATGTCGGTAAAATTGCCGGCATCTACGCAAACTGCGGTATCACAAGG‑AAAGCGAAAGCAACTGGCCTGGACGCGTTGTTTGACGCCACCATCAAAGAGATGAAAGAGGCGGGTGCTAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTCTGCTGCGGGTATGGAGCTGt                                                                                                                                                                                                                                                    >  2:301998/1‑251 (MQ=255)    gctgaAGAAAATGGTTGAAGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATGGCAGACATCGACGCCAACTGCGGTAACACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTGCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCg                                                                                                                                                                                                                                                 <  2:373389/251‑1 (MQ=255)      tgcAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                                                                                                                                                                                                                                               >  1:320015/1‑251 (MQ=255)         aGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAccgc                                                                                                                                                                                                                                            >  1:492654/1‑251 (MQ=255)         aGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAccgc                                                                                                                                                                                                                                            >  1:370290/1‑251 (MQ=255)                       gaAAAACTCGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGGGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCt                                                                                                                                                                                                                              <  1:459015/251‑1 (MQ=255)                       gaAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGc                                                                                                                                                                                                                               <  2:307761/250‑1 (MQ=255)                         aaaaCACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCtt                                                                                                                                                                                                                             >  2:457614/1‑249 (MQ=255)                           aacacgttcAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATg                                                                                                                                                                                                                          <  1:432244/243‑1 (MQ=255)                            acacGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCtgttgt                                                                                                                                                                                                                         >  2:157239/1‑251 (MQ=255)                                     ggtggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCtgctgc                                                                                                                                                                                                                >  1:351695/1‑251 (MQ=255)                                     ggtggatcatgacGGTAAAATTGCAGACATCTCCGACAACTGCGGTATAACCAAG‑ACAGACAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCCAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTC‑CCACACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTTTGGGGCTGTTCGACCGCCGTCTGCCGGCGCTGATGTCTCtgctgc                                                                                                                                                                                                                >  2:493516/13‑251 (MQ=255)                                              tGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGGGCTGTTCGACCGCCTTCTGCCGTGGCTGATGTCTCTTCTGCGCGagtacg                                                                                                                                                                                                       >  2:348502/1‑246 (MQ=255)                                                 cGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTctgct                                                                                                                                                                                                                  >  1:361863/1‑237 (MQ=255)                                                 cGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTctgct                                                                                                                                                                                                                  <  2:361863/237‑1 (MQ=255)                                                  ggTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAt                                                                                                                                                                                                   >  1:319837/1‑251 (MQ=255)                                                  ggTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAt                                                                                                                                                                                                   >  2:406045/1‑251 (MQ=255)                                                      aaaTTGCAGACATCTACGCCAACTGCGGTATCACCACG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCCAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTAGGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTg                                                                                                                                                                                               >  2:73568/1‑251 (MQ=255)                                                           gcagACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCt                                                                                                                                                                                          <  1:87975/248‑1 (MQ=255)                                                              gACAGCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACCCCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCAc                                                                                                                                                                                       <  2:319837/251‑1 (MQ=255)                                                                   cTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCActgctg                                                                                                                                                                                  >  1:373871/1‑251 (MQ=255)                                                                     aCGCCAACTGCGGTATCACCAAA‑AAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGAc                                                                                                                                                                                >  1:145314/1‑251 (MQ=255)                                                                        ccAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCCCTGCTGAcccc                                                                                                                                                                             >  2:390609/1‑249 (MQ=255)                                                                         cAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACCTCCTGCTCCTCACTGCTGACCACg                                                                                                                                                                            >  2:401039/1‑251 (MQ=255)                                                                            acgtcgttataCCAAAGAAAGTTAAAGCAACGTGCCTGGACGCGCTTTTTGACGCCACCATCAAAGAGAGGAAAGAGGCGGGTGATATCACCATCGTCTTCTCCA‑ACTGCTTTGACTTCGACGCTTCCTGGGGCCACCGGCGCGCCGTCGCCGGTTATGCTGCGGGTCTGGGGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGt                                                                                                                                                                          <  2:268854/241‑1 (MQ=255)                                                                                     cACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCg                                                                                                                                                                 <  2:373871/250‑1 (MQ=255)                                                                                      aCCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGAc                                                                                                                                                               <  1:73568/251‑1 (MQ=255)                                                                                                   aaaGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTAc                                                                                                                                                   >  1:295513/1‑251 (MQ=255)                                                                                                       caaCTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGAc                                                                                                                                               <  1:258683/249‑1 (MQ=255)                                                                                                       caaCTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGAc                                                                                                                                               >  1:79757/3‑251 (MQ=255)                                                                                                        aaCTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGAcc                                                                                                                                              >  1:461848/2‑251 (MQ=255)                                                                                                         aCTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACca                                                                                                                                             <  1:155146/251‑1 (MQ=255)                                                                                                         aCTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACca                                                                                                                                             <  1:134740/251‑1 (MQ=255)                                                                                                         aCTGGCCTGGACGCGCTGTTGGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACGTCGGTGACTTCGACTCTTCCGGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACca                                                                                                                                             <  2:86622/251‑1 (MQ=255)                                                                                                            ggCCTGGACGCGCGGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACcaca                                                                                                                                           <  1:425158/250‑1 (MQ=255)                                                                                                            ggCCGGGACGCGCTGTTTGACGCCCCCATCAAAGAGATGAAAGAAGCGGGTTATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTGCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACcac                                                                                                                                            <  2:351695/249‑1 (MQ=255)                                                                                                              ccgggACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCGTCACCA‑ACGTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGGTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACAcgc                                                                                                                                        <  1:154931/248‑1 (MQ=255)                                                                                                               cTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACCcgcg                                                                                                                                       >  2:400969/1‑251 (MQ=255)                                                                                                                        cTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACAt                                                                                                                              >  1:291553/1‑251 (MQ=255)                                                                                                                               ccGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTCCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTa                                                                                                                       <  2:130684/250‑1 (MQ=255)                                                                                                                               aCGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGt                                                                                                                        >  1:177868/1‑250 (MQ=255)                                                                                                                               aCGCCACCAGCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCTCCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGTCCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTTATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTa                                                                                                                       <  1:292158/251‑1 (MQ=255)                                                                                                                                 gCCACCATCAAAGAGATGAAAGAATCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACt                                                                                                                     <  1:457614/251‑1 (MQ=255)                                                                                                                                    accaTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGt                                                                                                                  <  1:406045/251‑1 (MQ=255)                                                                                                                                         cAAAGAGATGAAAGAAGCGGGTGATAACACCATCGGCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCGGGGGCCACCGGCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGtttacg                                                                                                             <  2:111492/251‑6 (MQ=255)                                                                                                                                         aaaaGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGtttac                                                                                                              <  2:241862/249‑5 (MQ=255)                                                                                                                                            agagATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGcc                                                                                                          <  1:400969/251‑1 (MQ=255)                                                                                                                                              agaTGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGGCTCTCCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCg                                                                                                        <  2:174773/251‑1 (MQ=255)                                                                                                                                                 ggAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACGTCGTTGACTGCGACGCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGaaa                                                                                                     <  2:461848/250‑1 (MQ=255)                                                                                                                                                        agCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCGTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGGCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAAc                                                                                              <  1:316473/251‑1 (MQ=255)                                                                                                                                                         gCGGGTGATAACACCATCGTCTCCACCA‑ACTTCGGTGACTGCGACTCTTCCGGGGGCCACCGCCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGCCCACTGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAAcc                                                                                             <  2:33991/251‑1 (MQ=255)                                                                                                                                                             gTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGc                                                                                         <  1:113761/251‑1 (MQ=255)                                                                                                                                                                    cacCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGCGCTGTTCGACCGCCGTCTGCCGGAGCTGTTGTCTCTGCTGCGCGCTGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGGACTGGTTTACGGCCAGaacagaaaacccgggacactgg                                                                                  >  2:16292/1‑231 (MQ=255)                                                                                                                                                                     acCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTggg                                                                                 <  1:98955/251‑1 (MQ=255)                                                                                                                                                                       cATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTTGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTAACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGAAAAACCGGGCTCACTGGGTc                                                                               >  1:46373/1‑251 (MQ=255)                                                                                                                                                                           gTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCaccg                                                                           >  2:425363/1‑248 (MQ=255)                                                                                                                                                                              ttCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGCTCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGGAAAACCGGGCTCACTGGGTCaccggg                                                                         >  1:293517/1‑245 (MQ=255)                                                                                                                                                                               tCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTg                                                                         <  2:492654/249‑1 (MQ=255)                                                                                                                                                                                cACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGGTCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGaa                                                                       <  1:441806/250‑1 (MQ=255)                                                                                                                                                                                cACCA‑ACTTCGTTGACTTCGACTCTTACTGGGGCCACCGTCACGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCaccgcgaa                                                                       >  1:388355/1‑243 (MQ=255)                                                                                                                                                                                  ccA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAAcc                                                                    >  1:241578/1‑251 (MQ=255)                                                                                                                                                                                          cGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTACAACCGGGCTCACTGGGTCACCGTGAAACCTTCCCgg                                                             >  2:325779/1‑251 (MQ=255)                                                                                                                                                                                           gTTGACTTCGCCTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTCTGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGGTGTGCTCCGACCTGGACCGGTACTGACCACACGCGTGAACACCTTCCGGTACTGGTTTACGGCCCGAAAGTCAAACCGGGCTCACTGGGTCACCAGAAAACCTTCGCggg                                                            >  2:203636/1‑250 (MQ=255)                                                                                                                                                                                            ttGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGAAAAACCGGGCACCCTGGGGCACCGTGAAACCTTCGCgg                                                             >  2:269262/1‑249 (MQ=255)                                                                                                                                                                                              gACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGatat                                                         >  1:408832/1‑251 (MQ=255)                                                                                                                                                                                                          tCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGAct                                               <  1:157239/249‑1 (MQ=255)                                                                                                                                                                                                          tCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGAc                                                <  1:192085/248‑1 (MQ=255)                                                                                                                                                                                                          tCCTGGGGCCACCGGCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGAct                                               <  1:348502/249‑1 (MQ=255)                                                                                                                                                                                                           ccTGGGGCCCCCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTg                                            <  1:67650/251‑1 (MQ=255)                                                                                                                                                                                                                    acagTCGCGACGTCGCCGTGTATGCGGCGGTTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCTCACGCGTGAACACATTCCGGTACTTGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCGGGCAAAatat                                   <  2:67907/248‑1 (MQ=255)                                                                                                                                                                                                                     ccGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGTTCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGTAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATAtt                                  <  1:220221/251‑1 (MQ=255)                                                                                                                                                                                                                      cGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATAttt                                 <  1:183883/251‑1 (MQ=255)                                                                                                                                                                                                                          gcgACGGCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGt                             <  1:314170/251‑1 (MQ=255)                                                                                                                                                                                                                           ggACGTCGCCGGTTATGTTGCGGGTCTGGAGCTGGTCGACCGCCGTCTGCCGGAGCTGATGTCGCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGGTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTa                            <  2:370290/250‑1 (MQ=255)                                                                                                                                                                                                                           cgACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTa                            <  1:305268/251‑1 (MQ=255)                                                                                                                                                                                                                                 cGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTTACCACGGTTGTGATCCGACCTGGACCGGTACGGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTg                      <  2:447211/251‑1 (MQ=255)                                                                                                                                                                                                                                 cGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGCAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTg                      >  1:474111/1‑251 (MQ=255)                                                                                                                                                                                                                                         aTGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGaa              <  1:425363/251‑1 (MQ=255)                                                                                                                                                                                                                                                     tGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGcc  >  1:215990/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAA‑AAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCC  >  NC_000913/4620325‑4620817

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