Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,614,448	C→T	5.1%	intergenic (+187/+233)	yjjV → / ← yjjW	putative DNase/putative pyruvate formate lyase activating enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,614,448	0	C	T	5.1%	307.7 / 3.8	94	intergenic (+187/+233)	yjjV/yjjW	putative DNase/putative pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base C (47/41);  new base T (2/3);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTT  >  NC_000913/4614169‑4614728                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                ttCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑Ga                                                                                                                                                                                                                                                                                       >  1:262200/1‑288 (MQ=255)                     gcgcCGGGAAACGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGt                                                                                                                                                                                                                                                                  <  2:635674/289‑1 (MQ=255)                           ggAAGCGGCGGAGGGGATTGGGCAAGCGGTGCGTAATAACACTGATGCGTGGGTTAAAGGGCCGGAGGCCGGATAAGGCGT‑TCACGCCGCATCAGGCAGTTGGCGCACAATGCCGGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                                                                                                                                             <  1:467891/288‑1 (MQ=17)                                      aTGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                                                                                                                                  <  2:788914/288‑1 (MQ=255)                                        gagaTTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt                                                                                                                                                                                                                                               <  2:68635/289‑1 (MQ=255)                                            ttGCGCAAGGGGTGCTTAATAAAACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCGACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCgg                                                                                                                                                                                                                                            <  1:386141/288‑1 (MQ=255)                                                         gCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑g                                                                                                                                                                                                                                                                                        >  1:757746/1‑231 (MQ=255)                                                         gCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑g                                                                                                                                                                                                                                                                                        <  2:757746/231‑1 (MQ=255)                                                         gCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGa                                                                                                                                                                                                                             >  2:565424/1‑289 (MQ=255)                                                            taataaCACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGCCCCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCCTAGGCCGGACAAGGCTTTCACGCCGCATCCGGCCGTCGGCGCA‑CAATGCCTGATg                                                                                                                                                                                                                           >  1:186355/1‑288 (MQ=255)                                                                    cGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTg                                                                                                                                                                                                                  >  2:136247/1‑289 (MQ=255)                                                                    cGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTg                                                                                                                                                                                                                  >  2:739080/1‑289 (MQ=255)                                                                    cGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTg                                                                                                                                                                                                                  >  2:670069/1‑289 (MQ=255)                                                                    cGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGAACAGGCGTTCACTCCGCATCCGGAAGTTG‑GCGCA‑CAATGCCTGATGCGACGCCTGTCCCGTCTTTTCCTGCCTACCAGCATCT‑GCC‑GACCCGTCGGCCGGATTAGGCGCTCCCACCCCTTCCGGCAGTCGgccgttaacacatcgagccacgcg                                                                                                                                                                                                                     >  1:335431/1‑265 (MQ=14)                                                                       ataCGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTc                                                                                                                                                                                                                <  2:764722/288‑1 (MQ=255)                                                                       ataCGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCCACGCTTGTc                                                                                                                                                                                                                >  2:368473/1‑288 (MQ=255)                                                                             tgttTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑g                                                                                                                                                                                                                                                                                        >  1:667525/1‑211 (MQ=255)                                                                             tgttTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑g                                                                                                                                                                                                                                                                                        <  2:667525/211‑1 (MQ=255)                                                                             tgttTAACGCGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCtt                                                                                                                                                                                                         >  2:580527/1‑289 (MQ=255)                                                                               tttAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCGTCGGCGCC‑AAATGCCTGATGCGACCCTTTTCGCGTCTTa                                                                                                                                                                                                        >  1:654603/1‑288 (MQ=255)                                                                               tttAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTa                                                                                                                                                                                                        >  2:132737/1‑288 (MQ=255)                                                                               tttAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTa                                                                                                                                                                                                        >  2:446334/1‑288 (MQ=255)                                                                               tttAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTAt                                                                                                                                                                                                       >  1:345375/1‑289 (MQ=255)                                                                               tttAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTAt                                                                                                                                                                                                       >  1:582018/1‑289 (MQ=255)                                                                               tttAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCCGCGCA‑TAATGCCTTATGCGACGCTTGTCGCGTCTTa                                                                                                                                                                                                        >  2:246014/1‑288 (MQ=255)                                                                                ttAACGTGCGGACGGCCGGATAAGGCGG‑TCACGCGGCGACCGGCAGTTGGCGCACAATGCCTGTTGCGACGCTTAACGCGTCTGATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATc                                                                                                                                                                                                      <  1:616112/289‑1 (MQ=255)                                                                                ttAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATc                                                                                                                                                                                                      <  2:462490/289‑1 (MQ=255)                                                                                ttAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGAACAAGGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATc                                                                                                                                                                                                      <  2:851962/289‑1 (MQ=255)                                                                                ttAACGGGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATc                                                                                                                                                                                                      <  1:565424/289‑1 (MQ=255)                                                                                 tAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGTTTAACGCGTCTTATCATGCCTACAGGTTTGTGAC‑GAACCGTCGGCCGGATAAGGCGTTCACGAAGCATCCGCCAGTTG‑GCGCA‑CAAACCCTGATGCGACGGTTGTCGCGTCTAATCATTCCAGCAAACCTTC‑GCC‑CAACCGTAGCCCGGATACGAGGACGACGCCCCATCCGGCCGCCAGCGCA‑TAAATCCTGAttccgcgcacctctcatcttaccg                                                                                                                                                                                                     >  1:634665/1‑266 (MQ=9)                                                                                  accatGCCGTAGGCGGGATAA‑GCGTGGCCTGTCGCATGCTGCGGTGGGCGCGCGAGGCCTGATGCGCCGCGGAAGGGGTCGTAGCGTGCCGCCAGGTGTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCGGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCAt                                                                                                                                                                                                    <  1:534736/285‑1 (MQ=11)                                                                                  aaCGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCA‑CAATGCCTGATGCGACGCTTGACGCGTCTTATCag                                                                                                                                                                                                    <  2:815808/289‑2 (MQ=255)                                                                                  aaCGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCAt                                                                                                                                                                                                    <  2:262200/289‑1 (MQ=255)                                                                                  aaCGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCAGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCAt                                                                                                                                                                                                    <  2:708439/289‑1 (MQ=255)                                                                                    cGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAAGCCGGATAAGGCGTTCACCCCGCATCCGGCCATCCGCCCA‑TAATGCCTCATGCGACTCTTGTCGCGTCTTATCCTg                                                                                                                                                                                                   >  1:248152/1‑288 (MQ=255)                                                                                       gCCGTAGGCCGGATAAGGCGT‑TCCCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTa                                                                                                                                                                                               <  2:590470/289‑1 (MQ=255)                                                                                        ccGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTATAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTa                                                                                                                                                                                               <  2:42369/288‑1 (MQ=255)                                                                                        ccGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTa                                                                                                                                                                                               <  1:739080/288‑1 (MQ=255)                                                                                         cGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TCATGCCTGATGCGCCGCTTGTCGCGTCTTATCCTTCCTAcc                                                                                                                                                                                             >  2:50988/1‑288 (MQ=255)                                                                                                         cGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAAcc                                                                                                                                                                             <  2:346539/288‑1 (MQ=255)                                                                                                              cACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCAt                                                                                                                                                                                                    <  2:174593/262‑1 (MQ=255)                                                                                                              cACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCAt                                                                                                                                                                                                    >  1:174593/1‑262 (MQ=255)                                                                                                              cACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAg                                                                                                                                                                         >  2:80543/1‑288 (MQ=255)                                                                                                               aCGCCGCATCCGGCAGTTGGCGAACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGc                                                                                                                                                                       <  2:455291/289‑1 (MQ=255)                                                                                                                    gcACCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGAt                                                                                                                                                                  <  2:824639/289‑1 (MQ=255)                                                                                                                                   ggacgaagtcaggaggggaccgttAACGCGGTTTATGATGCATACAGGTTTGTGCC‑GAACCGTAGGCAGGTGAGGGCGTGCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGCTGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTATGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGcgcc                                                                                                                                                  <  1:599549/266‑1 (MQ=17)                                                                                                                                    gcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGcgccg                                                                                                                                                 >  2:609322/1‑288 (MQ=255)                                                                                                                                      acaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCAt                                                                                                                                              <  2:345375/289‑1 (MQ=255)                                                                                                                                           gCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGc                                                                                                                                         <  2:454767/289‑1 (MQ=255)                                                                                                                                                 tGCGACGCTT‑ACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCA‑TCC‑GAGCACCATTGCCTGATGCGACGCTT‑ACCCGTCTTATCAGGCCTACAAACTTGT‑GCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCA‑TCCGAGCACCATTGCCTGATGCGACGCTT‑ACGCGTCTTCTCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCCGtcat                                                                                                                                   >  2:618431/1‑286 (MQ=18)                                                                                                                                                  gCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTgg                                                                                                                                   <  1:136247/288‑1 (MQ=255)                                                                                                                                                    gACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGg‑cg                                                                                                                                <  2:397579/288‑1 (MQ=255)                                                                                                                                                      cGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATCAGGCGTTCCCGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTCTCGCGTCTTCTCATGCCTACCGGTTTGTCCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCCCCGCATCCGGCAGTTGt‑cgccc                                                                                                                             >  1:720365/1‑287 (MQ=255)                                                                                                                                                      cGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGGCTTATCATTCCTACAGGTTCGTGCC‑GACCCCTTGGCCGGACAAG‑‑GCCCTTGCGCCCCAATCTGCCGTTTG‑CGcac                                                                                                                             >  1:147850/1‑289 (MQ=17)                                                                                                                                                      cGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCCTCCGGCAGTTGG‑CGcac                                                                                                                             >  2:575228/1‑289 (MQ=255)                                                                                                                                                       gCTTAACGCGTGTTATCAGGCGTAAAGGTTTGGGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGcaca                                                                                                                            <  1:368473/289‑1 (MQ=255)                                                                                                                                                                gTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGAt                                                                                                                   >  2:182597/1‑289 (MQ=255)                                                                                                                                                                  cTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGc                                                                                                                 <  2:582018/289‑1 (MQ=255)                                                                                                                                                                  cTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCCACGCTTGTCGCGTCTTATCATGCCTACCGGTTTGTGCC‑GAACCGTAGGCCGGATCAG‑‑GCGTTCCGCCCGCATCCGGCAGTTGG‑CCCACAATGCCTGAtt                                                                                                                  >  1:280813/1‑287 (MQ=255)                                                                                                                                                                   ttATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GCACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTCTCGCGTCTTATCATTGCTACATGTTTGTGCC‑CAACCGCACGCCCGATAAG‑‑GCGTCCGCGCCGCCTCCGTGAGCTGGTCGC‑CAATGGCTCATGCg                                                                                                                >  1:547676/1‑289 (MQ=17)                                                                                                                                                                             tACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGAcg                                                                                                      <  2:884762/289‑1 (MQ=255)                                                                                                                                                                                 ggTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GCACCGTAGGCCGGTTAAG‑‑GCGTTCGCGCCGCATCCGGCCGTTGG‑CCCCCAATGCCTGATGCGACGCTTGAcgcc                                                                                                    >  2:242712/1‑286 (MQ=255)                                                                                                                                                                                     tgtgCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCGTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCCGCAGTCGGCGCC‑TACTGCCTCATGCGGCGCTTGTCGCGGCTTATCATGCCTACCGCTTTGTGCC‑CCACCGTCGGCCGGATCCA‑‑GCCTTCGCGCCCCATCCGGCCGTTGC‑CGCACAGTACATGCTGCCAACCCTGCCGCGTCTTa                                                                                               >  1:543081/1‑288 (MQ=14)                                                                                                                                                                                      gtgCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGCTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCCGTTGG‑CGCACAATGACTGATGCGACGCTTTACGCGTCTTAt                                                                                              >  2:712745/1‑288 (MQ=255)                                                                                                                                                                                       tgCC‑GAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTgg                                                                                                                                   <  1:151521/251‑1 (MQ=255)                                                                                                                                                                                       tgCC‑GAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTgg                                                                                                                                   >  2:151521/1‑251 (MQ=255)                                                                                                                                                                                         cc‑GAACGGTAGGCGGGATAAGGCGTGCACGCAGCATCGGGCGGTGG‑GCGCA‑CAGTGCCGGATGAGATGCGTGTAGCGTCTTAGCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAg                                                                                           <  1:543395/288‑1 (MQ=14)                                                                                                                                                                                         cc‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAg                                                                                           >  2:869797/1‑288 (MQ=255)                                                                                                                                                                                          gccatccGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCA‑‑GATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTAC‑GGTTTATGGGC‑GAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTC‑GTGCG‑CCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGC‑GAAGTGTAGGCCTGATAAGACGCGTTAGCGTCGCa                                                                                                                                               >  1:654703/6‑238 (MQ=17)                                                                                                                                                                                          gccatccGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCA‑‑GATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTAC‑GGTTTATGGGC‑GAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTC‑GTGCG‑CCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGC‑GAAGTGTAGGCCTGATAAGACGCGTTAGCGTCGCa                                                                                                                                               <  2:654703/233‑1 (MQ=17)                                                                                                                                                                                             aaCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAgg                                                                                                                                                                                                                                                                               >  1:373540/1‑110 (MQ=40)                                                                                                                                                                                             aaCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGt                                                                                                                                                                                                                 <  2:869200/171‑1 (MQ=255)                                                                                                                                                                                                  tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAg                                                                                  >  2:824857/1‑289 (MQ=255)                                                                                                                                                                                                   aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGt                                                                                 <  2:186355/289‑1 (MQ=255)                                                                                                                                                                                                    ggCCGGTAAAGGCTTTCACGCCGCATCCGGCAGTTG‑GCGCA‑TAATGCCTGATGCGAGGCTTGTCGCGGCTGATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGt                                                                                 <  1:580527/288‑1 (MQ=255)                                                                                                                                                                                                       cGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTAt                                                                                                                                                                                                       >  1:869200/1‑171 (MQ=255)                                                                                                                                                                                                        ggATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCtgtg                                                                            <  2:330755/289‑1 (MQ=255)                                                                                                                                                                                                         gATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGACAAG‑‑GCGTTCGCGCCGCACCCGGCCGTTGG‑CGCCCACTGCCTGCTGCGACGCTTCCCGCGCCTTATCCGGCCCACAAGTCcgtg                                                                            >  1:233971/1‑288 (MQ=14)                                                                                                                                                                                                                         cgcATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTa                                                           >  2:500186/1‑289 (MQ=255)                                                                                                                                                                                                                            aTCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGCTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCCCCCCCTCCGCCAGTTTG‑CGCCCAATACCTGATGCCCCCGTTGACGCATCTTTCCCGGCCCACCAGTCTTCGCCGCACCCTCGGCCGtccc                                                         >  1:517679/1‑285 (MQ=14)                                                                                                                                                                                                                               cGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGc                                                     <  2:248152/289‑1 (MQ=255)                                                                                                                                                                                                                               cGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGc                                                     <  2:584089/289‑1 (MQ=255)                                                                                                                                                                                                                                                     cTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGCTCAG‑‑GCGTTCGCCCCGCACCCGGCAGTTGG‑CGCACAATGCCTGATCCGACGCTTGACCCGTCTTATCAGGCCTACAAGTCCGTGCCGAACCGCCGGCCGTCTCCGGCATGTCCCAACAAGAGCGCCg                                 >  1:86269/1‑289 (MQ=17)                                                                                                                                                                                                                                                      tGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCgg                                <  2:654603/289‑1 (MQ=255)                                                                                                                                                                                                                                                          gCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTTACGCGGCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGGCGTATCCGGCATGTCACCAATAGAGCGCCGGACat                            >  2:793940/1‑289 (MQ=255)                                                                                                                                                                                                                                                          gCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAat                            >  2:113170/1‑289 (MQ=255)                                                                                                                                                                                                                                                              cGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATTCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCCACGCTTGACGCGTCTTATCAGGCCTACCAGTCTGTGCCGCACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGACATAACa                        >  2:401445/1‑289 (MQ=255)                                                                                                                                                                                                                                                              cGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCCCAAATAGAGCGCCGGAAATATCa                        >  2:327199/1‑289 (MQ=255)                                                                                                                                                                                                                                                              cGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCa                        >  2:668115/1‑289 (MQ=255)                                                                                                                                                                                                                                                               gCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCaa                       <  2:280813/289‑1 (MQ=255)                                                                                                                                                                                                                                                                           ttATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCAccc            >  2:732914/1‑288 (MQ=255)                                                                                                                                                                                                                                                                               caggCCTACAGGTTGGT‑GCC‑GACCCGTAGGCCGGATAAGGCGTTCACGCCGCGTCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCGTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCgcgc       <  1:246014/286‑1 (MQ=255)                                                                                                                                                                                                                                                                                   ccTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCt   >  2:229355/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                    cTACAAGGCGGT‑GCC‑GAACCGTAGGGCGGAGAAGGCGGTCACGCCGCAACCGGCAGTCGGCGCA‑GAGTGCCTGATGCGACGCTTGTCGCGTCTTATCAGGCCTACAGGTTTGGGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGAGGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCtt  <  1:575228/289‑1 (MQ=18)                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                TTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTT  >  NC_000913/4614169‑4614728

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