Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,614,538	T→G	5.1%	intergenic (+277/+143)	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,538	0	T	G	5.1%	294.0 / 4.4	88	intergenic (+277/+143)	yjjV/yjjW	putative DNase/putative pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base T (43/39);  new base G (2/4);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTGCGCCTCGCCATACACGCCGTGGGCATGAAACGCGTTCAGG  >  NC_000913/4614250‑4614822                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                 accatGCCGTAGGCGGGATAA‑GCGTGGCCTGTCGCATGCTGCGGTGGGCGCGCGAGGCCTGATGCGCCGCGGAAGGGGTCGTAGCGTGCCGCCAGGTGTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCGGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCAt                                                                                                                                                                                                                                                                                                  <  1:534736/285‑1 (MQ=11) 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)                                                                                                                 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)                                                                                                                       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:584089/289‑1 (MQ=255)                                                                                                                                              cGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGc                                                                                                                                                   <  2:248152/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)                                                                                                                                                                                         cttaccgggctgccagagatgt‑gCG‑GAACCTTAGGCCGGTTAAGGCGTTCACGCCGCATCGGGCAGTCGGCGCA‑TAATGCGTGATGCTACGCTTGTCGCGTCTTATCAGGCCTACAGGTTTGGGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCAccc                                                                                                          <  1:414805/270‑1 (MQ=18)                                                                                                                                                                                          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)                                                                                                                                                                                                            tgt‑gCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATc                                                                                       >  2:856867/1‑289 (MQ=255)                                                                                                                                                                                                              t‑gCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTCGGCCGTATACGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCa                                                                                      >  2:228386/1‑288 (MQ=255)                                                                                                                                                                                                              t‑gCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGc                                                                                                                                                   <  1:585912/227‑1 (MQ=255)                                                                                                                                                                                                              t‑gCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGc                                                                                                                                                   >  2:585912/1‑227 (MQ=255)                                                                                                                                                                                                                        cGTAGGCCGGATAAGGCGTTCACGCCGCATACGGCGGGCGGCGCA‑GAATGCCTGATGCGAGGCTTGTCGCGTCTTATCATGCCTACAGGTTTGGGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGc                                                                             <  1:132737/289‑1 (MQ=255)                                                                                                                                                                                                                         ttgggCCGTATGAAGGGTTGACGCCGCATCCTGCAGTAGGGGCA‑TAATGCATGATGCGACGCTTTGCGCGTGTGATCATGCCTACAGGTTTGTGCC‑GAACGGTAGGCCGGATAAG‑‑GCGTTCGGGCCGCATCCGGCAGTTGG‑CGGACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCt                                                                            <  1:450637/286‑1 (MQ=255)                                                                                                                                                                                                                          tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCCTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCCCAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCCGCCacaggcc                                                                            >  1:577668/1‑282 (MQ=255)                                                                                                                                                                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCa                                                                          >  2:199994/1‑288 (MQ=255)                                                                                                                                                                                                                                        cGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGc                                                              <  2:147850/288‑1 (MQ=255)                                                                                                                                                                                                                                           tCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATCTCAACCGGCTCACCCCGCGCACCTTTAACGCAACCGCCAACGGCTCAACGGCCTCCGGccgt                                                           >  1:496988/1‑285 (MQ=255)                                                                                                                                                                                                                                            cACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTaa                                                                                             >  2:732047/1‑253 (MQ=255)                                                                                                                                                                                                                                            cACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGgc                                                         >  2:307525/1‑289 (MQ=255)                                                                                                                                                                                                                                            cACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCCTGgc                                                         >  2:646257/1‑289 (MQ=255)                                                                                                                                                                                                                                            cACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGgc                                                         >  2:148945/1‑289 (MQ=255)                                                                                                                                                                                                                                            cACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGCCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCCCCCACCTTTAAAGCATCAGCCAACCGCTCCACGTCTTCCGGCGGGgc                                                         >  2:591182/1‑289 (MQ=255)                                                                                                                                                                                                                                              cgccgcATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGgcgc                                                       <  2:363727/289‑1 (MQ=255)                                                                                                                                                                                                                                                 cgcATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCcgc                                                                                                       <  1:104732/238‑1 (MQ=255)                                                                                                                                                                                                                                                 cgcATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCcgc                                                                                                       >  2:104732/1‑238 (MQ=255)                                                                                                                                                                                                                                                           aGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCttt                                           >  2:9874/1‑288 (MQ=255)                                                                                                                                                                                                                                                             tCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTgc                                         >  2:107551/1‑288 (MQ=255)                                                                                                                                                                                                                                                              cGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTgcg                                        <  1:113170/288‑1 (MQ=255)                                                                                                                                                                                                                                                               ggCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTgcgc                                       >  2:160049/1‑288 (MQ=255)                                                                                                                                                                                                                                                                gcgcA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTGCGcc                                      <  2:496988/288‑1 (MQ=255)                                                                                                                                                                                                                                                                gcgcA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTGCGcc                                      <  1:80543/288‑1 (MQ=255)                                                                                                                                                                                                                                                                           ccgtaggCGACGCTGGTCGCGTCTTATCGTGCCTACAGGTTTGTGCC‑GAACCGTAGGCGTGAGAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCGTCCGGCGTGGCGCTCGCCCAGCTTTGCGCCTCGCCATac                             <  1:670069/281‑1 (MQ=255)                                                                                                                                                                                                                                                                            cTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTGCGCCTCGCCATacac                           >  2:780292/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                tGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTGCGCCTCGCCATACACGCCg                       <  2:86269/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                ggCGACGCTTGTCGCGTCTTTTCTTGCTTACGGGTTTGGGCC‑GAACCGTAGGCGGGATAAG‑‑CCGTTCGCGCCGCATCCGGCAGTTGG‑CGAACAAGGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTGCGCCTCGCCATACACGCCGt                      <  1:229355/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                 gCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTGCGCCTCGCCATACACGCCGTg                     <  1:199994/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                               gTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGACTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTGCGCCTCGCCATACACGCCGTGGGCATGAAACGCGt       <  2:547676/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                                 cTTATCAGGCCTACAGTTGTGTGCC‑GAACCTTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTGCGCCTCGCCATACACGCCGTGGGCATGAAACGCGTTc     <  1:182597/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                                    aTCGTGCCAACAGGTTGGTGCC‑GAACCGTAGGCCGGAGAAG‑‑GCGTTCGCGCCGCATCCGGCGGTTGG‑CGAACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCGGGGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTGCGCCNCGCCATACACGCCGTGGGCATGAAACGCGTTCAgg  <  2:335431/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                 AACGTGCCGTAGGCCGGATAAGGCGT‑TCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTG‑GCGCA‑CAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGT‑GCC‑GAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCA‑TAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCC‑GAACCGTAGGCCGGATAAG‑‑GCGTTCGCGCCGCATCCGGCAGTTGG‑CGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTAACGCATCAGCCAACGGCTCAACGTCTTCCGGCGTGGCGCTCGCCCAGCTTTGCGCCTCGCCATACACGCCGTGGGCATGAAACGCGTTCAGG  >  NC_000913/4614250‑4614822

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: ‑