Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,967,843	A→C	17.8%	V197G (GTG→GGG)	ptsP ←	PEP‑protein phosphotransferase enzyme I, GAF domain containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,967,843	0	A	C	17.8%	142.9 / 4.0	79	V197G (GTG→GGG)	ptsP	PEP‑protein phosphotransferase enzyme I, GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base A (37/21);  new base C (5/13);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCG  >  NC_000913/2967609‑2968078                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              cctgTTTTTCCCCCCCACCCTGCCCCCG‑GGCTTTAAACACCTCCGAAAAAAATTCCCCGACGGAGCGGGGTTTCCGAAGGCGGGTGGGAGAAAA‑GACCAAAAAGAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATacac                                                                                                                                                                                                                                              <  1:98722/233‑1 (MQ=255) aCCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTGGCGTTCCAGG‑GCCGGATCCAGCCT‑TGATGCCTGTTAcccccg                                                                                                                                                                                                                                           >  2:1293154/1‑237 (MQ=255)   cttttttcccccccccttcccccccGGGCCTTACTCACCCTGGGCAAACATTT‑CGGGGCCGCGCCGGGTTGCCGAAAGCCGGGGCGGAGAAA‑GCTCGAAAATTGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑Gtt                                                                                                                                                                                                                                        <  1:538403/216‑1 (MQ=255)             gCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTTGTTTGCCGCTTCTTCCACGGCCCCGGTCAA‑TCGTTCGCGGTCCTAG‑GCCGGCTCCACGCT‑TGTGCCCTGATAAAACC‑GTCCCATTaaaag                                                                                                                                                                                                                              >  1:1281508/1‑237 (MQ=255)             gCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCCTTTTCCGCTTCTTCCACCGCCCCGGGCAG‑TCGTTCGCGTTCCAGC‑GCCGGATCCCACGT‑GGATGCCCGATACAACT‑GGTCCATTAAAgg                                                                                                                                                                                                                              >  1:1062448/1‑239 (MQ=255)                             tctttttATCAACCCGCGAAAAAAATT‑CCGGGCGGCGCCGGGTGTCCGAAAGGAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGTTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGc                                                                                                                                                                                                              <  1:1293154/235‑1 (MQ=255)                             aGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGACAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGCAACCCTTTTCCCCCTACTTCCCGCGCCCCCGCCAG‑TCGTCCGCGTTCCCAC‑GCCCGGTCCCCTCT‑TGAGCCCTGATAAACCC‑GTTCCCTTaacgtgaccctgccaccctc                                                                                                                                                                                                               >  1:130456/1‑220 (MQ=255)                             aGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGCAACTCGTTTTCCGCTTCTTCCAGCGCCCCGGGCGG‑TCTCTCCCGTTCCCGA‑CCCGGCGCCA‑CGTGTGAAGCCTGGCACCTCT‑GCTCCCTTAAAAGGAAAGGGGcctcctt                                                                                                                                                                                                               >  1:549934/1‑237 (MQ=255)                             aGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCCTTAAAGGGAACGTGGCATCCTGc                                                                                                                                                                                                              >  2:294342/1‑239 (MQ=255)                             aGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTg                                                                                                                                                                                                               >  2:436198/1‑238 (MQ=255)                             aGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAACGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATAAACCT‑GTTCCATTAAAGGTAACGTGGCACCCTg                                                                                                                                                                                                               >  2:56631/1‑238 (MQ=255)                                    aTCAACCTGGGCAAACAATC‑CCGGACCAAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCtt                                                                                                                                                                                                       <  2:683396/239‑1 (MQ=255)                                        ccccgggcaacaacgt‑CGCGGCCCGGCGGGGTGGCCGAAAGCCGG‑GGGAATAAAGGACCAAAAATGCCCGCCGTTTCGTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTGGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGc                                                                                                                                                                                                   <  1:920164/224‑1 (MQ=255)                                          cTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGGAACGTGGCATCCTGCCAGCCTTCGGCa                                                                                                                                                                                                  >  2:464027/1‑238 (MQ=255)                                            cGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGCAACTCGTTTGCCGCTTCTTCCATCGCCCCGGTCAA‑GCGGTCACGTTCCAGA‑GCCGGGTCCAGCCT‑TTATGCCTGATAAACCT‑GGTCCATTAAAGACACCGAGGCAACCGGCCACCCTTCGACCATc                                                                                                                                                                                               >  1:406947/1‑239 (MQ=255)                                             ggCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAACCTCTTTTTCCGCTTCTTCCCGCGCCCCGGTCAG‑GCGTTCGCGTTCCAGA‑TCCGGATCCACCGT‑TGATGCCTGATACCCCC‑TTTCCATTCAAGGTACCGCGCCACCCTCCCAccccctggccacg                                                                                                                                                                                               >  1:85209/1‑228 (MQ=255)                                             ggCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCCGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATc                                                                                                                                                                                               >  2:53215/1‑238 (MQ=255)                                             ggCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATc                                                                                                                                                                                               >  2:217390/1‑238 (MQ=255)                                                                         gtCCGAAAGCGGGGGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCCCCGGCGGCAAAGCGTTTGCTGTAGCTGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATgcg                                                                                                                                                                    <  2:510801/238‑1 (MQ=255)                                                                          tCCGAAAGACGGTGCGAGTAAA‑GATCGAAAATAACCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGCAACTCGTGGTCCGTTTCTTACAGCGCCCCGGTCAG‑TCGTTCGCGAGCCAGG‑ACCGGATCCACCAG‑TTATGCCGTATACACCC‑TTTCCCTTAAAGGCCACGTGGAATCCTGCACACCTTCGGCAATCGCCACACCacgagcagcacgaaacgggcg                                                                                                                                                                  >  1:1059782/1‑219 (MQ=255)                                                                              actggggggctgctataa‑GGTCGAAAATTCGCGCCTTTTTTTTTGGTGCCCCGGGGGCAAAGGCTTTGCTGGTGGGGCGAACCCCGTTTCCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGGTCGCGTTCCAGG‑GCCGGATCCAGCGT‑TGATGCCTGATACCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGAtt                                                                                                                                                              <  1:1205989/223‑1 (MQ=255)                                                                                   gggggCGAGTAAG‑AATCGAAAATGGCCGCCGTTTCTTTTTGGGCCCCGGCGGCAAAGCCTTTGCTGGAGCGGCGAAACTGGTTTGCCGGTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                                                                                                                          <  1:536349/234‑1 (MQ=255)                                                                                    ggTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCC‑GTTCCCTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGCCAATCCCCACACCCGGTGCTCCCGGTAATTGCCGGATTctcg                                                                                                                                                          >  2:407001/1‑237 (MQ=255)                                                                                     gTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTcgcg                                                                                                                                                          >  2:177079/1‑236 (MQ=255)                                                                                            tAAA‑GTCCGAAAATCGCCCCCGTTTTTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGat                                                                                                                                                 <  1:294342/238‑1 (MQ=255)                                                                                               a‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGGGTTTGCTGTAGCGGCGAAACTGGTTTGCCCCTTCTTCCAGCGCCCCGGTCAG‑GCGTTCGCGGTTCAGG‑GCCCGATCCATCGT‑TGATGCCTGGTACACCT‑GTTACATTCAAGGGAACGTGGAATCCAGCCGGCCTTTGCCAATACCCCCCCCCGGGGCACCCGCTAAAGCGCGGAGTcgcccaaagcgacaat                                                                                                                                              >  1:209410/1‑225 (MQ=255)                                                                                                   gcaaaatttgcggcctttgttttttGGTCCCCGGCGGCAAAGCGTTTGTCTTGGCGCCGATACTCGTTTGCCGGTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                                                                                                                            <  1:746719/217‑1 (MQ=255)                                                                                                   ccGAAAATAGCCGCCGTTTCTTGTTGGGCACCGGCGGCAAAGCGTTGGCTGTAGCGGCGAAACTCGGTTGCCGCTGCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGGAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGc                                                                                                                                            <  2:85209/236‑1 (MQ=255)                                                                                                    cGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGGTCCAACGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGGGGCATCCTGCCAACCTTCGGCAACCCCCCCACCCGGCGATGCGGGTACTAGGGGGATTcggccctggggaacatgccca                                                                                                                                         >  1:493240/1‑220 (MQ=255)                                                                                                     gAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCCGCGAAACTCGTTTGTCGCTTCTTCCACGGCCCCGGTCAG‑TCGTTTCGGTTCCAGA‑GCCGGGTCCCGCGT‑TGCTGCCTAGTACACCA‑TTTTCCTTAAAGGTAACCGGGAAGACTGCAGGCCGTCGGCACACGCCCACCCAGGGCATACCGGAGACGCCGGGATTCGCGTCTgggagaatggcccaa                                                                                                                                        >  1:654057/1‑226 (MQ=255)                                                                                                           aGCCCCCGTTTCTTTTTGGGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAgg                                                                                                                                  <  2:1062448/239‑1 (MQ=255)                                                                                                                             tgCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCCGCGC‑TGATGCCTGAAACCCCT‑GTTCCATTAAAGGGAACGTGGCATCCTCCCAGCCTTAGGCAATCGCCACACCAGGTGCTGCCGGTAATACGCGGATTCGCCCATGGCGATATTGACCAAACaaagaaagcaaagggaaacaaa                                                                                                                >  1:159466/1‑219 (MQ=255)                                                                                                                              gcgcgggggggaaaGCGTTTGCTGTAGCGGCGAACCCCGTTTCCCTCTTCTCCCAGCGCCCCGTCCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCGGCGT‑TGAGGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCGGCCAGCCTTCGGCAATCGCCACGCCAGGTGCGGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTgaga                                                                                                              <  1:56631/228‑1 (MQ=255)                                                                                                                              gccccGGGGGCAAAGCGTTTGTTGTGGGGGCCAAATCCTTTTTCCGTTTCTTCCAGCGCCCCGGTGAG‑TCGGTCGGGTGCAAGA‑GCCGGATCCAGGGT‑TGATGCCTGATACACCA‑GTTCCATTAAAGGTACCGTGGCATCCTGCCAGCCTTCGGCAATTGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTgaga                                                                                                              <  1:1139700/236‑1 (MQ=255)                                                                                                                               cACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGGTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCCTCCTTCGGCCATCCCCCCACCAGGGGTGTCCGGTAATGCGGCGGTTCGCGTCTGGCGCAATTGCCCAAACAAgcaccgcaaccgcgactcaaaa                                                                                                              >  1:742631/1‑218 (MQ=255)                                                                                                                                         aaaGCGTTTGCTGGAGCGGCCAATCTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTGCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATCCCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCGTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCt                                                                                                    <  2:654057/238‑1 (MQ=255)                                                                                                                                              gtttgtgttcgcgacgAACATCTTTTGCCGCTTTTTCCGGCCCCCCGGTCAG‑CCCTTCGTGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACGCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTGCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATc                                                                                              <  1:562472/229‑1 (MQ=255)                                                                                                                                               tttGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCCGCGCCCCGGTCAG‑TCGTTCGAGGTCCCGA‑GCCGCATCCAGCGT‑TGCTGCCTGATACCCCT‑GTTCCATTAAAGGCAACGTGGCATCCTGCCAGCCTTCCGCCATCGCCACACCAAGAGCTGACGGTCATGCCGGGAATCTCGGCTGCCGATATTGCCCCACCAACGCAGACAACTGCT‑ACGGagaaagaaaaacagcat                                                                                                >  1:5396/1‑227 (MQ=255)                                                                                                                                               tttGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCCG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGTGATTGCCCAAACAAGCGAGTCAACTGCG‑ACTGAGAAAGAATAACTGCCa                                                                                                >  2:382798/1‑236 (MQ=255)                                                                                                                                               tttGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGGAAAGGGGCATCCTCCCAGCCTTCCGCAATCACCACACCCGGTGTTGCAGGCAATGAGGGGATTCGGGGCTGGCGATATTGCCCAACCCAGGCAGTC‑ACCGCGCACTGagaaagaaaaacggccc                                                                                                >  1:939845/1‑227 (MQ=255)                                                                                                                                                  gCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATCCACCT‑GTTCCATTAAAGGTAACGAGGCATCCTGCCAGCCCTCGGCAATAGCCACACCACGTGCTGCCGGTAATGCGCGGATTCGCGCCCGCCGTAATTTCCCAAACCAAGAGCTAAACTACG‑CCTGcgaaagagacgaatcccactcgg                                                                                          >  1:188223/1‑223 (MQ=255)                                                                                                                                                       aGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGGGGCa                                                                                     >  2:1219787/1‑239 (MQ=255)                                                                                                                                                          cgggacaccgcgttGCCTTTTCCTCCCGCGCCCCGGTCGG‑CTGTTCGCGTTCCAGA‑GACGGTACGAGCGT‑GTATCGCTGATACGACC‑GTTCCAAGAAAGGTATCGGGGCATGCTGCCGGCCTTCGGAAATCGACACACCAGGTGCTGCCGGGAATGAGCGGATTCGCGTCTGGCGCTATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCaaac                                                                                  <  1:1303039/227‑3 (MQ=255)                                                                                                                                                            cgAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCGGCAACCTGGGTGGCAAGGGt                                                                                >  2:134949/1‑239 (MQ=255)                                                                                                                                                                   cGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GGCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCagg                                                                         <  2:380559/239‑1 (MQ=255)                                                                                                                                                                     ttttCCGTTTCTCCCGGCGCCCCGGTCCG‑TCTGTCGCGTTCCGGA‑GCCGGATCCGGCGT‑TGAGGCCTGATCCCCCT‑GTTCAATTAAAGGTAACGTGGCAGCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaa                                                                       <  1:924538/239‑1 (MQ=255)                                                                                                                                                                     tttGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaa                                                                       <  2:1209132/239‑1 (MQ=255)                                                                                                                                                                      ttGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGCGAAAGACAAGCTGCCACCTGGGGGGCAGGCGCAACCaggaag                                                                      >  2:1104803/1‑239 (MQ=255)                                                                                                                                                                      ttGCCGCTTCTTCCAGCGCCCCGGTCAA‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaag                                                                      >  2:728491/1‑239 (MQ=255)                                                                                                                                                                        gCCCCTTCTTCCCCCACCCCGGTGAC‑CCCTTCGCGGTCCAGG‑CGCGGACCCCGAGT‑TTGTGCCTGA‑ACCCCTGGTTCGTTTAAAGGTAACGGGGCACTCTCCAAGCTTGCGGAAATCGGCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCCACGGGCGATATTGCCCAAACAAGGCATTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGTTGGCAAGCGTCACCaggaagga                                                                    <  1:1161628/239‑1 (MQ=255)                                                                                                                                                                              tcttcCAGCGCCCCGGTCAG‑TCGTTTGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGACGTCAACTCCG‑ACTGAGAAAGAGTAGCTTCCATCTGGGTGGCAACGGTGACCaggaaggtgtctc                                                               >  2:533069/1‑232 (MQ=255)                                                                                                                                                                                tcccAGCCCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCAc                                                            <  2:1112422/237‑1 (MQ=255)                                                                                                                                                                                  ccAGCGCCCCGGGCAG‑TCGTTCGCGTTCAAGG‑GCGGGATCCCAGGT‑TGAGGCCTGATACACCT‑GTCCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACtt                                                          <  1:464027/239‑1 (MQ=255)                                                                                                                                                                                   cAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACttt                                                         <  2:159466/239‑1 (MQ=255)                                                                                                                                                                                     gcgcCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACGCCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCg                                                       <  2:130456/239‑1 (MQ=255)                                                                                                                                                                                      ggcCCCCGGCAGTTCGGTCGGGTCCAAGA‑GCCGGGTCCGGTGT‑TGATGCCTGAAACACCC‑TGTCCATTAAAGGGACCGTGGCATCCTGCCGGCCTGCGGCAAGCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCg                                                       <  1:53215/238‑1 (MQ=255)                                                                                                                                                                                         cccGGTACG‑TCGTTGGCTTTCCAGG‑GCGGGATCCAGCGT‑TGAGTCCTGTGACCCCG‑GGTGCATTAAGGGTAACGTGGCATCCTGCCGGCCTCCGGCATACGCCACCCCAGGTGCTGGCGGTAATGCGCGGATTCGCGACTGGCGATATTGCCCAATCGAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTGCACTTTCGt                                                      <  1:382798/236‑1 (MQ=255)                                                                                                                                                                                                  acggttgccgTTTCAGAGGCAGGAACCCGCGT‑TAATGCCTGATCCCCCT‑GTTCCATTAAAGGGAACGGGGCACCCTGCCAGCCTTCTGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCGCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCTCCAGGAAGGATTCTTCACTTTCGTCATACTGgcg                                            <  1:709286/231‑1 (MQ=255)                                                                                                                                                                                                    cGTTCGGGTTCCAGA‑GCCGGATCCAGCGT‑GGATGCCTGATCCCCCT‑GTTCCATTAAGGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCa                                          <  1:882410/238‑1 (MQ=255)                                                                                                                                                                                                      ttCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGTTGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGTCCTCCTGCCAGCCTTCGGCACTCGCCACACCAGGTGCTGCCGGTAATGCGCGGCTTCGCGTCTGCCGCTATTGCCCAAACAAGGCAGTCAACTGCG‑CCTGAGAGAGAATGACGGCCCACTGGGTGGCAACGGTCACCGGGAAGATttctttcattccttcttccgagcgcaac                                        >  1:1178945/1‑215 (MQ=255)                                                                                                                                                                                                        cgcgTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGTCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCACCTGCG‑ACTGAGAAAGAATAGCTGCCATCTTGGTGGCAAGCGTCACCAGGAAGGCTTCTTCCCTTTCGTCATACTGGCCCACctc                                      >  2:1129020/1‑238 (MQ=255)                                                                                                                                                                                                        cgcgTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGATTAGCTGCCATCTGGGTGGCAAGCGCCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAActct                                     >  2:43456/1‑239 (MQ=255)                                                                                                                                                                                                         ccgcttccaagggcGGACCCAGGGT‑TGAGGACTGATACACCT‑ATCCCGTGAAAGGTAACGGGGCATCCGGCCAGCCTGCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAActct                                     <  1:407001/226‑1 (MQ=255)                                                                                                                                                                                                           gTTCCAGA‑CCCGGCTCCGGCTT‑TGATGCTTGATACCCCT‑GTTCCATTAAAGGTAACGGGGCATCCGGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGt                                  <  1:436198/239‑1 (MQ=255)                                                                                                                                                                                                              ccAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTAACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTg                                >  2:882410/1‑238 (MQ=255)                                                                                                                                                                                                               cgga‑gCCGGCTCCCGGTT‑TGATGCCTGATACCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCt                              <  2:509362/237‑1 (MQ=255)                                                                                                                                                                                                                     ccGGATCCAGCGT‑TGATGCCTGATACCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTac                          <  2:455317/238‑1 (MQ=255)                                                                                                                                                                                                                      cGGATCCAGCGG‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTac                          <  2:5263/237‑1 (MQ=255)                                                                                                                                                                                                                          tCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGGATAACTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGACATACTGGCGCAACCCCCGTTGCGGTacaaacac                    >  1:646037/1‑233 (MQ=255)                                                                                                                                                                                                                          tCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTAccacca                     <  2:742631/238‑1 (MQ=255)                                                                                                                                                                                                                           ccAGCGT‑TGATTCCTGATACCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAg                    <  2:406947/238‑1 (MQ=255)                                                                                                                                                                                                                            cAGCGT‑TGATGCCTTATACGCCG‑GTTCCATTAAAGTTAACGTGGCATCCTGCCAGCCTTCGGAAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTa                  <  2:1308562/239‑1 (MQ=255)                                                                                                                                                                                                                              gCGT‑TGGTGCCTGATACCCCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTaca                <  2:435556/239‑1 (MQ=255)                                                                                                                                                                                                                               cGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACCCCAGGTGCTGACGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGGG‑ACTGTGAAGGAATAGCTGCCCTCTGGGTGGCACGCTGCACCAGGAAGGATTCTTCACTTTCGCCATCCTGGCGAAACCCCCGttccggagccacagtatca                >  1:587112/1‑221 (MQ=255)                                                                                                                                                                                                                               cGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTaccc               >  2:85319/1‑237 (MQ=255)                                                                                                                                                                                                                               cGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGAGTTTTCCCTTTCGTCATACTGGCGCCACTCCCGTTGCGGGACCCCCACAacac               >  1:998392/1‑239 (MQ=255)                                                                                                                                                                                                                                     aTGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATCTTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAAAAGCTGCCATCTGGGAGGCAGGCGTCACCAGGAAGGATTCTTCACTTTCGTCATAAAGGGGCAACTCTCGTTGCCGTACCACCAGTA‑ACCaaacc         >  1:1172694/1‑236 (MQ=255)                                                                                                                                                                                                                                     aTGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAAct                                       >  2:946905/1‑210 (MQ=255)                                                                                                                                                                                                                                     aTGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTCCTGTACCACCAGACCACCAAGc          >  2:546545/1‑239 (MQ=255)                                                                                                                                                                                                                                          tGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGGCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGACCCACCAGAAAACCaacaaacgg     >  1:348963/1‑232 (MQ=255)                                                                                                                                                                                                                                            aTACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGCAATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCCTCTGGGGGGCAGGCGCCACCAGGAAGGGTTCTTCCATTCCGACATCCGGGCGCACCCCCCGTCGCTGTCCCACCCGtcacccaacaaacatcg   >  1:1053885/1‑223 (MQ=255)                                                                                                                                                                                                                                              ccacCT‑GTTCCATTAAGGGTAACGTGGCCTCTGGCCAGGCTTCGGGAAGCGCCACAACAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCATACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCg  <  2:796196/237‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              ACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATT‑CGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAA‑GATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAG‑TCGTTCGCGTTCCAGA‑GCCGGATCCAGCGT‑TGATGCCTGATACACCT‑GTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCG‑ACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCG  >  NC_000913/2967609‑2968078

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