BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,937,351	A→T	N20Y (AAT→TAT)	rbsK →	ribokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,937,351	0	A	T	99.5%	715.2 / NA	187	N20Y (AAT→TAT)	rbsK	ribokinase
Reads supporting (aligned to +/- strand): ref base A (0/1); major base T (94/91); minor base G (0/1); total (94/93)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

TGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCC  >  NC_000913/3937072‑3937635

tggagagAAAGGCGGCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATc                                                                                                                                                                                                                                                                                    <  1:343422/283‑1 (MQ=255)

   cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAACCGCAGGCAGCCTCGTTGTTCTTGGAAGCATTCATGCTGACCACCTTCTTTATCTTCACTCttt                                                                                                                                                                                                                                                                                 >  1:251050/1‑289 (MQ=255)

   cgcgAAAGGCGTCGAAACCGCAGATAAAGTGATGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCt                                                                                                                                                                                                                                                                                   <  1:84927/288‑1 (MQ=255)

       aaaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTcc                                                                                                                                                                                                                                                                              >  1:127945/1‑288 (MQ=255)

        aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCt                                                                                                                                                                                                                                                                             >  1:434860/1‑289 (MQ=255)

                  aaCCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCg                                                                                                                                                                                                                                                                   <  1:164483/288‑1 (MQ=255)

                       cAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGCTGTTCTTGGCAGCCTTAATGCTGACCACCTTCTTTATCTCCACTATTTTCCCACCCCAgcctcacccg                                                                                                                                                                                                                                                             >  2:164483/1‑280 (MQ=255)

                       cAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTT‑TGGACATCCCGAATATGCAAAACGCACGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTCACCACCTTCTTTATCTCCAATCTTTTCCTACTCCAGGCGAAccca                                                                                                                                                                                                                                                             >  2:264902/1‑288 (MQ=255)

                              aGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCgg                                                                                                                                                                                                                                                      <  1:345026/289‑1 (MQ=255)

                              aGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCgg                                                                                                                                                                                                                                                      <  1:836634/289‑1 (MQ=255)

                                agCTTAAAGGCGAGAAAGTTCAAGGTAAATGTCCGGTTTATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATc                                                                                                                                                                                                                                                                                    <  1:406059/257‑1 (MQ=255)

                                 gCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTGAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGt                                                                                                                                                                                                                                                     <  1:245351/287‑1 (MQ=255)

                                  cTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACCTTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTaa                                                                                                                                                                                                                                                   >  1:602744/1‑288 (MQ=255)

                                   tgtgaGGCGGGAAAGTTCAGGCGAAGTAGCGGGTTGGTCTTAAACAGGTTGTTAAGAAGTAGGTTGAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATAGAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACAGCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAAc                                                                                                                                                                                                                                                  <  2:743633/285‑1 (MQ=255)

                                        ggggagaAAGGTCAGGCGAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTa                                                                                                                                                                                                                                             <  2:434860/286‑1 (MQ=255)

                                          ggagaAAGTTAAGGCGAAGGATCCGGTGGATCTGAAACTGGTTGTTAAGCAGTAGTTTGAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTa                                                                                                                                                                                                                                             <  2:800691/286‑1 (MQ=255)

                                          cGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GAACATCCGGAATATGCACAACGCAGGCAGACTCGTTGTTCTTGGCAGCATTAATGCTGACCACACTCTTTATCTTCAATCTTTTCCTACTCCAGGCGTCACTCTCACCTGTACACACTATCa                                                                                                                                                                                                                                          >  2:231465/1‑289 (MQ=255)

                                                   tCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCAt                                                                                                                                                                                                                                  <  1:508133/289‑1 (MQ=255)

                                                                ccGTTTTATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGgc                                                                                                                                                                                                                    <  1:21477/289‑1 (MQ=255)

                                                                     tGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGaa                                                                                                                                                                                                               >  1:358307/1‑289 (MQ=255)

                                                                       aTCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATc                                                                                                                                                                                                                                                                                    >  2:344635/1‑218 (MQ=255)

                                                                          tGAAAAGGGTTGTTAAGCATGAGTTTTAATCAGGTTTTATGCACTGATGGTGACATAAAACCGGCATCGAAAGAGGAACGTGGAAGATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTATCCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAgg                                                                                                                                                                                                          <  2:787308/289‑1 (MQ=255)

                                                                              aatgtgtggtAAGGAGTAGTTTTAAGCAAGGTGTATGACCAGAAGGTGACAAAAATACGGCAACGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCtgt                                                                                                                                                                                                      <  2:251050/280‑1 (MQ=255)

                                                                                  gttgttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGCAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATAAGGCTGTGgcg                                                                                                                                                                                                  >  1:141652/1‑288 (MQ=255)

                                                                                  gttgttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTCTATCTTCAATCTTTTCCTATTCCAGGCGACACCGTAACCGGTAACCACTATCAGGTTGCTTTTTGCGGCAACGGCTCACATCAGGCTGTgggt                                                                                                                                                                                                  >  2:720454/1‑287 (MQ=255)

                                                                                   ttgttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAACGCAGCTAATCAGGCTGTGtct                                                                                                                                                                                                  >  1:784357/1‑288 (MQ=255)

                                                                                    tgttAATCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGgctg                                                                                                                                                                                                 <  1:421517/289‑1 (MQ=255)

                                                                                    tgttAAGCAGTAGTTTTAAGCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACTTCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGgctg                                                                                                                                                                                                 <  1:110562/288‑1 (MQ=255)

                                                                                        aaGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTgg                                                                                                                                                                                             <  1:561906/289‑1 (MQ=255)

                                                                                              tagtTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTa                                                                                                                                                                                        >  1:453738/1‑288 (MQ=255)

                                                                                               agtTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGc                                                                                                                                                                                      <  1:158758/288‑1 (MQ=255)

                                                                                               agtGTTAATCAGGTTTTATGAACTGATGGTGACATAAATACGTCATCGACAGAGGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGc                                                                                                                                                                                      <  2:54364/289‑1 (MQ=255)

                                                                                                           gTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCCCGCGCCACCTTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCCGACCACATTCTTTATCTCCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATc                                                                                                                                                                          >  1:636142/1‑288 (MQ=255)

                                                                                                           gTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGTTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCCGGGCGTAGCGGTGCGaagttc                                                                                                                                                                          >  1:399422/1‑284 (MQ=255)

                                                                                                           gTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCCAAACCGTAACCGGTAACCACTACCAGCTTGCCTTTGGCGGCAAAGGCGCGAATCAGGCTGTGTCTGCTGGGCGTTGCGGTTGGAATATc                                                                                                                                                                          >  2:222017/1‑288 (MQ=255)

                                                                                                            ttGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAatat                                                                                                                                                                           <  1:235467/287‑1 (MQ=255)

                                                                                                                  gACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTa                                                                                                                                                                   <  1:473618/289‑1 (MQ=255)

                                                                                                                     cTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGgctg                                                                                                                                                                                                 >  2:150609/1‑255 (MQ=255)

                                                                                                                     cTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGgctg                                                                                                                                                                                                 <  1:150609/255‑1 (MQ=255)

                                                                                                                     cTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCAAGTCGAAACCGTAACCGGTAACCACTATCAGGTTGCATGTGGCGGCCAAGGCTCGAACTAGGCTGTTGCTGCTGGGCGTAGCGGTGCGAGTATCGCGTTTATTg                                                                                                                                                                >  1:190710/1‑288 (MQ=255)

                                                                                                                      tGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAAACACTATCAGGTTGCCTTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCGGGGCGTAGCGGTGCGCATATCGCGTttgttg                                                                                                                                                                >  1:653759/1‑288 (MQ=255)

                                                                                                                        aTGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGGGGTGCGAATATCGCGTTTATTg                                                                                                                                                                >  1:79616/1‑286 (MQ=255)

                                                                                                                         tGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGt                                                                                                                                                           >  1:646473/1‑289 (MQ=255)

                                                                                                                            tGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACTCGCCACCCTAACACGGTGGCGCATGTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAAGCTTTTCCTACGCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTAc                                                                                                                                                         <  1:815325/288‑1 (MQ=255)

                                                                                                                              aCATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTTCTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACggg                                                                                                                                                      >  2:693796/1‑289 (MQ=255)

                                                                                                                                    aTACGTCATCGACCAATGAACGTGTAACATAAAGAAAACCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGtgatga                                                                                                                                                <  2:405661/289‑1 (MQ=255)

                                                                                                                                         tcagataCAGATGAACGGAGTATATCAAGCAAAGCAGGGGTCGCGCCACCCGAACACTGGGGCGCAGTTTA‑CGACAACACGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCATCAGTAATGCTGACCACATTCTTTATCGTCAATATTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCCCTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCa                                                                                                                                           <  2:61019/283‑1 (MQ=255)

                                                                                                                                          cATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCATGTTGCCTTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGtgatgtcggctt                                                                                                                                          >  2:705671/1‑282 (MQ=255)

                                                                                                                                                cAGATGAAATGGTAATATAAAGAAAAGCAGGGCACGGTCCACCCTAACCCGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCTTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTga                                                                                                                                    <  2:602744/289‑1 (MQ=255)

                                                                                                                                                          tgtAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCACTCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCCGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGTTTATTGCCGCATTGGTGAGAGCGTTCg                                                                                                                           >  1:675132/1‑288 (MQ=255)

                                                                                                                                                            tAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCAAATCAGGCTGTTGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGGGCGTTCGCca                                                                                                                        >  2:361306/1‑289 (MQ=255)

                                                                                                                                                             aagataAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCcag                                                                                                                       <  1:6706/286‑1 (MQ=255)

                                                                                                                                                             aacagaaaGAAGAGCACGGCACGCGCCACCCTATCACGGTGGCGCATTTTA‑GGACATCCCGAATATACAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACCCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCcag                                                                                                                       <  1:323747/284‑1 (MQ=255)

                                                                                                                                                                    aGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAACGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGTTGACGTCATTGGTGAGAGCTTTCGCAcgcagccc                                                                                                                 >  2:18454/1‑286 (MQ=255)

                                                                                                                                                                       aaaGCAGGGCACGCGCCACCCTAACACGGTGGCGCAGTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCg                                                                                                                <  1:303904/287‑1 (MQ=255)

                                                                                                                                                                         aGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATCTGGCGGCAAAGGCGCGAATCAGGTTGTGGCTGCTGGGCGGAGCGGTGCGAATAGCGAGTTTCTTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCTCCAGCTGCTTGCCat                                                                                                            >  1:873216/1‑288 (MQ=255)

                                                                                                                                                                             cggCACTCGCCACACTATCACGGTGGCGCAGTTTATGGACATCCCGAATATGCACAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGAt                                                                                                        <  1:537862/288‑1 (MQ=255)

                                                                                                                                                                                cACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGAGAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCCCTGCTAAc                                                                                                     >  2:302904/1‑289 (MQ=255)

                                                                                                                                                                                 aCGCGCCACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACAt                                                                                                   <  2:358307/289‑1 (MQ=255)

                                                                                                                                                                                 aCGCGCCACCCGAACACTGTGGAGCATTTTAT‑GACACCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCCGCATTAATGCTGACCACATTCTTTATCTTCAACCTTTTCCTACGCCAGGCGAAACCGTAACCGGTAACCACTCTCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACAt                                                                                                   <  2:487163/289‑1 (MQ=255)

                                                                                                                                                                                 aCGCGCAAACCTAACTCGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCGTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACAt                                                                                                   <  1:7524/289‑1 (MQ=255)

                                                                                                                                                                                      ccACCCTAACACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTTCTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGata                                                                                              <  1:693796/289‑1 (MQ=255)

                                                                                                                                                                                        aCCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACAt                                                                                                   >  1:110202/1‑283 (MQ=255)

                                                                                                                                                                                             aaCACGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCCGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGCCAGCATTGGTGAGGGCGTTCCCACGCAGCTCGCCACTGATAACATTGATTTTACTc                                                                                        >  2:817524/1‑288 (MQ=255)

                                                                                                                                                                                               cacGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACCGCATTGGTGAGAGCGTTCGCCAGAAGCTAGCCACGGGTAAGATTGATAGTACTCCg                                                                                      >  2:276583/1‑289 (MQ=255)

                                                                                                                                                                                               cacGGTGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGtgatga                                                                                                                                                >  2:21013/1‑230 (MQ=255)

                                                                                                                                                                                                   gTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGgctgct                                                                                                                                                                                               >  1:198999/1‑180 (MQ=255)

                                                                                                                                                                                                    tGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTa                                                                                                                                                          <  2:31818/216‑1 (MQ=255)

                                                                                                                                                                                                    tGGCGCATTTTA‑GGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCCGCGGCTCGCCCCGTATAACTTTGAATGTACTCGGGTCAGc                                                                                >  2:264119/1‑289 (MQ=255)

                                                                                                                                                                                                                tGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTa                                                                                                                                                                   >  1:172641/1‑195 (MQ=255)

                                                                                                                                                                                                                 ggACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCCCGCCACTGATAACATTGATATTACTTCGGTCAGCGTGATCAAAGGc                                                                    >  1:399922/1‑289 (MQ=255)

                                                                                                                                                                                                                  gcaaTCCGGAACATGCAAACCGAAGGCGGCCTTGTTGTTCTGGGCACAATTAAGGCTGACAACATGCTTTATCTTTAATCTTTCCCTACGCCAGGCGAAACAGTAACCGGTAACCACTAGCAGGTGGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGc                                                                    <  2:789852/285‑1 (MQ=255)

                                                                                                                                                                                                                   aCATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGc                                                                    <  1:705671/287‑1 (MQ=255)

                                                                                                                                                                                                                    cATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGATTATCGCGCTTATTGCCTGTACGGGCGATGACGGCTTTGGGGGAAGCGTTCGCCGGCAGCCAGCCACTGATTACATTGTTATTCATCCGGTCCGCGGGACTAAAGGCgcc                                                                 >  2:168027/1‑287 (MQ=255)

                                                                                                                                                                                                                       cccGAATATGCATAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTTTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTGATTGCCTGTACGGGTGATGACGGATTTGGTGATAGCGTTCGCCAGCCGCTCCCCACTGATAACAGTGATATTACTCCGGTTAGCGTGATCAACGGCCAATCa                                                              >  2:715529/1‑289 (MQ=255)

                                                                                                                                                                                                                       cccGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCa                                                              >  1:389079/1‑289 (MQ=255)

                                                                                                                                                                                                                           aaTATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAg                                                          >  1:556742/1‑289 (MQ=255)

                                                                                                                                                                                                                                 cAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTGTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAAAATTGATATTACTACGGTCAGCGTTGTAAAAGGCGAAACAACAgggaggg                                                    >  2:771561/1‑284 (MQ=255)

                                                                                                                                                                                                                                           gcagCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGtgat                                                                                                                                                  >  1:264129/1‑185 (MQ=255)

                                                                                                                                                                                                                                                cTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGCCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTAGGGGTGATGCCAGCATTGGTGAGCGCGTTCGCCAGCAGAGCGCGACTGAGACCATTGAGATTACTCGGGTCAGCGTGATCAAAGGCACAACAACAGGTGTGTAGTTGATTTGTGttc                                     >  2:64815/1‑288 (MQ=255)

                                                                                                                                                                                                                                                       ttCTTGGCAGCATTAATGCCAACCACATTCGTTATCTTCAATCTTTTCCTACGCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGa                              <  1:720454/289‑1 (MQ=255)

                                                                                                                                                                                                                                                       ttCTTGGCAGCATTAAGGCTGACCACCTTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGa                              <  1:331677/289‑1 (MQ=255)

                                                                                                                                                                                                                                                       ttCTTGGCAGCATTAAATCTGACCACAGTCTTTATCTTCAATCTTTTCCTACGCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGa                              <  1:771561/289‑1 (MQ=255)

                                                                                                                                                                                                                                                         cTTGGCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAg                            >  1:688700/1‑289 (MQ=255)

                                                                                                                                                                                                                                                            ggCAGCATTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTAGCCACTGGTAACATGGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACATGTGAGGCGCTGCTTTTTGGTAATGGCTAAGGTg                         >  1:525333/1‑289 (MQ=255)

                                                                                                                                                                                                                                                               agcaTTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTgag                       <  1:550069/288‑1 (MQ=255)

                                                                                                                                                                                                                                                                 caTTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCATCTCGCCACTGCTAACATTGATATGACTCCGGTCACCGTCATAAAAGGCGAATCAACAGGGGTGGAGCTGATTTTGTGTACTGTCGAAGGTgaggag                    >  2:639104/1‑286 (MQ=255)

                                                                                                                                                                                                                                                                 caTTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATCTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTCTACGGGTGATGACAGCATTGTTGAGAGCGTCCGCCCGCAGCTCGCCACTGATCACATTGCTATTAACACGGTCAGCATGATCCAAGCCGATCCATCAGATTTCGCGCTTCTTTCTGTTCATGGCGatagttaacat                    >  2:554496/1‑279 (MQ=255)

                                                                                                                                                                                                                                                                  aTTAATGCTGACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGaa                     <  1:89805/287‑1 (MQ=255)

                                                                                                                                                                                                                                                                   tacaTGCTGACCACATTACTGCTCTTTACTGGTTTCCTCGTACATGCGAACCCGTAACCGATAACCACTGTCAGGTCGCATTTGGCGGCAAAGGGGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTGTTGCCGGTACGGGTGAAGACAGCATTGGTGAGAGCGTTCGCTAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGt                  <  2:636142/286‑1 (MQ=255)

                                                                                                                                                                                                                                                                           gACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGAt                                                                                                        <  1:283485/195‑1 (MQ=255)

                                                                                                                                                                                                                                                                           gACCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGGTTTTTGTTAATGGCGAAGGTGCGAATGTCATCGGTa          >  1:857613/1‑289 (MQ=255)

                                                                                                                                                                                                                                                                            aCCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGc                                                                                                                                                                                                                                                                    >  1:614833/1‑38 (MQ=255)

                                                                                                                                                                                                                                                                            aCCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGc                                                                                                                                                                                                                                                                    <  2:614833/38‑1 (MQ=255)

                                                                                                                                                                                                                                                                            aCCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGCCAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCAGGTCAGAGTGATAAAAGGCGAATCAACAGgt                                                        >  2:853691/1‑242 (MQ=255)

                                                                                                                                                                                                                                                                            aCCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGCTAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGAGAAGGTGAGAATGTCATCGGTAt         >  1:579713/1‑289 (MQ=255)

                                                                                                                                                                                                                                                                            aCCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCcgcg                                                                               >  2:685369/1‑219 (MQ=255)

                                                                                                                                                                                                                                                                            aCCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTTATCAAAGTCGACTCAACAGGGGTTGGTCTGATTTTTGTTAATGGCGAAGGTCAGGATTGCATTGGTAt         >  2:133116/1‑289 (MQ=255)

                                                                                                                                                                                                                                                                            aCCACATTCTTTATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTAt         >  1:618919/1‑289 (MQ=255)

                                                                                                                                                                                                                                                                                   gcTTTATCTTGAATCTTTTCCTTCTCCAGGCGAAACCGTAACCGGTAACCACTTGCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGcc  <  2:877035/288‑1 (MQ=255)

TGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAGCATTAATGCTGACCACATTCTTAATCTTCAATCTTTTCCTACTCCAGGCGAAACCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCC  >  NC_000913/3937072‑3937635

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