BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	68,476	A→G	C525R (TGC→CGC)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	68,476	0	A	G	98.5%	393.3 / ‑2.9	129	C525R (TGC→CGC)	araB	ribulokinase
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (64/63); total (66/63)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

Only 100 of 129 total aligned reads displayed.

ATTAACGACGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGCTCGCTGCACGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGAGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGATCGGTAAAGCA  >  NC_000913/68230‑68719

aTTAACAACGTGCTAGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGt                                                                                                                                                                                                                                                  >  1:836179/1‑251 (MQ=255)

  tAACAACGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGTCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGtt                                                                                                                                                                                                                                                 <  1:687678/250‑1 (MQ=255)

    acaacgTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTCTAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGcg                                                                                                                                                                                                                                              >  2:339354/4‑250 (MQ=255)

      aacgTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCCTTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAgg                                                                                                                                                                                                                                            >  1:809632/2‑251 (MQ=255)

              cGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTtc                                                                                                                                                                                                                                      <  1:609424/249‑1 (MQ=255)

               ggCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGCAGTTGGAAGCTAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCGAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTctct                                                                                                                                                                                                                                   >  2:391813/1‑251 (MQ=255)

                 cATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTAc                                                                                                                                                                                                                                 >  2:369897/1‑251 (MQ=255)

                  aTGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTAcc                                                                                                                                                                                                                                <  1:815588/251‑1 (MQ=255)

                       gggTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCAc                                                                                                                                                                                                                            <  1:39335/250‑1 (MQ=255)

                        ggTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCACAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTg                                                                                                                                                                                                                          >  1:478470/1‑251 (MQ=255)

                         gTGACCTGACTCAGGGTTTCCGGGCCATACAGAGTCTGGCTGCCAATGCCAAACCACTCTTCATAATGATCAAAAATCGTCATGATCTTGGCCTTATAGAGTCGCAACGTCCTGGGCTGTCTGTGCCGGGTCGGAAGTGGGAAGTGAGTGTTGTTCGGCGCTCATCGCCCAGTGCTGCTAGCGGCGATAAAGCTGTTCAAAGCGTTGTTCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTgg                                                                                                                                                                                                                         <  2:267219/251‑1 (MQ=255)

                           gACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGcc                                                                                                                                                                                                                       >  2:724736/1‑251 (MQ=255)

                            aCCTGACGCCGGGTTTCCGGTCCCTACAGATGCTGGCTGCAAATGACACACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGCAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGGTAAAGCTGTTCAATGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTCACGCCCTGGccc                                                                                                                                                                                                                      >  2:547003/1‑250 (MQ=255)

                              cTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCAtt                                                                                                                                                                                                                    >  1:40008/1‑251 (MQ=255)

                                gACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTTTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTGCGGCGCTCAGCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCAtttt                                                                                                                                                                                                                  <  1:754/251‑1 (MQ=255)

                                  cGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCAtttttt                                                                                                                                                                                                                <  1:371447/251‑1 (MQ=255)

                                      ggtgttCCGTGCCATACAGATGCTGGATGCCAATGACAAACCTCACTGCATAATTATCAAAAATCGTTATTATCGTTTCCTTATAGAGGCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTgctg                                                                                                                                                                                                            <  2:584288/247‑1 (MQ=255)

                                        gTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCAGTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTgctggg                                                                                                                                                                                                          <  1:547003/251‑3 (MQ=255)

                                           tCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTgctgtgct                                                                                                                                                                                                        <  1:810334/250‑1 (MQ=255)

                                           gccGGACCATACAGTAGCCGGCTGCCAATGACAAACCACACTTCAAAATTATCAAAAATCGTCATTATAGAGTCCTTATAGATTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAGGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTTCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTTCAAGGCTTTCTCTACCGCACTGGCCATTTTTTgctgggct                                                                                                                                                                                                        <  2:673723/249‑1 (MQ=255)

                                             cGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTgctgggctgtt                                                                                                                                                                                                     >  1:733598/1‑249 (MQ=255)

                                                 ccATACAGATGCTGGCTGCCAAGGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGa                                                                                                                                                                                                 <  2:625363/251‑1 (MQ=255)

                                                    tACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGt                                                                                                                                                                                              >  1:535114/1‑251 (MQ=255)

                                                       aGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTg                                                                                                                                                                                           >  2:165524/1‑251 (MQ=255)

                                                       aGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTg                                                                                                                                                                                           <  1:980818/251‑1 (MQ=255)

                                                               gctgCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACt                                                                                                                                                                                   >  1:511731/1‑251 (MQ=255)

                                                                ctgCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCg                                                                                                                                                                                                                               >  2:962452/1‑206 (MQ=255)

                                                                ctgCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCg                                                                                                                                                                                                                               <  1:962452/206‑1 (MQ=255)

                                                                  gCCAATGACAAACCTCACTTCATAATTATCAAAAAGCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACttt                                                                                                                                                                                 <  1:534866/250‑1 (MQ=255)

                                                                   ccAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTcg                                                                                                                                                                               <  1:509562/251‑1 (MQ=255)

                                                                   ccAATGACAAACCAAACTTCATAATTCTCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGACTATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGAAAAAGCTGTTCACAGCGTTGTGCCTGTTCGCTGCGCGGGTGCAGGGTTTTCTCTACCGCAATGACCCTTTTTTGCTGGGCTGTTGGGATGGCTGCGTGCACTTTcg                                                                                                                                                                               >  2:391614/1‑251 (MQ=255)

                                                                       tGACAATCCACACTTCATAATTATCAAAAAACGTCATTACCATGTCCTTATAGAGTCGCAACAGCATGGGCAGTCTGTGACGGGGCGGAAGTTGGCAGATAGTGTTGTTCAGCGCTCATCGCCCCTTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTTCGCGGTTGCCGGGTTTTCTCTTCCGCACTGGCCATTTTTTGCTGAGCTACGGGGATGTCTTCGTGCACTTTcgacgc                                                                                                                                                                           >  2:623323/1‑247 (MQ=255)

                                                                        gACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCgcggcg                                                                                                                                                                          <  1:718186/251‑1 (MQ=255)

                                                                         aCAAACCTCACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGTACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTcgcg                                                                                                                                                                             <  2:510350/247‑1 (MQ=255)

                                                                          cAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCAGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAACGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTCCCGCACTGGCCATTTTTTGCTGGGCTGATGGGCTGTCTGCGTGCACTTTCGCGGCGAc                                                                                                                                                                        >  2:472739/1‑251 (MQ=255)

                                                                               cacaCTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGgc                                                                                                                                                                     >  1:116357/1‑249 (MQ=255)

                                                                                acacTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGgcagc                                                                                                                                                                  <  1:142865/251‑1 (MQ=255)

                                                                                      cATAATGATCAAAAAGCGTCATTATCGTGTCCTTAAAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAATTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGGGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCaaaaa                                                                                                                                                             <  2:622081/250‑1 (MQ=255)

                                                                                       aTAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTCCACTTTCGCGGCGACGGCAGCAAAAATc                                                                                                                                                           >  1:963955/1‑251 (MQ=255)

                                                                                       aTAATTATCAAAAATCGTCAGTATCGTGTCCTTAAAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGAGGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATc                                                                                                                                                           <  1:312941/251‑1 (MQ=255)

                                                                                       aTAATTATCAAAAAACGTCATTATCTTGTCCTTATAGAGTCGCAACGGCCTGGGCATCCTGTGCCTGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCGCATCGCCCATTGCTGAGAGCGGCGATAAAGCTGTTCAAAGCTTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGGTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATc                                                                                                                                                           <  1:121378/251‑1 (MQ=255)

                                                                                            taacAAAAATCGTCATTTTCGTGTCCTTATAGAGGCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGGTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCGCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATcgacg                                                                                                                                                       <  2:66451/247‑4 (MQ=255)

                                                                                             ttcaaaaagtgacATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCGGATAGCGGCGATAAAGCTGGTCAAAGCGTTGTGCCTGTTCGCTGCGCGTTGGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGATGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCTACGGCAGCACAAATcgccgc                                                                                                                                                      <  2:425478/238‑1 (MQ=255)

                                                                                               cAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCAcc                                                                                                                                                   >  1:586918/1‑251 (MQ=255)

                                                                                               aaaaaaGCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGAACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCAcc                                                                                                                                                   <  1:994189/250‑1 (MQ=255)

                                                                                                aaaaaTCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGAGAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCg                                                                                                                                                  <  2:58366/251‑1 (MQ=255)

                                                                                                 aaaaTCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGa                                                                                                                                                 <  1:496510/251‑1 (MQ=255)

                                                                                                  aaaTCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATACAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGATTTCCCTACCGCACTGGCCATTTTTTGCTGGGCTGTTGGGATGTATGCGTGCAATTTCGCGGCGACGGCAGCAGAAATCGCCGCAcctgg                                                                                                                                                >  2:784809/1‑248 (MQ=255)

                                                                                                   aaTCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAgc                                                                                                                                               >  1:440756/1‑251 (MQ=255)

                                                                                                    aTCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAgcg                                                                                                                                              >  1:898294/1‑251 (MQ=255)

                                                                                                          aTTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCCGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGCTGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGc                                                                                                                                        >  1:553825/1‑251 (MQ=255)

                                                                                                          aTTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCCTTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCCACGGCAGACAAAATCGCCGCACCGAgcgaccagg                                                                                                                                        >  2:921549/1‑245 (MQ=255)

                                                                                                          aTTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGTCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGc                                                                                                                                        >  1:163609/1‑251 (MQ=255)

                                                                                                            tATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCCTTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCCCTGCGACGGCAGCAAAAAACGACGCACCGAgggccccgcac                                                                                                                                      >  2:419241/1‑241 (MQ=255)

                                                                                                              tCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTg                                                                                                                                    >  1:857365/1‑251 (MQ=255)

                                                                                                                 tgtCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTc                                                                                                                                 >  1:143746/1‑251 (MQ=255)

                                                                                                                  gcccTTTTAGAGCCGTAACGGCCTGGGCAGCCTGGTCCGGGGGGGAAGTTGGCAGATATTGGTGTGCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTc                                                                                                                                 <  2:591003/248‑1 (MQ=255)

                                                                                                                        ataGAGTCGCAACTGCCTGGGCTGCCTGTGCCTGGGGGGAAGTGGGAAGATAGTGTTGCGCGGCGCTCATCGCCCATTGCTGATGGCGGCGATAAGGCGGTTCAAAGCGTTGTTCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTGCCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTGCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACGGGTCagaagt                                                                                                                           <  2:217323/250‑4 (MQ=255)

                                                                                                                          agagTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACCGCACTGGTCAGAGGcaac                                                                                                                        >  1:219912/1‑251 (MQ=255)

                                                                                                                            agTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCCCAGCACTGGTCAGAGGcaaca                                                                                                                       >  2:291223/1‑250 (MQ=255)

                                                                                                                             ggcGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTACCTGTTCGCTACGCGGTTGCAGGGCTTTCTCTACCGAACTGGCCATTTTTTGCTGGGCTGATGGGATGGCTGCGTCCACTTTAGCGGCGACGGCAGCAACAATCGCCGCACCGAGCGCACAGCAcacgtcaacagcaaccat                                                                                                                     >  2:472000/3‑234 (MQ=255)

                                                                                                                                    gggactgggaggtgtgtgCCGGTGGGGAAGGTTGAAGAGAGTGTTATTCTGCGCTCAGCGTCCTTTGTGGATAGCGGCGATAAAGCTGGTCTAAGCGTTGTGCCTTTTCGCTGCGCGGTTGCAGTGTTTTCTCTACCGCACGGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAAACGCTGCACCGAGCGCACAGCACTGTTCAAAGTCAACAATTTG‑CAGc                                                                                                             <  2:117890/237‑1 (MQ=255)

                                                                                                                                    gggCCGGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATTGTTTTGTTCGGCGCGCACCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTTTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGc                                                                                                             <  2:897347/250‑1 (MQ=255)

                                                                                                                                     ggCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCCCAACGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCg                                                                                                            >  2:893263/1‑251 (MQ=255)

                                                                                                                                                              gAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGAGATTCAGCACGTCGCAGCAGGc                                                                                   >  1:221061/1‑251 (MQ=255)

                                                                                                                                                                         aTAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGa                                                                        <  1:286162/251‑1 (MQ=255)

                                                                                                                                                                         aTAGTGTTGCTCGGCGCTCAGCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGa                                                                        <  1:212614/251‑1 (MQ=255)

                                                                                                                                                                              gttgttCGGCGCTCATCGCCCATTGCTGATAGCGGCGATACAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGGTTCAGCACGTCGCAGCAGGCCTGCATAATGACCTgg                                                                   >  1:318521/1‑251 (MQ=255)

                                                                                                                                                                              gttgttCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTg                                                                    <  1:37853/250‑1 (MQ=255)

                                                                                                                                                                              gtgtgtcgtcgcTCATCGCCTAGTGCTGATAGCGGCGAGAAATCGGGTTAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGTGGTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGTGATGTCTGCGGGCACGTTCGCGGCTACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGATGCAACATGTTG‑CAGCGGGCGATTCAGCACGTCTCAGCAGGCCTGCATAATGAccagg                                                                   <  2:267773/242‑4 (MQ=255)

                                                                                                                                                                               ttgttCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGt                                                                  >  2:962990/1‑251 (MQ=255)

                                                                                                                                                                                tgttCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGtt                                                                 >  1:733407/1‑251 (MQ=255)

                                                                                                                                                                                     ggCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTcc                                                            <  1:925292/251‑1 (MQ=255)

                                                                                                                                                                                        gcgcATCGCCCATTGCTGAGAGCGGCGATAAAGCTGGTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCCCGCCCGCACTGGCCATTTTGTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTTAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCcgcg                                                         <  2:400963/248‑1 (MQ=255)

                                                                                                                                                                                          acATCGCACATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGATGCGCGGTTGAAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCAACGAGCGCACAGCACTGGTCAGAGGCCACAATTTGACA‑CGGGCGATTCAGCACGACGCAACAGGCCAGCATAATGACCTGGTTTTTCcgcgcg                                                       >  2:331422/2‑251 (MQ=255)

                                                                                                                                                                                              cGCCTATTGCTGATAGCGGGGATAAAGCTGTCCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCGGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATgc                                                   <  2:716995/251‑1 (MQ=255)

                                                                                                                                                                                               gCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATgcc                                                  >  1:776241/1‑251 (MQ=255)

                                                                                                                                                                                                 ccATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGCGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCCGGCCTGCATAATGACCTGGTTTTTCCGCTCGATgccgc                                                >  2:428799/1‑251 (MQ=255)

                                                                                                                                                                                                 ccATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATgccgc                                                >  1:442508/1‑251 (MQ=255)

                                                                                                                                                                                                     tGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCCGGCCTGCATAATGACCTGGGTTTTCAGCGCGATGcccccccg                                            >  2:966892/1‑249 (MQ=255)

                                                                                                                                                                                                     tGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGccc                                              >  1:450640/1‑249 (MQ=255)

                                                                                                                                                                                                         gATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGcc                                        >  1:240572/1‑251 (MQ=255)

                                                                                                                                                                                                          aTAGCGGCGATCAAGCGGTTCAAAGCGTTGTGCCTGTACGCTGAGCGGTTGCAGGGTTTTCTCTAACGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCAGTTCAGCACGTCGCAGCAGGCCTGCTTAGTGACATGGTTTTTCCACGCGATgccgccacgagccc                                       >  2:231829/1‑243 (MQ=255)

                                                                                                                                                                                                          aTAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCa                                       <  1:221807/251‑1 (MQ=255)

                                                                                                                                                                                                               ggcgATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGCTGCCGCCCAGCGCCCTCACg                                  >  1:174677/1‑251 (MQ=255)

                                                                                                                                                                                                                  gATAAAGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGtta                               <  1:962990/251‑1 (MQ=255)

                                                                                                                                                                                                                     aaaGCTGTTCAAAGCGTTGTGCCTGTTCGCGGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCAGTTTTTGCTGGGCTGATGGGAGGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCCCCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGAGATTCAGCACGTCGCAGCAGTCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGttattg                            <  1:48343/251‑2 (MQ=255)

                                                                                                                                                                                                                     aaaGCTGTTCAAAGCGTTGGGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCACACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCCGCACGTCGCAGCAGGCCTGAAAAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCCACACGttttt                             >  2:652197/1‑247 (MQ=255)

                                                                                                                                                                                                                      aaGCTGTTCAAAGCGTTGTGCCTTTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGttattga                           <  2:768888/251‑3 (MQ=255)

                                                                                                                                                                                                                       aGCTGTTCAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCCGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGttattgac                          >  1:361901/1‑248 (MQ=255)

                                                                                                                                                                                                                              cAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGAt                    <  2:558497/250‑1 (MQ=255)

                                                                                                                                                                                                                              cAAAGCGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATc                   <  1:724736/251‑1 (MQ=255)

                                                                                                                                                                                                                                 aGCGTTGCGCCTGTTCGCTGCGCGGTTGCAGGGTTTGCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGTCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGCCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATcccc                <  1:530425/251‑1 (MQ=255)

                                                                                                                                                                                                                                  gcgttggtccTGTTCGCTGCGCGGTTGCAGGGGTTTCGCTTCCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCTACGGCAGCAAAAATCGCCTCACCGAGCTCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGGCGCAGCAGTCCTGCATAATGACCTTGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCt               <  1:887335/243‑1 (MQ=255)

                                                                                                                                                                                                                                   cGTTGTGCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTg              >  2:326062/1‑251 (MQ=255)

                                                                                                                                                                                                                                       gtgCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGCTGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCg          >  2:242556/1‑251 (MQ=255)

                                                                                                                                                                                                                                       gtgCCTGTTCGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATc           <  1:291223/250‑1 (MQ=255)

                                                                                                                                                                                                                                            tgttCGCTGCGAGGTTGCAGGGTTTTCTCTAACGCACTGGCCATTTTTTGCTGGGCTGGTGGGATGTATGCGTGCACTTTCGCGGCGACGGCAGCAACAATCGCCGCACCGAGCGCACAGCAATGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGCCGCAGCAGGCCTGCATAAAGACCTGGGTTTTCCGCGCGATGCGGCCCAGCGCCATCACGTTATTGACCGCGATCCCCTGCTCggcac      >  2:122442/4‑247 (MQ=255)

                                                                                                                                                                                                                                              cTCGCTGCACGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGAGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGATCGGTAAAGc   >  1:286342/1‑251 (MQ=255)

                                                                                                                                                                                                                                                cGCTGCGCGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGGGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGCGCCATCACGTTATTGACGGCGATCCCCTGATCGGTAAAGCa  <  1:770721/250‑1 (MQ=255)

ATTAACGACGTGCTCGGCATGTTGGGTGACCTGACGCAGGGTTTCCGGGCCATACAGATGCTGGCTGCCAATGACAAACCACACTTCATAATTATCAAAAATCGTCATTATCGTGTCCTTATAGAGTCGCAACGGCCTGGGCAGCCTGTGCCGGGGCGGAAGTTGGAAGATAGTGTTGTTCGGCGCTCATCGCCCATTGCTGATAGCGGCGATAAAGCTGTTCAAAGCGTTGTGCCTGCTCGCTGCACGGTTGCAGGGTTTTCTCTACCGCACTGGCCATTTTTTGCTGAGCTGATGGGATGTCTGCGTGCACTTTCGCGGCGACGGCAGCAAAAATCGCCGCACCGAGCGCACAGCACTGGTCAGAGGCAACAATTTG‑CAGCGGGCGATTCAGCACGTCGCAGCAGGCCTGCATAATGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGATCGGTAAAGCA  >  NC_000913/68230‑68719

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