Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	9,668	C→T	G121G (GGC→GGT)	mog →	molybdopterin adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	9,668	0	C	T	97.5%	252.8 / NA	81	G121G (GGC→GGT)	mog	molybdopterin adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (41/38);  minor base G (0/2);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCG  >  NC_000913/9426‑9906                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               gAACTGGAAACC‑GCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCg                                                                                                                                                                                                                                        <  2:293509/250‑1 (MQ=255)     tGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCaa                                                                                                                                                                                                                                     >  2:339280/1‑250 (MQ=255)        aaaCCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAg                                                                                                                                                                                                                                 <  1:241964/251‑1 (MQ=255)          aCCCGATTAATCCCCGATGAGCAAGCGATCATCGAGCACAAGTTGTGTGAGCTAGTGGCTGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTACCGCCCGGTGCGACGCTGGCAGTAGCGGACAGCGAGATGCCGGGCTT‑TGGTGTACCGATGCGCCAGATAAACCTGCATTTGGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAgac                                                                                                                                                                                                                               >  2:249565/1‑249 (MQ=255)              gCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTg                                                                                                                                                                                                                           >  2:520273/1‑251 (MQ=255)               cTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTg                                                                                                                                                                                                                           >  2:10752/1‑250 (MQ=255)                 tAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGAt                                                                                                                                                                                                                         <  1:452844/250‑1 (MQ=255)                      cccGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGACATGAGTTGCCATCTGGTGCTCCCCACGGGCGGAACTGGACCGGCGCGCCGTGACGTAACGCCCGATGCGCCGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGCTAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTa                                                                                                                                                                                                                    >  1:592982/1‑250 (MQ=255)                        cGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACt                                                                                                                                                                                                                 <  1:140029/251‑1 (MQ=255)                        cGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACt                                                                                                                                                                                                                 >  2:192175/1‑251 (MQ=255)                           tGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTAc                                                                                                                                                                                                              >  1:100363/1‑251 (MQ=255)                                aGGCGATCATAGAGAAACCGTTGTGTGAGATAGTGGATGAAATCAGTTGCCATATGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGCCGTCACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCCGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGt                                                                                                                                                                                                         >  2:85740/1‑251 (MQ=255)                                 ggCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTc                                                                                                                                                                                                        >  2:39808/1‑251 (MQ=255)                                   cGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAg                                                                                                                                                                                                      <  2:267396/251‑1 (MQ=255)                                    gATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTACCTTACCCGGTCAGc                                                                                                                                                                                                     >  2:26291/1‑251 (MQ=255)                                     atcatcGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGcc                                                                                                                                                                                                    <  1:339280/251‑1 (MQ=255)                                           gAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTc                                                                                                                                                                                              >  1:257045/1‑251 (MQ=255)                                              cAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAt                                                                                                                                                                                           >  1:220056/1‑251 (MQ=255)                                                  cGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTaa                                                                                                                                                                                        >  1:419594/1‑250 (MQ=255)                                                      gtgtgAGCTAGTGGATGAAATGCGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAag                                                                                                                                                                                      >  2:257675/1‑248 (MQ=255)                                                          gagctcgaggATGAAATGAGTTGCCATCTGGCGCTCACCCCGGGCGGAACTGGACCGGAGCGACGTGACGTAACGCCAGATGCGACGATGGCAGTAGCGGACCGCGAGATGCCTGGCCT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGCTACTTTCGCGTCAGGTTGGTGTGCTTCGCAAACAGGCGCTGATCCTTAACTTCCCCGGTCAGCCGAAGTCTATTAAAGAGACGCt                                                                                                                                                                               >  2:419196/9‑251 (MQ=255)                                                               agtggATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGCCGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGa                                                                                                                                                                            >  1:233099/2‑249 (MQ=255)                                                                    aTGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCTCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGtgtg                                                                                                                                                                     <  2:386005/251‑1 (MQ=255)                                                                          tGAGTTGCCACCTGGAGCTCACCACGGGCGGAACTAGCCCGGCGCGTCGTCACGTAACGCCCACTGCGACGCTGGCAGGAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACGGCGCTCCTTTCGCGTCAGGTTGGTGTGATTCGCTAACAGGCGCTGATCCTTACCTTACCCAGACAGCCGAAGTCTATTAAAGAGGCGCTCGAAGGTGAGAAggcc                                                                                                                                                               >  2:62583/1‑249 (MQ=255)                                                                           gAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCGGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGAc                                                                                                                                                               <  2:481117/250‑1 (MQ=255)                                                                              ttGCCATCTGGTGCTCACCCCGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTg                                                                                                                                                           >  1:460749/1‑251 (MQ=255)                                                                                  cATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCggg                                                                                                                                                       >  1:539025/1‑251 (MQ=255)                                                                                   atcgggTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCGGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGt                                                                                                                                                      <  2:378128/247‑1 (MQ=255)                                                                                   aTCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCCGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGATGCGGGt                                                                                                                                                      >  2:354974/1‑251 (MQ=255)                                                                                       ggTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAAc                                                                                                                                                   <  1:192175/250‑1 (MQ=255)                                                                                       ggTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAAc                                                                                                                                                   <  1:520273/250‑1 (MQ=255)                                                                                       ggTGCTCACCACGGGCGGAACGGGCCCGGCGCGTCGGGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCGGGAAGGTGTGAAGGACGCTGCGGGTAAc                                                                                                                                                   <  2:312779/250‑1 (MQ=255)                                                                                        gTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGt                                                                                                                                                 >  1:41754/1‑251 (MQ=255)                                                                                            tCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTtgtg                                                                                                                                             >  1:351445/1‑251 (MQ=255)                                                                                             cacctgggtcGGAACTGGCCTGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGGAGCGGCTCGCGAGATGCCTGGCTG‑TGGTGATCAGATGCGCCAGATAAGCCTGCATTTTGGACCAACTGCGAGCCTGTCGCGTCATGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCATCCGAAGTCTACTAAATAGACGCTGGAAGGTGGGAAGGATGCTGCGGGTAACGTtgt                                                                                                                                              <  2:618164/240‑1 (MQ=255)                                                                                               ccacGGGCGGAACTGGCCCGACGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGt                                                                                                                                           >  1:334288/1‑250 (MQ=255)                                                                                                 ccGGGCGGAACTGGCCCGGCGCGTCGGGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAATCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCGTATCCTTAACTGACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTAGGGGTAACGTTGTGGTaca                                                                                                                                        <  2:56757/250‑1 (MQ=255)                                                                                                  cGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTacac                                                                                                                                       >  1:366398/1‑251 (MQ=255)                                                                                                             ggcgccgtgcgTCGTGAAGTAACGCCCGATGCGACGTTGGCAGTCGCGTACCGCGCTATGCCTTGCTA‑TTGTGAACAGATGCGCCAGATACTCCTGCATTTTGTACCAACTGCGATCCTTTGGCGTCAGGGGGGTGTGATTCGCAAAGAGGCGCTGATCCGTAATGTACCCGGTCAGCTTAAGTCTAGTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTCCACGGTAttttt                                                                                                                              <  2:556853/241‑1 (MQ=255)                                                                                                                  ggCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGcg                                                                                                                       >  1:507486/1‑251 (MQ=255)                                                                                                                      cgtcgtGACGTAACGCCCGATGCGAGGCTGGCAGTAGCTGACCGCGAGATGCTGGGCTTGTGG‑GAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGGGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTATAGAGACGCTGGAAGTTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGcgtacc                                                                                                                   <  2:579811/251‑1 (MQ=255)                                                                                                                       gtcgtGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGcgtaccg                                                                                                                  >  1:26212/1‑251 (MQ=255)                                                                                                                            gacttAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCATGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATGCGCAAACAGGCGCTGATCCTTAACGTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTGGCCAGCGTACCGTACt                                                                                                              <  2:142216/246‑1 (MQ=255)                                                                                                                              cGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGGACCGTACTGCa                                                                                                           <  1:147045/251‑1 (MQ=255)                                                                                                                                         aTGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACGGCATTCAGTTGCTg                                                                                                <  2:100363/251‑1 (MQ=255)                                                                                                                                             gACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGGCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAg                                                                                            <  1:85740/251‑1 (MQ=255)                                                                                                                                                     cAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAc                                                                                    <  1:257675/251‑1 (MQ=255)                                                                                                                                                         agcggacatctagaTGCCTGGCTT‑TGGTGAACAGGTGCGCCAGATAAGCCTGCATTTTGTCCCTACTGCGATCCTTGCGTGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCGTAACTTGCCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGGGGTACACGGTATTTTTGCCAGCGTACCGTACTGCAGTCAGTTGCTGGAAGGGCCATACGTTg                                                                                <  2:398220/240‑1 (MQ=255)                                                                                                                                                          gCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGa                                                                               <  1:563832/251‑1 (MQ=255)                                                                                                                                                           cGGACCGCGAGATGCCAGGCTT‑TGGTGAACAGATGCGCCAGATAACCCCGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGATGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGCATTTTTGCCAGCGTACCGTACTGCATTCCGTTGCTGGAAGGGCCTTACGTTGa                                                                               >  2:436489/1‑250 (MQ=255)                                                                                                                                                                      aTGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCAccaga                                                                   <  1:503500/251‑4 (MQ=255)                                                                                                                                                                       tGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCAccaga                                                                   >  2:320665/1‑247 (MQ=255)                                                                                                                                                                          cTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCCGAAGTg                                                               >  1:384890/1‑251 (MQ=255)                                                                                                                                                                            ggCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGt                                                             >  1:483911/1‑251 (MQ=255)                                                                                                                                                                             gCTT‑TGGTGAACCGATGCGCCCGATAAGCCTGCATTTTTTACCCACTGCTATCCTTTCTCGTCAGGTTGGTGTGATTCGCCAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTTAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCTTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGtt                                                            <  2:483297/251‑1 (MQ=255)                                                                                                                                                                               tt‑tGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGACGTGGTTgc                                                          >  2:103021/1‑251 (MQ=255)                                                                                                                                                                               tt‑tGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTgc                                                          >  1:455944/1‑251 (MQ=255)                                                                                                                                                                                         cAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAg                                                 <  1:103021/251‑1 (MQ=255)                                                                                                                                                                                         cAGATGCGCCAGATAAGCCTGAATGTTGGGCCAACTGCGGTCCTTTCGCGTCAGGTTGGTGGGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACAGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCGGCATTCAg                                                 <  2:170262/251‑1 (MQ=255)                                                                                                                                                                                          aGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGa                                                <  2:507486/251‑1 (MQ=255)                                                                                                                                                                                            aTGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGAcc                                              >  2:102524/1‑251 (MQ=255)                                                                                                                                                                                              gcgcCAGAGAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCg                                             <  2:41754/250‑1 (MQ=255)                                                                                                                                                                                                gcCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGcgtacc                                                                                                                   >  1:31435/1‑178 (MQ=255)                                                                                                                                                                                                gcCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGcgtacc                                                                                                                   <  2:31435/178‑1 (MQ=255)                                                                                                                                                                                                  cAGATAAGCCTGCAAGTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGGGTGAGTCGCAAACAGGCGCTGATTCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAATGTTGTGGTACACGGTAGTTTTGCCAGCGTACCGTACTGCATGCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGTTTGCTGCATTCAGACCGAagag                                        <  2:30689/251‑1 (MQ=255)                                                                                                                                                                                                     ataaGCCTGCATTTTGTTCCATCTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCTGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACTTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGc                                     <  1:420235/248‑1 (MQ=255)                                                                                                                                                                                                       aaGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCa                                    >  1:507196/2‑250 (MQ=255)                                                                                                                                                                                                          ccTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGAc                                <  2:1314/251‑1 (MQ=255)                                                                                                                                                                                                           cTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGAcg                               >  1:312972/1‑251 (MQ=255)                                                                                                                                                                                                             gCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTTCCCGG‑CAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGa                            <  1:348320/251‑1 (MQ=255)                                                                                                                                                                                                             gCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGAcg                               <  2:460749/249‑1 (MQ=255)                                                                                                                                                                                                                   tGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGtt                        <  1:353149/250‑1 (MQ=255)                                                                                                                                                                                                                    gTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGCTTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGAGGTTGCTGCATTCAGACCGCAGAGTGCAAGACGCGACGtt                        >  2:581030/1‑249 (MQ=255)                                                                                                                                                                                                                             gCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATaaa              <  1:26291/250‑1 (MQ=255)                                                                                                                                                                                                                              cGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTGACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATaaaaa            <  1:480779/251‑1 (MQ=255)                                                                                                                                                                                                                               gATCCTTTCGCGTCAGGTTTGTGTGATTCGCAAACAGGCGCTTATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCGGCGGGTAACGTTGTGGTACACGGTATTTTTGCCATCGTACCGTTCTTCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGTGTGCTGCAATCAGACCGAAGAGTGCAATACGCGACTGTAGCGAATaaaaaa           <  2:89175/251‑1 (MQ=255)                                                                                                                                                                                                                               gATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGCTCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTAAACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATaaaaa            >  2:623371/1‑250 (MQ=255)                                                                                                                                                                                                                                   cTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCAATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGACGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATcc       >  1:80811/1‑251 (MQ=255)                                                                                                                                                                                                                                    tttCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATccc      >  2:496398/1‑251 (MQ=255)                                                                                                                                                                                                                                      tCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATcccc     <  1:542528/250‑1 (MQ=255)                                                                                                                                                                                                                                         cgTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCAGGTCAGCCGAAGTCTATTAAAGAGACGATGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCg  >  1:359587/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               GAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCG  >  NC_000913/9426‑9906

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