Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,516,378	G→A	N100N (AAC→AAT)	fecA ←	ferric citrate outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,516,378	0	G	A	97.8%	267.1 / ‑4.2	93	N100N (AAC→AAT)	fecA	ferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (33/56);  total (35/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTG‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCA  >  NC_000913/4516138‑4516622                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         caatgtactCGCCAGGGCGGGGCTGCCGGTGCCGTTGGTTTCCGGGGGGCAGACGACAGGGAGGCGGTTAAGTACCTCA‑CGCATGGGGGTGGCGCCGGTTTTGGCGAAGTCCTCGCGGCGGAGCGCGTCACGCGCGCCAGCCTGTTCAAATACGGC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTGTTGGTGCGGGCGCGGGCTCCAGCGTGTAGCTA‑TTATTTCCCa                                                                                                                                                                                                                                               <  2:58532/243‑1 (MQ=255)    gttcatcccatgggcggggcgtccGGTGCCGTGGGTTCTCGGCTCGCTGACGCCAGTGTTGCGGTTTAGT‑CTCCAGCGCATGGTGTTTGCTCTGGGTTTGTCGAAACGCTCTCGGCGGATGTCGTCACGCGCGCCGGCATTTTCAAATACGTC‑GGTTTCACGCGCTTCACCCAGCCAGTCTCCGACCACGGTCAGGGCATTTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCTGCTA‑TTATTTCCCAGCg                                                                                                                                                                                                                                            <  2:377988/229‑1 (MQ=255)     ttCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGCTGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCgg                                                                                                                                                                                                                                           >  1:615151/1‑251 (MQ=255)           gCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCCCGCCCGCCCACAAGTTCAAATACGTCTGTTTCCA‑GCCCATCACCCAGCCAGGCGCCGAACCCGGTCCGGGCATCTTCTTTTGGGGCGGGGGCGCGGTCCACCGGCCCGCTA‑TTATTTTACAGCGGGTTTat                                                                                                                                                                                                                                     >  2:42390/1‑250 (MQ=255)            ccAGGTGGTGGCTGCCGGGGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTCAGTACCGCA‑CGCATGGTGGTTGTGCCGGTTTTGGCGAAATCCTTACGGCGGATCACGTCTCGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCATTCGCCGACCACGGGCAGGGCATCTTCTTTGGGTGCGGGCGCGGGGTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTGTTAcc                                                                                                                                                                                                                                    <  2:49924/251‑1 (MQ=255)             cAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCt                                                                                                                                                                                                                                   <  1:507699/251‑1 (MQ=255)                 gcggtgCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAg                                                                                                                                                                                                                               <  2:106147/246‑1 (MQ=255)                        tgtcGGTGCCTTTGTTGTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTGGCGCCGGTTTTGGCGAAAGCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCt                                                                                                                                                                                                                         <  2:165907/247‑1 (MQ=255)                           cGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGcc                                                                                                                                                                                                                      <  1:124985/250‑1 (MQ=255)                           cGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCCCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCCTCACCCCGCCAGTCCCCGGCCCCGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGCCCCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCttccc                                                                                                                                                                                                                     >  2:124985/1‑247 (MQ=255)                            ggTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGt                                                                                                                                                                                                                    <  1:168510/251‑1 (MQ=255)                               gCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCca                                                                                                                                                                                                                 <  1:120225/251‑1 (MQ=255)                                        ttCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAg‑ttgt                                                                                                                                                                                                        <  1:3601/250‑1 (MQ=255)                                         tCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAg‑ttgt‑t                                                                                                                                                                                                      <  1:520468/250‑1 (MQ=255)                                                  cTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCCCGCGCGCCAGCATGTTCAAATACGTC‑GTTTTTAAGCGCCTCACCCCGCCCGTCTCCGCACAAGGTCCAGGCATCTTTTTTTGGTGCGGGCGCGGGCGCCAGCGGCCAACTA‑TTATTTCCCAGGGGTTTTACCTTCAGTCCGCCGCCATCCAGCCG‑T‑GTGTGCAGGCCG‑cc                                                                                                                                                                                           >  2:525254/1‑250 (MQ=255)                                                      cGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGAGCACGTCACGGGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCCGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGGCCGCTGCCGTCCCCCAG‑TTGT‑TGCAGGCCG‑Ctctcc                                                                                                                                                                                       >  2:329412/1‑250 (MQ=255)                                                            ggATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTAACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGAcgtcg                                                                                                                                                                                 >  1:279266/1‑251 (MQ=255)                                                                 cGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCca                                                                                                                                                                                                                 >  1:292221/1‑217 (MQ=255)                                                                     tAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCCGCCCGTCGCCGACCACGGGCAGGGCATCCTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCCA‑TTATTTCCCAGCGGTTTTACCTGCCGTCCGCTGCGGTCCCGCAG‑TTTT‑TTCAGGCCG‑CTCTCGCCGTCGTAATcacct                                                                                                                                                                        >  2:392252/1‑247 (MQ=255)                                                                       aGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTg                                                                                                                                                                      <  1:433205/251‑1 (MQ=255)                                                                          aCCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAg                                                                                                                                                                   <  1:42390/251‑1 (MQ=255)                                                                            cTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCCGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCTGGc                                                                                                                                                                 >  1:476992/1‑251 (MQ=255)                                                                            cTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCATCACGGCGGATCACGCAACGCGCGCCAGCATGGTCCAATACGTC‑GTTTTCAAGCCCATCACCCAGCCCGAAGCCGACCACGGTCACGGCAACTTCTTTTGGGGCGGGCCCAGGCTCCAGCGTCCAGCCA‑TTATTTACCAGCGGGTTTACCATCAATCCGCCGCCGTCCAGCAA‑TTGT‑TTCCGGCCGGCT‑TCCACCTCGTAATCGCCGTGCatgc                                                                                                                                                                 >  2:656625/1‑248 (MQ=255)                                                                            cTCA‑CGAATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCCCCAGCATTTTCAAATACGTC‑GTTTTCACGCGCCTCACCCAGCCAGTCGACGACCACGGTAAGGGCATCCTCTTTTTGGGCGGGCGCGGGGTACCGCGGCCAGCTA‑TTATTTCCCCGCGGTTTTACCCGCAGTCCGCTGCCGTCCACGAG‑TTGT‑GGCAGGCCG‑ATCTCGACGGCGAAATCGCCgagccggg                                                                                                                                                                 >  2:627394/1‑244 (MQ=255)                                                                              cA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCg                                                                                                                                                               <  1:497128/251‑1 (MQ=255)                                                                                       gtggtTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATTTTCACATACGTC‑GGTTTCCCGCGCCTCCCCCAGCCCGTCACCGGCCCCGGTCCGGGCCACTTCTTTTGGTGCGGGCGCGGGCTCCAGACCCCCGATA‑TTATATCCCAGCGGTTTTACCTGCCATCCGCTGCCGTCCAGCCGCTTGT‑‑GCACGCCG‑CACTCGACGGCCGAATCGCCCAGCACGCCGGTGCTCTg                                                                                                                                                       >  2:215049/1‑251 (MQ=255)                                                                                             gcgcCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGCC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGcc                                                                                                                                                 <  1:67450/251‑1 (MQ=255)                                                                                                ccGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCGTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTGACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCcgcg                                                                                                                                              <  1:389457/251‑1 (MQ=255)                                                                                                       tGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCCGCAG‑TTGT‑TGCAGGCCG‑CTCTCGAc                                                                                                                                                                                     >  1:216917/1‑205 (MQ=255)                                                                                                       tGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCAGGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGATGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGAc                                                                                                                                                                                     <  2:216917/205‑1 (MQ=255)                                                                                                            aaaTCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCACGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGc                                                                                                                                  <  1:59072/251‑1 (MQ=255)                                                                                                                ccTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCCCCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTTCCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTTGCGTCa                                                                                                                              <  1:440487/251‑1 (MQ=255)                                                                                                                ccTCACGGCGGATCACGGCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCa                                                                                                                              <  1:493783/251‑1 (MQ=255)                                                                                                                     cggcggATCACGTCACGCGCGCCAGCATGTTCAAATACGGC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAAc                                                                                                                            <  1:166170/248‑1 (MQ=255)                                                                                                                       gcggATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTG‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGCCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGGCAACCgatg                                                                                                                       >  2:178311/1‑249 (MQ=255)                                                                                                                       gcggATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCgaga                                                                                                                       <  1:409903/251‑1 (MQ=255)                                                                                                                         ggATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAgg                                                                                                                     <  1:550427/251‑1 (MQ=255)                                                                                                                            tcacgtcacgCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCtt                                                                                                                                                    >  1:366027/1‑217 (MQ=255)                                                                                                                                  cacgcgcgccggcaggtgcaatgccgtc‑gtTTTCCCGCGCAGCACCCATCCAGTTGCCGACCACGGGCAGGGCATCTTCGTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCTTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCGCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCg                                                                                                            <  2:543408/227‑1 (MQ=255)                                                                                                                                  cacgCGCGCCCGCATGTTCAAATACGTC‑GTTTTCAAGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGACGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCg                                                                                                            >  1:513600/1‑251 (MQ=255)                                                                                                                                  cacgCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCGGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCATGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCg                                                                                                            <  1:32567/251‑1 (MQ=255)                                                                                                                                   acgCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGGGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCACCCGAGAGGGGAAATCCGc                                                                                                           >  2:183592/1‑251 (MQ=255)                                                                                                                                    cgcgcgCCATCATTGTCAAATACGTC‑GTTTTCACGCTCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCTTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGcg                                                                                                          <  2:268468/251‑2 (MQ=255)                                                                                                                                    cgcgcgCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGAcgt                                                                                                                                                                                   >  1:652460/1‑178 (MQ=255)                                                                                                                                    cgcgcgCCAGCATGTGCAAATACGTC‑GTTTTCACGCGCATCCCCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCTCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTGTACCTGCAGTCCGCTGCCGTCCAGCAG‑CTGT‑TGCAGGCCG‑CTCCCGAcgt                                                                                                                                                                                   <  2:652460/178‑1 (MQ=255)                                                                                                                                     gcgggcCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCCCCCAGCCAGTCGCCGACCACGGTCAGGGCACCTTCTTGTGGTGCGGTCGCGGGCTCCAGCGGCCAGCTA‑GTATTTCCCAGCGGTTTGACCTGCGGTCCGCTGCCGGCCAGCTG‑TTGT‑TGCAGGCCG‑GGCTCTACGTCGTAATCGCCGTGGAGGGCGTTGCTCTGCTTGCCGCGCGGCAGGCTGGCGTCAACCGAGAGGGTAAATCCGc                                                                                                           <  1:635725/245‑1 (MQ=255)                                                                                                                                       gcgcCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTCCCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgtg                                                                                                       <  1:88322/251‑1 (MQ=255)                                                                                                                                        cgcCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgtgt                                                                                                      <  1:348751/251‑1 (MQ=255)                                                                                                                                        cgcCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgtgt                                                                                                      <  1:535649/251‑1 (MQ=255)                                                                                                                                         gcCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgtgtg                                                                                                     <  1:451447/251‑1 (MQ=255)                                                                                                                                         gcCAGCATGTTCAAATACGCC‑GTTTTCCCGCGCATCACCCAGCCAGTCGCCGACCCCGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgtgtg                                                                                                     <  2:165392/251‑1 (MQ=255)                                                                                                                                           cAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCg                                                                                                   >  1:440677/1‑251 (MQ=255)                                                                                                                                           cAGCATGGTCAAATACGTC‑GTTTTCACGCGCATCACCCCGCCAGTCGCCGCCCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTGCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCTTCAACCGAGAGGGTAAATCCGCGGTGTGCg                                                                                                   <  2:440677/251‑1 (MQ=255)                                                                                                                                                tggtCAAATACGTC‑GTTTCCACGCGCATCACCCAGCCAGGCGCCGACCACGGTCAGGGCATCTTCGTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATGTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TTCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCa                                                                                                <  2:615151/246‑1 (MQ=255)                                                                                                                                                tGTTCAAATACGTC‑GTTTTCACGCGCATCCCCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGctct                                                                                                                                                        <  2:328299/193‑1 (MQ=255)                                                                                                                                                tGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGctct                                                                                                                                                        >  1:328299/1‑193 (MQ=255)                                                                                                                                                tGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCAt                                                                                               >  1:266109/1‑250 (MQ=255)                                                                                                                                                 gtgcAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCGCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATa                                                                                              <  1:347742/247‑1 (MQ=255)                                                                                                                                                    ataaTACGTC‑GTTTCAACGGTCATCACCCAGCCAGTCGCCGGCTCTGGTAAGGTCACTTTCTTTTTGTGCGGGCGCGGCCTCCAGCTTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTTCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTTTGCTTGCCGCGCGTCAGGCTGGCTTCAACCGAGAGGTTAAAGCCGCTGTTTGCGGCATACTg                                                                                           <  2:271235/248‑1 (MQ=255)                                                                                                                                                      aaTACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGAtt                                                                                        >  1:153474/1‑251 (MQ=255)                                                                                                                                                        aaCGTC‑GTTTTCCCGCGCATCACCCAGCCAGTCGCCGCCCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCAGACTGATTGa                                                                                      <  2:97235/250‑1 (MQ=255)                                                                                                                                                         ccgtc‑ggTTTCACGCGCACCACCCAGCCAGGCGCCGTCCACGGTCAGGGCATCTGCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCGGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAg                                                                                     <  2:546397/250‑1 (MQ=255)                                                                                                                                                           gtc‑gtTTTCACACGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTC‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TAGT‑TGCAGGCCG‑CTCTCTACGTCGTAATCGCAGTGCTGGCCGTTGCTCTGCTTGCCGCGCGCGTGGCTGGCGTCAACCGAGCGGGTAGTTCCGCTTTGTGCGGAATACtagttaagcg                                                                                   >  2:26030/1‑242 (MQ=255)                                                                                                                                                                       gcgcATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAAACGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGc                                                                        <  1:525254/251‑1 (MQ=255)                                                                                                                                                                       gcgcATCACCCAGCCAGTCGCCGACCACGGTCAGGGCACCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGc                                                                        <  2:573488/251‑1 (MQ=255)                                                                                                                                                                            tCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATcc                                                                   >  1:42684/1‑251 (MQ=255)                                                                                                                                                                             cccccAGCCGGTCGCCGACCACGGTCAGGGCATCTTCTGTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATcc                                                                   <  2:153474/248‑1 (MQ=255)                                                                                                                                                                                  agccagTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTACTCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCCACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGGTCCCggggc                                                             >  2:36220/1‑248 (MQ=255)                                                                                                                                                                                  agccagTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAACTCCGCTGTGTGCGGCATACTGTTTGATCGCTTTATCGGGCGATCCCGGTGc                                                             >  2:311642/1‑251 (MQ=255)                                                                                                                                                                                    ccagTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGCTTGAGCGCTTTATCGAGCGATCCCGGTGCgg                                                           >  1:147406/1‑250 (MQ=255)                                                                                                                                                                                         cGCCGACCACGGTCAGGGCATCTTCTTTTGGTGAGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGCCAACCGGGAGGGAAAATCCGCTGAGTGCGGCCAACTGATTGAGCGCTTTGTCGAGCGATCCCGGTGCGATAtt                                                       >  2:349064/1‑250 (MQ=255)                                                                                                                                                                                           ccGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGCCGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGCCAGGCTGGCGTCACCCGAGAGGGTATATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGTTATTAAc                                                    >  1:555782/1‑251 (MQ=255)                                                                                                                                                                                                cacgggtaGGGCATCTTCTTTTGGGGCGGGCGCGGGCTCTAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTCCATGCAGGCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTGGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAAGCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCtgtg                                               <  2:577507/244‑1 (MQ=255)                                                                                                                                                                                                cACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGa                                                                                      <  2:520567/212‑1 (MQ=255)                                                                                                                                                                                                cACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGa                                                                                      >  1:520567/1‑212 (MQ=255)                                                                                                                                                                                                   ggTCCGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTGCAGCTA‑TTATTTCCCAGCGGTTTTACCGGCAGTCCGCTGCCGGCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTTAGCGCTTTATCGAGCGAGCCCGGTGCGATATTAACCTGTgcag                                            <  1:513876/251‑1 (MQ=255)                                                                                                                                                                                                    gTCAGGACCTCTTCATTTGGTGCCGGCGCGGGCTCCAGCGTCCAGCTA‑TTCTTTCCCAGCGGTTTTAACTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCCGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTTCTCTGCTTACCGCGCGTCAGGCTGGCGTCAACAGAAAGGGTAAATCCGCTATTTGCGGCCTACTGATTGAGAGCTTTATCGAGCGATCCCGTTGCGATATTAAccagggc                                              >  2:591904/1‑243 (MQ=255)                                                                                                                                                                                                        gggccttttcttTGGGTGCGGGCGCGGGCTCCAGTGTCCAG‑TGATTATTTCCCAGCGTTTTTACCTGCAGTCCGCCGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCGGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCaaaa                                       <  2:147406/244‑1 (MQ=255)                                                                                                                                                                                                         ggCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAg                                      <  1:392252/251‑1 (MQ=255)                                                                                                                                                                                                         ggCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAg                                      <  1:107356/251‑1 (MQ=255)                                                                                                                                                                                                         ggCATCTTCTTTTGGGGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCATCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAg                                      <  2:202337/251‑1 (MQ=255)                                                                                                                                                                                                          gCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGc                                     >  1:679730/1‑251 (MQ=255)                                                                                                                                                                                                                 cgtgttgggcgggcgcgggcTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTGACCTGCAGGCCGCTGCCGGCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGCAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTTGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGaaaa                               <  2:460074/242‑1 (MQ=255)                                                                                                                                                                                                                           gggcgcgggcTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCgg                    >  1:658978/1‑251 (MQ=255)                                                                                                                                                                                                                              cgcgggcTCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCa                 <  1:17727/251‑1 (MQ=255)                                                                                                                                                                                                                                cgggcTCCCGCGGCCAGATT‑TTAGTTCCCAGCGGTTTTTCCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CCCTCGACGTCGGAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGTTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGc               <  2:289928/251‑1 (MQ=255)                                                                                                                                                                                                                                    cTCCAGCGACCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCCCTGTGTGCGGCATCCTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATCTTAACCTGTTCAGCAAAAGCGGACACCGAGAACCCGGCCacaggcc           >  2:572905/1‑245 (MQ=255)                                                                                                                                                                                                                                     tCCAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGgc            >  1:296369/1‑249 (MQ=255)                                                                                                                                                                                                                                      ccAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGgcagc         <  1:546163/251‑1 (MQ=255)                                                                                                                                                                                                                                      ccAGCGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGgca           <  1:329412/249‑1 (MQ=255)                                                                                                                                                                                                                                       caccGTCCAGCTA‑TTAGTTCCCGGCGGTTTTACCTGCAGGCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGAGCAGGCCGTTGCTCTGCTGGCCGCGCGTCAGGCGGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCTGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATAGTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGTTCAGCGgcagca        <  2:333707/248‑1 (MQ=255)                                                                                                                                                                                                                                          cGTCCAGCTA‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGc     <  1:666322/251‑1 (MQ=255)                                                                                                                                                                                                                                             ccAGCTA‑TTATTGCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCa  <  1:194643/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         AAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCTCA‑CGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGCGCCAGCATGTTCAAATACGTC‑GTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTG‑TTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAG‑TTGT‑TGCAGGCCG‑CTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCA  >  NC_000913/4516138‑4516622

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