Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,037,225	A→C	5.7%	intergenic (+8/‑35)	ileT → / → alaT	tRNA‑Ile/tRNA‑Ala

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,037,225	0	A	C	5.7%	196.9 / 4.3	94	intergenic (+8/‑35)	ileT/alaT	tRNA‑Ile/tRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base A (39/49);  new base C (4/2);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAG  >  NC_000913/4036950‑4037499                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                   tGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCa                                                                                                                                                                                                                                                                               >  1:518312/1‑281 (MQ=32)   ggTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCaaa                                                                                                                                                                                                                                                                             >  1:5277/1‑281 (MQ=32)   ggTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCaaa                                                                                                                                                                                                                                                                             >  2:1061524/1‑281 (MQ=32)   ggTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCaaa                                                                                                                                                                                                                                                                             >  2:446978/1‑281 (MQ=32)    gTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTATGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACgg                                                                                                                                                                                                                                                                                 <  1:668844/276‑1 (MQ=17)      tGCAAACGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACCCAGGCCTACAAAATTTGCACGGCAAAttt                                                                                                                                                                                                                                                                          >  1:842784/1‑281 (MQ=17)      tGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTATGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTTCATTCAGGCTTCCCAAATTTCCACGGCa                                                                                                                                                                                                                                                                               >  2:668844/1‑276 (MQ=12)        cAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGCCTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACTAAATTCGCACGGCAAAATTGa                                                                                                                                                                                                                                                                        >  1:248629/1‑281 (MQ=17)        cAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCCCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGGACACCTGATAAGGGTGAGGTCGGTGGGTTAAGTGCGATCAGGCCTAACAAATTTGCACGGCAAATTTGa                                                                                                                                                                                                                                                                        >  2:172099/1‑281 (MQ=255)          aaaGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGaa                                                                                                                                                                                                                                                                       <  2:895930/280‑1 (MQ=32)          aaaGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGaa                                                                                                                                                                                                                                                                       <  1:516952/280‑1 (MQ=32)           aagaagTGGGTGGCTTAACCTTGCGGAGGGCGATTAACACTTTGTGATTCATGACTAGGGTGAAGTCGTAACAACTTAACCGTAGGAGAACATGCGGTTGGATCACCTCCTAACCTTAAAGAAGCGTTCTGTGAAGTGCTCACACAGATTGTCTGAGGAAAA‑TGAGCAGTAAAACCTCTACGGGCTTGTAGCTCAGGTGGTTAAAGCGCACCCCTGGTAAGGATGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAag                                                                                                                                                                                                                                                                      <  2:724648/280‑1 (MQ=11)           aagaagTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAag                                                                                                                                                                                                                                                                      <  1:446978/280‑1 (MQ=32)           aagaagGTGGTGGCTTAACCGTCGGGAGAGGGCTTACCTCTTGGTGATTCATTACTGGGGTTAAGTCGTAACGAGGTAACCGTAGGGGAACCTGCGGTTTGATCACCTCCTTACCTTAAAGAAGCGTTCTGTGAAGTGCTCACCCGGGTTGTCTGGTGAAAA‑TGAGCAGTAGAACCTCGACGGGCTTGTAGCTCAGGTGGTTAGTGCGCACCCCTGATAGGGGTGAGGTCAGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAaga                                                                                                                                                                                                                                                                     <  2:326483/281‑1 (MQ=11)                       gcttgtcccgcGGGTGGGCGCTTTCGCAATTGTGATTGATGTCTGAAGTGAAGTCGTAACAAGGTAACCGTTGGGGAACCTGCGGTTAGATCACCTTCTTACCTTAAAGAAGGATTCTTTGAAGTGCTGCCACAGATTTTCTTATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTATCTATGGTGGTTAGAGAGCACCCCTTATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATGTGCACGGCAAATTTGAAGAGGTTTTAACTAc                                                                                                                                                                                                                                                         <  2:915301/271‑1 (MQ=11)                       gCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTa                                                                                                                                                                                                                                                          <  2:1002826/280‑1 (MQ=32)                          tAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACa                                                                                                                                                                                                                                                        <  2:757014/279‑1 (MQ=32)                              cctctgggaggCGCTTAGCTCGTGTTGATTCCTGACTGGGGTGCAGTCCTAACAAGGTAACCGTAGGAGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAACTGCTCACACAGATTGTCTAATGAAAA‑TGAACAGTAAAACCCCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGAGAAGGGTGTGGTCGGTGGTTCAAGTCCACTCAAGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTAt                                                                                                                                                                                                                                                  <  2:1165074/272‑1 (MQ=11)                              cTTCGGGGGGGCGCTTACCACTTTGTGTTTCGTGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTa                                                                                                                                                                                                                                                   <  2:774564/280‑1 (MQ=17)                              cTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTa                                                                                                                                                                                                                                                   <  1:554187/280‑1 (MQ=32)                               ttgcggATGGCGCTTACCACTTGGTGATTCAGGACTGGGGTGAAGTCGTAACAAGGTAACCGTATGGTAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTa                                                                                                                                                                                                                                                   <  2:406659/275‑1 (MQ=12)                                tctggAGGGCGCTTACCACTTTGTGATTCATGACTGGGATGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACa                                                                                                                                                                                                                                                        <  1:573686/270‑1 (MQ=17)                                 tttgagagagtgacgaacaagttgcggtgttcagtgtttgaatattcagctttatccagatttttaaagagcaaaacttcgcagtgaacctttgcaggtacactctgaagtattttatgttgaatcactacagagatggtggagctatgcgggatcgaaccgcagacctcctgcgtgcaaagcaggcgctctcccagctgagctatagcccgataacatgtagttaaaacctcttCAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATgg                                                                                                                                                                                                                                                >  2:946593‑M1/236‑281 (MQ=255)                                 cGGGAGGGCGCTTACCACTTTGTGCTTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGCGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCTTTCTTTGATGTGCTCACACAGATTCTCTGATGAAAA‑TGAGTAGTAAAACCCCTACAGGCTTTTAGCTCAGGTGGTTAGAGCGCACTTCTGTTAAGGGTGAGGCCGGTGGTTCAAGTCCACTCCGGCCTCCTCAATTTGCACGGCAACTTTGCAGAGGCTTTAACTCAATGTCATg                                                                                                                                                                                                                                                 >  2:448630/1‑279 (MQ=11)                                 cGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATg                                                                                                                                                                                                                                                 >  1:407256/1‑279 (MQ=32)                                  gggAGGGCGCTTACCACTTTGTGATTTATGACTGGGGTGAAGTCGTAATAAGGTAACCGTAGGGGAACCTGTGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTCTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTATAGGCTTTTAGCTCAGGTGGTTAGAGTGCATCCCTGAGAAGGGTGAGGTAGGTGGTTCAAGTCTACTCAGGCCTACCAAATTTGCACGGCAAAGTTGAAGAGGCTGTAACTAAATGTTAtcg                                                                                                                                                                                                                                                >  2:573686/1‑277 (MQ=11)                                  gggAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATgggg                                                                                                                                                                                                                                              >  1:406659/1‑281 (MQ=32)                                           tcacgaacaactttcgttgttcagtgtttcaattttcagctttattcagatttttaaagagcaaaacttcgcagttatccttttcaggtacactctgaaggattgtttatttaagcactacagagatggtggagctatgcgggatcgaaccgcagacctcctgcgtgcaaagcaggcgctctcccagctgagctatagccccataacatgtagttaaaacctcttCAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCt                                                                                                                                                                                                                                      >  1:545544‑M1/226‑281 (MQ=255)                                                    ttGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCTGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGgag                                                                                                                                                                                                                            <  1:172099/281‑1 (MQ=17)                                                          attgttcagtgtttcaattttcagcttgatccagatttttaaagagcaaaacttcgtagtgaacctttgcaggtacactctgaagtattttttatttaaccactacagagatggtggagctatgcgggatcgaaccgcagacctcctgcgtgcaaagcaggcgctctcccagctgagctatagccccataacatgtagttaaaacctcttCAAATTTGCCGTGCTAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAgcgc                                                                                                                                                                                                                       >  1:866460‑M1/211‑281 (MQ=255)                                                                   ggggTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCa                                                                                                                                                                                                             <  2:5277/281‑1 (MQ=32)                                                                   ggggTGAAGTCGTAACAAGGTAACCGTAGAGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTTAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGATGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCa                                                                                                                                                                                                             <  1:482313/281‑1 (MQ=17)                                                                    gggTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCAc                                                                                                                                                                                                            >  1:147491/1‑281 (MQ=32)                                                                    gggTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTCCAAGTCCACTCAGGCCTCCCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCGTTGCAc                                                                                                                                                                                                            >  1:240253/1‑281 (MQ=17)                                                                       tGACGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCCCCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCCCCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCa                                                                                                                                                                                                             <  1:498672/277‑1 (MQ=17)                                                                         aagtcttgaacaGGTAACCGTAGGGGAACCTGCGGTTTGATCACCTCCTAACTTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCag                                                                                                                                                                                                        <  2:573391/269‑1 (MQ=12)                                                                           gTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCATGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATg                                                                                                                                                                                                                                                 >  1:381496/1‑237 (MQ=17)                                                                                    aGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGAGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGACTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTAAATGTTATGGGGCTATAGCGAACCTGGGAGACAGCCGGCTTTGCACGCAGGGGTTCTGCgg                                                                                                                                                                                             >  2:81197/1‑280 (MQ=11)                                                                                        aaCCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGa                                                                                                                                                                                        <  1:1061524/281‑1 (MQ=32)                                                                                             tAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATccc                                                                                                                                                                                    >  1:719225/1‑280 (MQ=32)                                                                                               gggtaaccgcgggTTGGGTAACCTCCTTACCTTAAAGAGTCGTTATTGTAAGTGCTCTCAAAGGTTGTCTGATGAAAA‑TGAGGAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAAAGCGCACGCGTAATGGGGGTGAGGTCGGTGGTTCGGGTGCACTCAGGCCTACCAAATTTGCACGACAAATTTTAAGATGTTTTAACTACATGTTATGGGGCTGTAGCTCAGCTGGGGGAGCGCCTGCTTTGCACACAGGAGGTCTGCGGTTCGATCCAGCa                                                                                                                                                                                 <  2:1072827/270‑1 (MQ=11)                                                                                               ggggAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGc                                                                                                                                                                                  >  1:990482/1‑280 (MQ=32)                                                                                               ggggAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCATACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACTCAGGAGGTCTGCGGTTCGATCCCGc                                                                                                                                                                                  >  1:102189/1‑280 (MQ=17)                                                                                                gggAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCa                                                                                                                                                                                 <  1:1159679/280‑1 (MQ=32)                                                                                                gagacccTGCGGTTGGATCACCTGCTTACCTTAAAGATGCGTTCTTGGAATTGCTGACACCGGTTGTTATATGAAAT‑TTAGCAGTAAAAACTCTACAGGGTTGTAGCTCAGGTGGTTAGGTCTTGCCCGTAGTAAGGGTGAGGTAAGTGGTTCAAGTCCACTCAGGCCTACCAAATGTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATTGCACGCAGGAGGTCTGAGGTTCGAACCCGCa                                                                                                                                                                                 <  2:248629/275‑1 (MQ=11)                                                                                                 ggAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATa                                                                                                                                                                               <  1:19005/281‑1 (MQ=32)                                                                                                   aagcTGCGGTTGCATCACCTCCCTACCTTAAAGAAGCGTTCTTTGAAGTGATCACTCAGATTGTCTGATGAAAT‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTATATCGCACCCCTGATAAGGGTGAGGTAGTTGGTTCAAGTCCACTCGGGCCTACCAAATTTTTACGGCAAATTTGCAGAGGTTTTAACTACATGTTATGGGGCTTTAGCTCAGCATGAAGAGCGCCTGCTTTTCACGCAGAAGCTCTTCGGTTTCTTCCCGCAtgc                                                                                                                                                                              >  2:1192745/4‑278 (MQ=11)                                                                                                   aaCCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCCAGTCCACTCAGGCCTACCAAATTTGCACGGCACATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGATCGCCTGCTTTGCACGCAGGCGGTCTTCGGTTCGTTCCGGctgag                                                                                                                                                                              >  2:444674/1‑276 (MQ=255)                                                                                                     ccTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATGGCt                                                                                                                                                                            >  2:199466/1‑280 (MQ=17)                                                                                                     ccTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCt                                                                                                                                                                            >  2:1018093/1‑280 (MQ=32)                                                                                                     ccTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCt                                                                                                                                                                            >  1:251511/1‑280 (MQ=32)                                                                                                      cTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTc                                                                                                                                                                           <  1:199466/280‑1 (MQ=32)                                                                                                        gctgTTGGGTCACCTCCTTACCATAAAGAAGCGTTTTTTGAAGTGCGCACACAGATGGTCTGATGAAAA‑TGAGCAGTAAAACCTCGACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGTTAAGGGTGGGGTCGGTGGTCCAAGTCCACTCAGGCCTACCAAATTTACACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTccac                                                                                                                                                                        <  2:240253/278‑1 (MQ=11)                                                                                                         tgtttGGATCACCTCCTTACCTTAAAGAAGCTTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTAAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGGCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTccac                                                                                                                                                                        <  1:739509/277‑1 (MQ=17)                                                                                                         gggTTGAGTCACCCCCGTACCCTAAAGAAGGTTTCTGTGAAGTGGTGACAAGGATTGTCTTATCAAAA‑TGGACAGTAAAGCATTTAAAGGTATGTAGCGCAGGTGTTTAGAGCGCCACCCTGATAAGGGTGGGGTCGGTGGTTCTAGTCCACTCAGGCCTACCAAAGTTGCACGGCAAATTTGGAGAGTTTTTAACTACATGTTATGGAGCTATAGCTCAGCTGGGAGGGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTccac                                                                                                                                                                        <  2:102189/279‑1 (MQ=11)                                                                                                         cggatagATCACCTCCGTACCTGAAAGAAGCGTTCTTTGAAGTGCTCACACTGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGAGCCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCGCTCAGGCCGACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGGGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTccac                                                                                                                                                                        <  1:432264/274‑1 (MQ=11)                                                                                                         cGGTTGTATCAGCACCTTACCTTAAAGAAGGCTTCTTTAAAGTGCTCACACAGTTTGTCTGATGAAAA‑TGAGCAGTAAAACCTCAACGAGCTTTTAGCCCAGGTAGTTAGAGCGCACCCCGGATCAGGGTGAGGTCGGGGGTTCAAGGCCACTCAGGCGTACCAAATTTGCACGTCAAATTTGAAGGGGTTTTAACGACATGTTATGGGGCTATATCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGGGGTCTGCGGTTCGATCCCGCATAGCTccac                                                                                                                                                                        <  2:1143738/280‑1 (MQ=11)                                                                                                         cGGTTGGATGACCTCCTTAACTTCAAGAAGCGTGATTTGAAGTGCTCAAACGGATTAGCTGATGTAAA‑TGGGCAGGAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCGTTATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTcgac                                                                                                                                                                        <  2:842784/280‑1 (MQ=11)                                                                                                         cGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTccacc                                                                                                                                                                       >  1:389831/1‑281 (MQ=32)                                                                                                         cGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTccac                                                                                                                                                                        <  2:389831/280‑1 (MQ=32)                                                                                                         cGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTccac                                                                                                                                                                        <  1:667516/280‑1 (MQ=32)                                                                                                         cGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTccac                                                                                                                                                                        <  1:81153/280‑1 (MQ=32)                                                                                                         cGGTTGGATCACCTCCTTACCTGAAAGAGGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAACAGTAAGACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTGTAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTccac                                                                                                                                                                        <  1:437283/280‑1 (MQ=12)                                                                                                         cGGTTGGATCACCTCCTTACCGTAAAGAAGCGGTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAAGTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTccac                                                                                                                                                                        <  1:241711/280‑1 (MQ=17)                                                                                                           gTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTCTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGCAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGCCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGTCTGCTCTGCACGCTGGAGGTATGCGGTTCGTTCTCGCATAGCCcccccaa                                                                                                                                                                     >  1:1143738/1‑280 (MQ=11)                                                                                                            ttGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCAt                                                                                                                                                                     <  2:407256/280‑1 (MQ=32)                                                                                                                aTCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCGACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGCGAGCGCCTGCTTTGCCGGAAGGAGGTCTGCGGTTCGATTCGGCCTAGACACCCCAtctcg                                                                                                                                                                 >  2:241711/1‑279 (MQ=11)                                                                                                                        cttaccttaAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCATCCATGATAAGGGTGAGGTCGGTGGTGCAAGTCCATTCAGGCCTACCAAATTTGCACGGCAAATATGAAGAGGTTTTAACTACATGTTATGTGGCTTTAGCACAGCTGGCAGAGCGCATGCTTTGCACGCCGGAGGTTTGGGGTTCGTGCCCGCATAGCTTCAACATTTCTGTAgggttt                                                                                                                                                        >  2:507017/1‑276 (MQ=11)                                                                                                                        cttaccttaAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCAtct                                                                                                                                                                   <  1:339317/270‑1 (MQ=32)                                                                                                                        cttaccttaAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGAtt                                                                                                                                                        >  2:805748/1‑281 (MQ=32)                                                                                                                          taccttaAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCa                                                                                                                                                                                                             <  1:1043121/226‑1 (MQ=32)                                                                                                                          taccttaAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCa                                                                                                                                                                                                             >  2:1043121/1‑226 (MQ=32)                                                                                                                          taccttaAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATtaa                                                                                                                                                      >  2:384895/1‑281 (MQ=32)                                                                                                                          taccttaAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCCCTGTAGTGATtaa                                                                                                                                                      >  1:718548/1‑281 (MQ=17)                                                                                                                          taccttaAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGCTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGTTATAGCTCAGCTGGGAGAGCGCCTGATTTGCACGCACGAGGTCTCGGGTTCGATTCCGCATAGCCCCACCACCTCTGTACTGAttct                                                                                                                                                      >  2:339317/1‑279 (MQ=11)                                                                                                                                 ggtacactctgaagtattttttatttaatcactacagagatggtggagctatgcgggatcgaaccgcagacctcctgcgtgcaaagcaggcgctctcccagctgagctatagccccataacatgtagttaaaacctcttCAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAgcgc                                                                                                                                                                                                                       <  2:326993‑M1/71‑1 (MQ=255)                                                                                                                                 ggtacactctgaagtattttttatttaatcactacagagatggtggagctatgcgggatcgaaccgcagacctcctgcgtgcaaagcaggcgctctcccagctgagctatagccccataacatgtagttaaaacctcttCAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAgcgc                                                                                                                                                                                                                       >  1:326993‑M1/140‑210 (MQ=255)                                                                                                                                     agCGTTCTTTGAAGTGCTCACTCAGATTGTCTGATGAAAATTGAG‑AGTAAAACCTCTGCAGGCTGGTAGCTCAGGTTGTTAGAGCGCACCCCTTATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGTCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCa                                                                                                                                                                                 <  2:558493/243‑1 (MQ=11)                                                                                                                                       cGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATa                                                                                                                                                                               >  1:558493/1‑243 (MQ=32)                                                                                                                                         ttttttGAATTGCTCACACAGATTGTTTGAGGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGGGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTc                                                                                                                                        <  2:1130742/277‑1 (MQ=9)                                                                                                                                           cTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTc                                                                                                                                        >  1:1088079/1‑278 (MQ=31)                                                                                                                                                                                                                                                               cACTCAGGCCTACCAAATTTGTCCAGACAATTTTAAGGGGTTTTAACTACATGTTGTGGGGCTATAGCTCAGCTGGGGAAGTGCCTGCTTTGCACACAGGGGTTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATTAAAAATCCTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAGAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAaca                   <  2:611396/280‑1 (MQ=11)                                                                                                                                                                                                                                                               cACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAaca                   <  1:1101101/280‑1 (MQ=11)                                                                                                                                                                                                                                                               cACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGAATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAaca                   <  1:1053994/280‑1 (MQ=11)                                                                                                                                                                                                                                                                aCTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAacac                  <  1:48186/280‑1 (MQ=11)                                                                                                                                                                                                                                                                 cTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGAAACACga                >  2:220161/1‑281 (MQ=14)                                                                                                                                                                                                                                                                        cTACCAAATTTGCACGGCAAATTTGGAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa         <  2:990482/281‑1 (MQ=17)                                                                                                                                                                                                                                                                               ATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaa    >  2:251421‑M1/1‑279 (MQ=255)                                                                                                                                                                                                                                                                                TTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCg                                                                                                                                                                                   >  2:1200895‑M1/1‑103 (MQ=255)                                                                                                                                                                                                                                                                                TTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCg                                                                                                                                                                                   <  1:1200895‑M1/103‑1 (MQ=255)                                                                                                                                                                                                                                                                                TTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAACTTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaag  <  1:546503‑M1/280‑1 (MQ=255)                                                                                                                                                                                                                                                                                TTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGTTGAATCgaaag  <  1:508186‑M1/280‑1 (MQ=255)                                                                                                                                                                                                                                                                                    CACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACtt                                                                                                                                         >  2:784845‑M1/1‑141 (MQ=255)                                                                                                                                                                                                                                                                                    CACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACtt                                                                                                                                         <  1:784845‑M1/141‑1 (MQ=255)                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                   TGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTCACACAGATTGTCTGATGAAAA‑TGAGCAGTAAAACCTCTACAGGCTTGTAGCTCAGGTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGTGGTTCAAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAG  >  NC_000913/4036950‑4037499

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