Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,721,136	C→T	L38F (CTT→TTT)	ydhI →	DUF1656 domain‑containing protein YdhI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,721,136	0	C	T	91.2%	217.9 / 9.2	80	L38F (CTT→TTT)	ydhI	DUF1656 domain‑containing protein YdhI
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base T (32/41);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

GCTAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCG  >  NC_000913/1720894‑1721372                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             gCTAAAACAACTGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGt                                                                                                                                                                                                                                       >  2:421619/1‑250 (MQ=255)    aaaaCAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTGTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTacac                                                                                                                                                                                                                                   <  2:354669/251‑1 (MQ=255)    aaaaCAAATGTATTGCTGCATTTGGGTACCGTCAGACATATTTTTCAGAAATGGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCTGAAAGTAATCTACTTTGTGAAGTGTATG‑CTCAA‑GGCAACAGGATTGCCCTTACAAGATCTG‑TTGTTCGGTGCGTCCGTCTACTTTCCGCCGTGTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTacac                                                                                                                                                                                                                                   <  1:649554/251‑1 (MQ=255)    aaaaCAAATGTAGTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAAGTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGAGCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTac                                                                                                                                                                                                                                     <  1:80652/249‑1 (MQ=255)     aaaCAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACAcc                                                                                                                                                                                                                                  >  1:421968/1‑251 (MQ=255)       aCAAATGTATTGCTGCATTTGGTTACCGGCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTACGTCTACTTTACTCCGTTTTTCAAAGCATTCGCGTGTGGATTCGTCATCTGGTTTGTCGTACACCGc                                                                                                                                                                                                                                >  2:800173/1‑251 (MQ=255)                       aTTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTgg                                                                                                                                                                                                                 >  1:770253/1‑250 (MQ=255)                          ttgTTACCGTCAAACATATTTTTCAGAAATTGCGCGTAAATTGTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATGCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGTATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCTTCTTCTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATc                                                                                                                                                                                                              <  2:235935/248‑1 (MQ=255)                                ccGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATCCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTTGATTCGTCATCTGGTTTGTCGTACTCCGCCTGCTTCGTGGCTGGATCTACGCCg                                                                                                                                                                                                       <  1:143071/251‑1 (MQ=255)                                  gTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGt                                                                                                                                                                                                     <  1:448338/251‑1 (MQ=255)                                   tCAGACATATTTTTCAGAAATTGGGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCCTCTGGCTTGTCGTACACAGTCTGCTTCGTGGCTGGCTCTACGcagttg                                                                                                                                                                                                    >  2:318231/1‑246 (MQ=255)                                   tCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTg                                                                                                                                                                                                    >  2:48542/1‑251 (MQ=255)                                     agtatattttttCATAAAGTGCGCGTAAATTGTTCGCACTATAAGAATATTTAGTAATCTAACGCAATATCGTCGGTCGTAAAGGAATCTAATTTGTGAATTTTAGG‑CTCAT‑TGCAACATGATTGCCCTTACAAGATCTG‑GTTTTCGTGGCGTCCGTCTACTTTCCTCCTTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGCCGTACCCCGCCTGCTTCGTGGCTGGATCTACGCCGGAGAc                                                                                                                                                                                                  <  2:562284/244‑1 (MQ=255)                                     aGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGAc                                                                                                                                                                                                  <  2:770253/251‑1 (MQ=255)                                            tttttCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTCATCTAACGCAATATATTCGGTCGAAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAT‑TGCAAAAGGATTGCCCTTACCAGATATG‑GTGTTCGGTGCGTCCGTCTACTTTCTGCCGTTTGTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTCCACCCCCTGCTTCGTGGCGGGATCTACCCCGGGGCCATCTGGc                                                                                                                                                                                           >  2:447430/1‑251 (MQ=255)                                              tttCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCa                                                                                                                                                                                          <  1:288069/250‑1 (MQ=255)                                                tCAGAAAATGCGCGAAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGGCAGCTGGTGTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATc                                                                                                                                                                                        <  1:132120/250‑1 (MQ=255)                                                cctgcacatGCGCATAAATTTTTCGTAGTTAAAGAATATTTATTAATCTATCGCAATATAGTCGGTCGTAAAGGAATCGACTATGTGATGTTTATG‑CTCAA‑TTCAACAGGATTGCGGTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTCATTGCCTCCGTTTTTCAAAGCATTCGCGTTTGGATGCGTCATCTGGTTTGTCGTCCACCGCCTGCTTCGTGGGTGTATCTACGCCGTTGACATCTGGCATcc                                                                                                                                                                                       <  2:223282/243‑1 (MQ=255)                                                 cAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATccc                                                                                                                                                                                      <  1:753679/251‑1 (MQ=255)                                                   gAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACACGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGAGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATACCtt                                                                                                                                                                                    >  2:228056/1‑251 (MQ=255)                                                               aaaTTTGTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGATTCTGCTTTGTGAA‑TTTTTGTCTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTGCGGTGCGTCTGTCTCCTTACCACCGTGTTTCAAAGCAGCCGCGTTTGGATTCGCCTTCTGGTTTGTCGTACACCGCCGGCTGCGTGGCTGGATCTACGCCGTTTACATCTGGCATACCTTGGTAATGGAtt                                                                                                                                                                         <  2:413194/250‑1 (MQ=255)                                                                 aTTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGTTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTAt                                                                                                                                                                      <  2:880734/251‑1 (MQ=255)                                                                  tttttCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCCTCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGGTCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGGTTTATc                                                                                                                                                                     >  2:294765/1‑251 (MQ=255)                                                                   ttttCGCACTTAAAGAATATTTATTAATCTAACGCAATTTATGCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTGCGGTGCGTCCGTCTACTTTCCGCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATc                                                                                                                                                                     <  2:477226/250‑1 (MQ=255)                                                                        gCACTTAAAGAATTTTTATTAATCTAACGCAATATATGCGGGCGTAAAGGAATCTTCTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGAACTG‑GTGGTCGGTGCGTCCGTCTTCTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCGGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGtt                                                                                                                                                               <  2:596615/251‑1 (MQ=255)                                                                        gCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGAGCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGtt                                                                                                                                                               <  1:421619/251‑1 (MQ=255)                                                                          aCTTAAAGAATATTTATTAATCTAACGAAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGGTTTTG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTGGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTg                                                                                                                                                             <  2:283598/251‑1 (MQ=255)                                                                                aGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGAttt                                                                                                                                                       <  1:570944/251‑1 (MQ=255)                                                                                 gAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGCTTTATCGCTGTTTGCGATTTg                                                                                                                                                      >  2:904579/1‑251 (MQ=255)                                                                                           taATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCtt                                                                                                                                            >  1:170514/1‑251 (MQ=255)                                                                                            aaTCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCtt                                                                                                                                            <  1:385054/250‑1 (MQ=255)                                                                                            aaTCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTg                                                                                                                                           <  1:776814/251‑1 (MQ=255)                                                                                              tCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGct                                                                                                                                         >  2:953921/1‑251 (MQ=255)                                                                                                tAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGctct                                                                                                                                       >  2:291593/1‑251 (MQ=255)                                                                                                aaaCGCAAAAAATTCGGTCGGAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGAGCTG‑GTGTTCGGTGCGCCCTTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGTTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGctct                                                                                                                                       <  2:581512/250‑1 (MQ=255)                                                                                                 aaCGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTg                                                                                                                                      >  1:438506/1‑251 (MQ=255)                                                                                                 aaCGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGTCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGct                                                                                                                                         <  1:318231/248‑1 (MQ=255)                                                                                                 aaCGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGct                                                                                                                                         <  1:48542/248‑1 (MQ=255)                                                                                                     cAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCa                                                                                                                                   >  1:798417/1‑250 (MQ=255)                                                                                                     cAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCa                                                                                                                                   >  1:792987/1‑250 (MQ=255)                                                                                                      aaTATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGAAAt                                                                                                                                 >  1:750102/1‑251 (MQ=255)                                                                                                        tatatTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATAc                                                                                                                               >  1:330444/1‑251 (MQ=255)                                                                                                               ggTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGc                                                                                                                        <  1:276537/251‑1 (MQ=255)                                                                                                                     aaaGGAATCTTCTTTGTGAAGTTTATG‑CTCAA‑GGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCTTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGt                                                                                                                   <  1:720232/250‑1 (MQ=255)                                                                                                                           acttacgttGTGAAGTTTATG‑CTCAA‑TGCAAAAGGATTGCCCTTACAAGATTTG‑GTGTTCGGTGCGTCCGTCTACTTTCCCCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATg                                                                                                            <  2:860997/244‑1 (MQ=255)                                                                                                                               ccTTTGTGAAGTTTATG‑‑TCAAATTCAACAGGATTGCCCTGTCAATGT‑TGTGTGTTCGGTTCGTCCGGCTACCTTCCTCCGTTTTTCAAAGCATTCTCGTTTGGATTCGTCATCTGGTTTGTCGTTCCCCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCa                                                                                                         <  2:246172/249‑1 (MQ=255)                                                                                                                               aCTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGGAACTATGTCa                                                                                                         >  1:706678/1‑250 (MQ=255)                                                                                                                               aCTTGGTGAAGGTTATG‑CCCAA‑TGCAACAGGATTGCCCTTACAAGATCGT‑GTGTGTGGTGCGTCCGTCTAATTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTTCACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCa                                                                                                         <  2:593519/250‑1 (MQ=255)                                                                                                                                    gtgAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTaa                                                                                                    <  1:163043/250‑1 (MQ=255)                                                                                                                                     tgAAGCTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGCTTGCGTGGTAACTATGTCAATTaaa                                                                                                   >  1:172933/1‑250 (MQ=255)                                                                                                                                      gAAGTTTATT‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCTTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTTGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTaaaa                                                                                                  <  1:962246/250‑1 (MQ=255)                                                                                                                                       aaGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑TTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTTGCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCTATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACa                                                                                                <  2:103045/251‑1 (MQ=255)                                                                                                                                          tttATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAAtt                                                                                             >  1:66857/1‑251 (MQ=255)                                                                                                                                           ttATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAAtt                                                                                             >  1:926450/1‑250 (MQ=255)                                                                                                                                                         cAACATGATTGCCCTTACACAATCTG‑GTGGTCGGTGCGTCCGGCTACTTTCCTCCTTTTTTCAAAGCATTCGCTTTTGGATTCGTCATCTGGTTTTTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGTCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCGcaa                                                                                <  2:262516/251‑1 (MQ=255)                                                                                                                                                         cAACAGGATTGCCCTTACAAACTCTG‑GTGTTCGGTGCGTCCGACAACTTTCCTACGATTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCACTGTTTGAGCTTTGCGTTTGCCTTGCTATGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAAACATTAAGTCTTTATca                                                                                 >  2:947506/1‑250 (MQ=255)                                                                                                                                                           aCAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTcaac                                                                               >  1:465900/1‑250 (MQ=255)                                                                                                                                                                  tGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCAtt                                                                        <  1:429364/250‑1 (MQ=255)                                                                                                                                                                        tACAAGATCTG‑GTGTGCGGTGCGTCCGGCTTCTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTATTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCg                                                                  <  2:214664/250‑1 (MQ=255)                                                                                                                                                                         aCAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTAATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCggg                                                                >  1:609878/1‑250 (MQ=255)                                                                                                                                                                          cAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGt                                                                >  1:901321/1‑250 (MQ=255)                                                                                                                                                                                     tgtTCGGTGCGTCCGTCTACTTTCCTCCCTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGct                                                                                                                                         >  2:653884/1‑167 (MQ=255)                                                                                                                                                                                     tgtTCGGTGCGTCCGTCTACTTTCCTCCCTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGct                                                                                                                                         <  1:653884/167‑1 (MQ=255)                                                                                                                                                                                      gtttggggcgcccgtcTACTTTCCTCCGTGTTTCAAAGCATGCGCGTTTGGACTCGTCATCTGGTTTGTCGTACTCCGCCTTGTTCGTGGCTGGATCTTCGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCt                                                    <  2:637055/240‑1 (MQ=255)                                                                                                                                                                                      gtTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGCTTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCCTTGTAGCGGTAGTTGCGGtttt                                                    >  2:980101/1‑249 (MQ=255)                                                                                                                                                                                      gtTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATAATTGTAGAGGTAGTTGCGGTTCt                                                    >  1:121078/1‑251 (MQ=255)                                                                                                                                                                                      gtTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGGATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTc                                                     <  1:371289/250‑1 (MQ=255)                                                                                                                                                                                               gtccgtcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGAtg                                           >  1:614248/1‑251 (MQ=255)                                                                                                                                                                                                  cgtcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGAtggtg                                        >  1:213601/1‑251 (MQ=255)                                                                                                                                                                                                     cTACTTTCCTCCGTTTCTCAAAGCATTAGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTACAACAATTAAGTCTTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGCTAGTGGCt                                     >  2:592364/1‑251 (MQ=255)                                                                                                                                                                                                       aCTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCtgt                                   >  2:157977/1‑251 (MQ=255)                                                                                                                                                                                                              tccGTTTTTCAACGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTTTTGCCGGATGGTGGCTGTGGAattt                            >  2:61022/1‑250 (MQ=255)                                                                                                                                                                                                              tccGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTCGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAatt                             >  2:72009/1‑250 (MQ=255)                                                                                                                                                                                                               ccGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATTTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAattat                           <  1:904579/251‑1 (MQ=255)                                                                                                                                                                                                                 gTTTTTCAAAGCATTCGCGTTTGGATGCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTTGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAAttatta                         <  1:716299/251‑1 (MQ=255)                                                                                                                                                                                                                        aaaGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCa                   >  2:608364/1‑250 (MQ=255)                                                                                                                                                                                                                           gCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAAtct               <  2:943521/251‑1 (MQ=255)                                                                                                                                                                                                                             aTTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTcc             <  1:711212/251‑1 (MQ=255)                                                                                                                                                                                                                                      gggATTCGTCATCTGGTTTGTCGTTCACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGtcac     <  1:931420/249‑1 (MQ=255)                                                                                                                                                                                                                                        gATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAAGTATTACATGCAATCTCCGTGGACAcgcg  <  1:228056/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             GCTAAAACAAATGTATTGCTGCATTTGGTTACCGTCAGACATATTTTTCAGAAATTGCGCGTAAATTTTTCGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATG‑CTCAA‑TGCAACAGGATTGCCCTTACAAGATCTG‑GTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGCAATCTCCGTGGACACGCG  >  NC_000913/1720894‑1721372

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