Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,558,907	T→A	K207M (AAG→ATG)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,558,907	0	T	A	94.1%	279.6 / ‑5.4	85	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base A (48/32);  total (48/37)
Fisher's exact test for biased strand distribution p-value = 1.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

ACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCT  >  NC_000913/4558663‑4559152                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      aCCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTa                                                                                                                                                                                                                                                  <  1:237826/250‑1 (MQ=255)   cGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTg                                                                                                                                                                                                                                               <  2:198607/251‑1 (MQ=255)          gTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                                                                                                                                                                                                                         >  1:190399/1‑250 (MQ=255)           gcGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAt                                                                                                                                                                                                                                       <  2:336447/250‑1 (MQ=255)              gAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                                                                                                                                                                                                                                    >  1:85105/1‑251 (MQ=255)                    aGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAacac                                                                                                                                                                                                                              >  1:157199/1‑251 (MQ=255)                      ggTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATTCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACgg                                                                                                                                                                                                                            <  1:293490/251‑1 (MQ=255)                           cgcgTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACg                                                                                                                                                                                                                       >  1:219947/1‑251 (MQ=255)                           cgcgTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTACTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACg                                                                                                                                                                                                                       >  1:260761/1‑251 (MQ=255)                                 ccAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                                                                                                                                                                                                                 >  1:445175/1‑251 (MQ=255)                                   aGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                                                                                                                                                                                                               <  2:121046/251‑1 (MQ=255)                                        aaTACCCGCCTGCACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGCTTAccgc                                                                                                                                                                                                          >  1:97336/1‑251 (MQ=255)                                            cccGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAccgccg                                                                                                                                                                                                        >  2:398336/1‑249 (MQ=255)                                                ccTGAACGGCGCGGGCAATTCCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTGCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaaa                                                                                                                                                                                                  <  2:170041/251‑1 (MQ=255)                                                ccTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa                                                                                                                                                                                                   >  1:349531/1‑250 (MQ=255)                                                   gAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTACCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAccc                                                                                                                                                                                               >  2:274497/1‑250 (MQ=255)                                                     aCGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGGGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGATGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTTACGCCAGGTTTACCGCCGAGCAAACCGcc                                                                                                                                                                                             >  2:191944/1‑251 (MQ=255)                                                      cGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc                                                                                                                                                                                             >  2:78686/1‑250 (MQ=255)                                                         cgcgGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcg                                                                                                                                                                                         <  1:192702/251‑1 (MQ=255)                                                            gggCAATACCTTCGGCAGGGGCGACCGGTTCGTAGATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCGAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCGGAGCAAACCGCCAACGCgg                                                                                                                                                                                      <  2:214081/251‑1 (MQ=255)                                                               cAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGAt                                                                                                                                                                                   >  1:516347/1‑251 (MQ=255)                                                                                gCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAACCAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGCGCAAACCGCCACCGCGGGATTCTTCCGCCATATTGGc                                                                                                                                                                  >  2:444016/1‑251 (MQ=255)                                                                                      ggTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACACCGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTg                                                                                                                                                            <  1:409175/251‑1 (MQ=255)                                                                                          cGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGGGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATa                                                                                                                                                         <  1:314922/250‑1 (MQ=255)                                                                                           ggaaaTGCTGCTGTTGATATTGATGGTTGCGCCTTTGCGCGCGAACGCAAGCGCCGCCTCAAACAACGGTCTGTGGCTGTAAACGGGGGACGGCAGCCGCTTGTTGCTCGGCACCCTGCAGTTTTCCGGAGGGTAATAAATAGACTGTAACGCCATTTTACTGTCGCCCATGTGGCACACGGTTACGCCAGGTTTACCGCCGAGAAAACCGCAAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATaaa                                                                                                                                                       <  2:317254/248‑1 (MQ=255)                                                                                           gTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATaa                                                                                                                                                        >  1:231618/1‑250 (MQ=255)                                                                                            tCAATGCTGCTCGTGATATCGATGGTGCCGCCTTTGCACGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCCGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATCGGCCAGGTGATAAAc                                                                                                                                                      >  2:443631/1‑251 (MQ=255)                                                                                             cAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGATTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACg                                                                                                                                                     >  1:220512/1‑251 (MQ=255)                                                                                               aTGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTc                                                                                                                                                   >  2:426842/1‑251 (MQ=255)                                                                                                         tGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAAAACGGTGCGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGACATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCAAGGTTTACCGCCGCGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGACAGGGGATACACGTCCGGTGCGGct                                                                                                                                         >  2:19492/1‑250 (MQ=255)                                                                                                          gATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAg                                                                                                                                        <  1:402229/251‑1 (MQ=255)                                                                                                             atCGAGGGTGCCGCCTTTGCGCGCGAACGCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTGGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACGCGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAAc                                                                                                                                     <  2:137036/251‑1 (MQ=255)                                                                                                                 aTGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTg                                                                                                                                 >  1:416234/1‑251 (MQ=255)                                                                                                                  tGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTg                                                                                                                                 >  1:319321/1‑250 (MQ=255)                                                                                                                      gccgccTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCCGCAGCTTGCTGCTCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGGGATAAACGTCCGGTCAGGCAGAACGGGGATc                                                                                                                              >  2:92649/1‑249 (MQ=255)                                                                                                                      gccgccGTTGCGCGCGAACTCCAGCGCCGGCTCAAACAACGGTGCGTTGCGTTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCGGGTAATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCGGGTTTCCCGCCGAGCAACCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCag                                                                                                                            <  2:56557/251‑1 (MQ=255)                                                                                                                         gccTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGCTTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGAGATCaacgg                                                                                                                         >  2:346671/1‑247 (MQ=255)                                                                                                                                   gggAACACCAGCGCTTGCTCAAACAAGGGTGCGTTGCGGTTAACGTGTGGCGGCAGCAGCTTGCTGACCGGCACATGGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACTCCATTTTACTGTCGCCCATGTGGAGCACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCAttt                                                                                                               <  2:490442/249‑1 (MQ=255)                                                                                                                                   gcgAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCAttt                                                                                                               >  2:320022/1‑251 (MQ=255)                                                                                                                                   gcgAACCCCAGCGCCTGCTCAAACAACGGTACGTTGCGTTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCTCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCAttt                                                                                                               <  1:282102/251‑1 (MQ=255)                                                                                                                                     gAACTCCAGCGCCTGTTCAAACACCGGTAGGTTGCGTTTAACGGGGGTCGGCAGCAGCTGGCTGATCGGCACATCGCAGTTTTCCAGCAGTTCATAAATGGGCTGTACCGCCATTTTACTGTCGCCCATGGGGAACACGGTGGCGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCa                                                                                                             <  2:85105/251‑1 (MQ=255)                                                                                                                                      aaCTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCAc                                                                                                            >  1:431279/1‑251 (MQ=255)                                                                                                                                          ccAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGcc                                                                                                         >  1:152731/1‑250 (MQ=255)                                                                                                                                            aGCGCCTGCTCAAACAACGGTACATTGCGGTTACCGCGGGTAGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCAAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCACCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTACGGTGCGGCAGAACGGTGATCAGAGATGGCGAATTTCACGCCAAt                                                                                                      >  1:153491/1‑251 (MQ=255)                                                                                                                                               gccggCTCAAACAACGTTGCGTGGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGGGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATcac                                                                                                   <  2:264443/247‑1 (MQ=255)                                                                                                                                                ccTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCCGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGCGCTGGCGCATTTCACGCCAATcacc                                                                                                  >  2:277224/1‑250 (MQ=255)                                                                                                                                                    cTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgat                                                                                              >  1:228203/1‑251 (MQ=255)                                                                                                                                                      cAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTTATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGACCACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCGCGCCAATCACAcgatcg                                                                                            <  1:394576/251‑1 (MQ=255)                                                                                                                                                       aaaCAACGGTACGTTGTGGTTAAGGTGGGTGGGCGGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGAAGGTCATAAATAGGCGGGACCGCCATTTTACGGCCGCCCAGGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatcga                                                                                           <  2:250783/251‑1 (MQ=255)                                                                                                                                                       aaaCAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTgcc                                                                                                                                                                             >  1:57109/1‑169 (MQ=255)                                                                                                                                                       aaaCAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTgcc                                                                                                                                                                             <  2:57109/169‑1 (MQ=255)                                                                                                                                                       aaaCAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatcg                                                                                            >  1:483914/1‑250 (MQ=255)                                                                                                                                                         aaaacggtcggtgggggtaaaCGTGGGACGCAAGCACCTTGCTGACCGCCCCATCGCAGTTTTCCAGCAGGTCAAAAAGAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGACCACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAGCGGTGATCAGAGATGGCGCATTTCACGCCATTCACACGATCGATa                                                                                         <  2:152731/232‑1 (MQ=255)                                                                                                                                                            acGGTACGTGGCGGTTAACGGGCGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCCTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCCCGCCAATCACACGATCGATAATc                                                                                      >  2:462949/1‑251 (MQ=255)                                                                                                                                                             cGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCg                                                                                     <  1:10930/251‑1 (MQ=255)                                                                                                                                                                 aCGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCAc                                                                                 >  2:515350/1‑251 (MQ=255)                                                                                                                                                                             aCGTGGGTCGGCAGCAGCTTGCTGATCGGCCCATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCAc                                                                     <  1:402297/251‑1 (MQ=255)                                                                                                                                                                               gTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACgg                                                                   >  1:461797/1‑251 (MQ=255)                                                                                                                                                                                    gcGGCAGCGGCTTGCTGATCGGCACAGCGCAGTTTTCCAGAGGGTAAAAAATGGGCTGTAACCCCATTTTACTGTCGCCCATGTGGAACACGGGGACGCCGGGTTTCCCGCCGAGAAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAAccc                                                              <  2:335113/250‑1 (MQ=255)                                                                                                                                                                                      ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGt                                                            >  1:111999/1‑251 (MQ=255)                                                                                                                                                                                       gcagcaGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACGGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTa                                                           <  1:87925/251‑1 (MQ=255)                                                                                                                                                                                            agcTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGt                                                      <  1:121537/251‑1 (MQ=255)                                                                                                                                                                                                 gCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAACTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGACCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATACTTGACACGTCTTTTTCCACGGAACCCGTAATGGCGCggg                                                 >  2:276430/1‑251 (MQ=255)                                                                                                                                                                                                      acgccaccgcGCAGTTTTCCAGCAGGTCATAAAAAGGAGGTAACTCCTTTGTACCGCCGCCTATGTGGAACACGGGGAAGCAAGGTTTACGGCCGAGCAAACCGCAACCGCGGTATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGTGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAggg                                            <  1:325156/242‑1 (MQ=255)                                                                                                                                                                                                          cacaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACAt                                        >  2:281620/1‑251 (MQ=255)                                                                                                                                                                                                          accagcgtaGTTTCCCAGCGGGTCAAAAATGGGCTGAACCGCCTTTTTACTGTCCCCCAGGTGGACCACGGTGCGGCCGGGTTTCCCGCCGAGCAAACCGCCACCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAAGCGCCACGTCTTTTGCCACGGAACCCGTAATGGTGCGGGAAGGGACAt                                        <  2:98284/243‑1 (MQ=255)                                                                                                                                                                                                           acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATg                                       >  1:96856/1‑251 (MQ=255)                                                                                                                                                                                                             aTCGAAGTTTCCCAGCAGGTAAAAAATAGGCTCAACCTCCTTTTTACTGCGGCCCATGGGGAGCCCGGGGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGAt                                     <  2:431279/251‑1 (MQ=255)                                                                                                                                                                                                               cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATa                                    <  1:398336/250‑1 (MQ=255)                                                                                                                                                                                                               cGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTGACTGTCGCCCATGTGGAACACGGGGACGCCAGGTTTCCCGCCGAGCAAACCGCCACCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATa                                    <  2:331643/250‑1 (MQ=255)                                                                                                                                                                                                                   gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAgcgc                               <  1:373557/251‑1 (MQ=255)                                                                                                                                                                                                                     tttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGcc                              <  1:320022/250‑1 (MQ=255)                                                                                                                                                                                                                           gcagGTCATAAATAGGGTGTAACGCCATTTTACTTGCGCCCATGTGGAACACGGTGGCGCCAGGTTTTCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGAGAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCCCGCCAATCACACGAGCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGc                       <  1:443631/251‑1 (MQ=255)                                                                                                                                                                                                                             agGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCAt                     >  1:65744/1‑251 (MQ=255)                                                                                                                                                                                                                                 aaaaaatgggttGTAACGCCATTTTCCTGGCCCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAg                 <  2:157199/240‑1 (MQ=255)                                                                                                                                                                                                                                   tAAATCGGCTGTAACGCCATTTTACTTTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGGGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGGAATGGTGCGGGAAGGGCCATGATAAGCGCCGGTCAGCATCCAGgc               >  2:325036/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCCGAGATGGCGCATTTCACGCCAATCACACGACCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGACGGGACCTGATAAGCGCCGGTCAGCATcacggcc              >  2:196867/1‑245 (MQ=255)                                                                                                                                                                                                                                     aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGgcgc             >  2:37334/1‑251 (MQ=255)                                                                                                                                                                                                                                      aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCt            >  1:409557/1‑251 (MQ=255)                                                                                                                                                                                                                                      aTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCt            >  1:514134/1‑251 (MQ=255)                                                                                                                                                                                                                                             gTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGCCGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGc       <  1:78686/249‑1 (MQ=255)                                                                                                                                                                                                                                             gTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCct     >  1:269258/1‑251 (MQ=255)                                                                                                                                                                                                                                               aaCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCcttc   >  1:3387/1‑251 (MQ=255)                                                                                                                                                                                                                                               aaCGCCATTTTACTGTCGCCCATGTGGAAAACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGAGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCCCGTCTTTTTCCACGGAACCCGTAATGGAGCGGGACGGGACATGATAAGCGCCGGTCAGCATCCATGCGATGATGCcttc   >  2:304761/1‑251 (MQ=255)                                                                                                                                                                                                                                                aCGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGAAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATTGCCACGTCTTTTTCCACGGAACCCGTACTGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGCTGCcttct  >  2:224411/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      ACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCT  >  NC_000913/4558663‑4559152

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