BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,269,785	G→A	3.9%	T9I (ACT→ATT)	ldrB ←	toxic polypeptide, small

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,269,785	0	G	A	3.9%	369.4 / 2.1	116	T9I (ACT→ATT)	ldrB	toxic polypeptide, small
Reads supporting (aligned to +/- strand): ref base G (49/61); new base A (3/2); total (53/63)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 116 total aligned reads displayed.

TTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGCTTGGTGTGGAT‑ACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCC  >  NC_000913/1269506‑1270062

tttGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCaa                                                                                                                                                                                                                                                                                   >  1:88546/1‑289 (MQ=255)

tttGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCaa                                                                                                                                                                                                                                                                                   >  2:21165/1‑289 (MQ=255)

  tGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCaaa                                                                                                                                                                                                                                                                                  >  2:573845/1‑288 (MQ=255)

       tgaatggatcATTCCGGTTTTTCAGCAGCGGTTAATGTCATTAAATGAGTGCTCTGGATGAGGCTGCAGGATCAGCTTGCCAAAATTTACTTGTGGAAGGGTTAAAGCCAAAAAGGGCCGCATGCTGTTACCCCGGCGGCCGGTTTCGGTCGCATATTGCCATCACGGAAGCCTGACGCCCGTTTTAACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTgcgc                                                                                                                                                                                                                                                                            <  1:713316/280‑1 (MQ=11)

       tGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTgcgc                                                                                                                                                                                                                                                                            <  1:849938/289‑1 (MQ=255)

       tGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTgcgc                                                                                                                                                                                                                                                                            <  2:427612/289‑1 (MQ=255)

              tGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTCCTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCCCGGGCGt                                                                                                                                                                                                                                                                      >  2:409411/1‑288 (MQ=255)

                   tCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCAtata                                                                                                                                                                                                                                                                >  1:74575/1‑289 (MQ=39)

                   tCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCCGTGCCGCCAGGGCGTTCCAGAAAGTCATGGCAAACTGCGCGAGCGTctacta                                                                                                                                                                                                                                                                >  1:260820/1‑284 (MQ=255)

                      ggcttttCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAgc                                                                                                                                                                                                                                                              >  2:57409/4‑288 (MQ=1)

                        tttttCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCCGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTAACATAgcgg                                                                                                                                                                                                                                                            >  1:479336/1‑286 (MQ=9)

                        tttttCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAgccg                                                                                                                                                                                                                                                            >  1:132112/1‑288 (MQ=18)

                                     gTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGa                                                                                                                                                                                                                                              <  1:688368/289‑1 (MQ=38)

                                           gTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACAt                                                                                                                                                                                                                                        >  2:598576/1‑289 (MQ=37)

                                             cATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACAttt                                                                                                                                                                                                                                      >  1:31349/1‑289 (MQ=37)

                                                ttCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGc                                                                                                                                                                                                                                   >  1:101858/1‑289 (MQ=35)

                                                   aTGAGTGCTTTGGATGAGGGTTCCAGCGCGGGTTGCCAATATTTACTTTTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGc                                                                                                                                                                                                                                <  1:436083/289‑1 (MQ=25)

                                                   aTGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGc                                                                                                                                                                                                                                <  2:31349/289‑1 (MQ=35)

                                                   aTGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGc                                                                                                                                                                                                                                <  1:715466/289‑1 (MQ=35)

                                                   aTGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGc                                                                                                                                                                                                                                <  2:603999/289‑1 (MQ=35)

                                                    tGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGc                                                                                                                                                                                                                                <  1:554732/288‑1 (MQ=35)

                                                        tGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAATGACAAAAATGGCCGCCGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTTCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACGTCCGGTTACGCCACCAGCCGACCATCGCTGCGGTAATAATTCCCGCCAGGATCGGTCCCGCCAGGTCGTGCCAGCAAGTCATGGCAAACTCCGCCAGCGTCACATAGCCGCCTTGTTGTAATGCCACCATTTAGCGGCTAttt                                                                                                                                                                                                                           >  1:811239/1‑288 (MQ=11)

                                                        tGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATCATTCCTGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATTTAGCCGCCTTGTTGTAATGCAAACATTTTTCGGCTAt                                                                                                                                                                                                                             >  1:316701/1‑287 (MQ=255)

                                                        tGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTAt                                                                                                                                                                                                                             >  2:317082/1‑287 (MQ=35)

                                                        tGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTAt                                                                                                                                                                                                                             >  2:554260/1‑287 (MQ=18)

                                                         gCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTc                                                                                                                                                                                                                           <  1:57409/288‑1 (MQ=18)

                                                          ctctGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCt‑t                                                                                                                                                                                                                        >  1:155779/1‑289 (MQ=14)

                                                          ctctGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCCGACGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTTCCAGAAAGTCATGGCAAACCGCGCGAACGTCCTAGAGCCGACTTGTTGTAATGACAACAATTTCAGGCGATTCt‑t                                                                                                                                                                                                                        >  1:215458/1‑289 (MQ=11)

                                                          ctctGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGc                                                                                                                                                                                                                                <  2:7029/282‑1 (MQ=35)

                                                          ctctGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGGCGTGCAAGAAAGTCCTGGCAAACTGCGCGAGTGCCATATAGCCGCCTTGTTTTAATGTCAACATTTTCTGGCTcttct‑t                                                                                                                                                                                                                        >  1:433119/1‑289 (MQ=11)

                                                                   aGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTc                                                                                                                                                                                                               <  2:74575/288‑1 (MQ=35)

                                                                    ggCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTg                                                                                                                                                                                                             >  2:673929/1‑289 (MQ=34)

                                                                    ggCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCCTCCTATAGCCGCCTTGTTCTAATGCCCAACTTTTGCGGCTATTCT‑TGAAGCTGTCt                                                                                                                                                                                                              >  1:7029/1‑289 (MQ=255)

                                                                    ggCTGCCAGCCGGGGTAGCCAAAATTTACTGTTGGAAGGGAAAAAGCAAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCt                                                                                                                                                                                                              <  1:573845/288‑1 (MQ=11)

                                                                     gCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTTCCAGAAAGTCATGGCAAACTGCCCGAACGTCATATAGCCGCCTTGTTGTAATGACAACCTTTTGCGCCTATTAC‑TGAAG‑TGTCTgg                                                                                                                                                                                                            >  1:489686/1‑289 (MQ=255)

                                                                               cGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGcgcg                                                                                                                                                                                                                                                                           >  2:11044/1‑218 (MQ=35)

                                                                               cGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGcgcg                                                                                                                                                                                                                                                                           <  1:11044/218‑1 (MQ=35)

                                                                                     gCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGccc                                                                                                                                                                                             >  2:714677/1‑288 (MQ=34)

                                                                                         atatTTACTTGTGGAAGTGATAAATACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCTGTTTCGGGCGCATATTGCCATCACGGCAGCCTAACGCCCGTTTTCACCTTACGTCCGGTTACGCAACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGtt                                                                                                                                                                                       <  2:853132/289‑1 (MQ=19)

                                                                                                cTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCtgttgt                                                                                                                                                                                <  2:406779/289‑1 (MQ=35)

                                                                                                cTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCtgttgt                                                                                                                                                                                <  1:598576/289‑1 (MQ=35)

                                                                                                    tgGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCTTTCGGCTGGTTTCGGCCGCATCTTGCCACCACGGCAGCCTGACGCACCTTTGCACCTTACTTCCGGTTACGCCACCACCCGACAATCGCTGCGGTAATAATTCCCGCCAGAACCGGTGCCGCCAGGTCGTGCCAAAAAGTCATGGCAAACTGCGCGAGCGCCATATAGCCCCCTTGTTGTACTCACCACATTTTTCCGCTACTCT‑TGAAG‑TTTCAGGTCTCAACATTAGCCCCACGTTCGTTGGTCAg                                                                                                                                                                             >  1:468381/1‑289 (MQ=255)

                                                                                                    tgGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAgg                                                                                                                                                                            >  2:294476/1‑289 (MQ=35)

                                                                                                        aGTGATAAAGACAAAAATGGCCCCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGtttt                                                                                                                                                                        <  2:572049/289‑1 (MQ=25)

                                                                                                         gTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGtttt                                                                                                                                                                        >  1:762337/1‑288 (MQ=35)

                                                                                                           gATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTAc                                                                                                                                                                      <  2:88546/288‑1 (MQ=34)

                                                                                                            aTAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTAcc                                                                                                                                                                     >  2:16428/1‑288 (MQ=34)

                                                                                                                aGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCa                                                                                                                                                                <  1:294476/289‑1 (MQ=34)

                                                                                                                aGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCa                                                                                                                                                                <  2:540928/289‑1 (MQ=34)

                                                                                                                   cAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACg                                                                                                                                                             <  2:433119/289‑1 (MQ=34)

                                                                                                                      aaaTGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCCCGGAAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGACGGGACGGGTGCCGCCAGGTCGTGCCAGAAATTCAGGGCAACCGGCGCGCGCGGCATATAGCCGCCTGTTTGTAATGAACACTTTTTGCGGC‑CTTCTCTGACT‑TGTCTGGTTTTTAGGTTACCCCC‑AGTTCTGTTGTCGGGttttccttatccagttgt                                                                                                                                                          >  1:636703/1‑275 (MQ=11)

                                                                                                                       aaTGCCCCCAGGCTGTTACCCCTGCGGCGGGTTTCGGGCACCTATTGCCATCACGGCAACCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCg                                                                                                                                                         <  1:409411/289‑1 (MQ=255)

                                                                                                                           gCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGAtt                                                                                                                                                                                                  <  1:519840/245‑1 (MQ=33)

                                                                                                                                            cTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTc                                                                                                                                                                                                               >  1:376169/1‑215 (MQ=34)

                                                                                                                                            cTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCa                                                                                                                                     >  2:637021/1‑288 (MQ=34)

                                                                                                                                             tGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTa                                                                                                                                                                                                                              <  1:352381/201‑1 (MQ=34)

                                                                                                                                              gCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGccc                                                                                                                                                                                             >  2:719177/1‑231 (MQ=33)

                                                                                                                                              gCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGccc                                                                                                                                                                                             <  1:719177/231‑1 (MQ=33)

                                                                                                                                              gCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                                                                                                                                  >  1:781838/1‑289 (MQ=34)

                                                                                                                                              gCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                                                                                                                                  >  2:164647/1‑289 (MQ=34)

                                                                                                                                              gCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAAcc                                                                                                                                  >  2:734537/1‑289 (MQ=34)

                                                                                                                                                ggCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCaa                                                                                                                                <  2:591229/289‑1 (MQ=34)

                                                                                                                                                   cGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATg                                                                                                                              <  2:215458/288‑1 (MQ=34)

                                                                                                                                                      tttCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcca                                                                                                                           >  2:217596/1‑288 (MQ=34)

                                                                                                                                                        tCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccacca                                                                                                                        <  2:883307/289‑1 (MQ=34)

                                                                                                                                                                atatTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATaaa                                                                                                                <  1:317082/289‑1 (MQ=34)

                                                                                                                                                                           aCGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAAc                                                                                                     >  2:158886/1‑289 (MQ=37)

                                                                                                                                                                                 gCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑gcg                                                                                              <  1:164647/289‑1 (MQ=37)

                                                                                                                                                                                  ccggACGCCCGGTTTCCCCTTACTTCCGGTTACGCCACCAGCCGCCAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑gcgc                                                                                             <  1:637021/286‑1 (MQ=255)

                                                                                                                                                                                   cTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑gcg                                                                                              <  1:714677/287‑1 (MQ=37)

                                                                                                                                                                                   cTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑gcg                                                                                              <  2:781838/287‑1 (MQ=37)

                                                                                                                                                                                                tCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATaacg                                                                               <  2:316701/289‑3 (MQ=25)

                                                                                                                                                                                                     ttACTTCCGGTTACGCCACCAGCCTACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGc                                                                          <  2:636703/289‑1 (MQ=39)

                                                                                                                                                                                                       aCTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGCtt                                                                        >  2:110918/1‑289 (MQ=39)

                                                                                                                                                                                                       aCTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGCtt                                                                        <  1:21165/289‑1 (MQ=39)

                                                                                                                                                                                                       aCTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGACCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCAGCAACATT‑GCGCCTCCCCGGATAATTCCGGCtt                                                                        >  2:596754/1‑289 (MQ=255)

                                                                                                                                                                                                          tCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGCTTGgt                                                                     <  1:673929/289‑1 (MQ=39)

                                                                                                                                                                                                          tCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGCTTGgt                                                                     <  2:811239/289‑1 (MQ=39)

                                                                                                                                                                                                          tCCGGAGACGCCACCAGCCGACAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCAAGCGTCATATAGCCGCCTTGTTGTAATGACAACATGTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGTGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGCTTGgt                                                                     <  1:596754/289‑1 (MQ=255)

                                                                                                                                                                                                              gTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGtt                                                                                                                                                                                       >  1:113861/1‑172 (MQ=35)

                                                                                                                                                                                                              gTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGtt                                                                                                                                                                                       <  2:113861/172‑1 (MQ=35)

                                                                                                                                                                                                                    ccaccaGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGCTTGGTGTGGAt‑acta                                                          <  1:16428/289‑1 (MQ=39)

                                                                                                                                                                                                                               cAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGCTTGGTGTGGAT‑ACTACGTCTCAAtt                                                <  2:101858/288‑1 (MQ=39)

                                                                                                                                                                                                                               cAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGC‑ACATTGGCGCCCCCCCGATTAAGTCCGGC‑TGGTGGGGATCACTCCGTCACATTTc                                               >  1:15949/1‑289 (MQ=11)

                                                                                                                                                                                                                                  tCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCa                                                                                         <  1:415725/245‑1 (MQ=37)

                                                                                                                                                                                                                                  tCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCa                                                                                         >  2:415725/1‑245 (MQ=37)

                                                                                                                                                                                                                                     cTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccacca                                                                                                                        <  2:354065/212‑1 (MQ=34)

                                                                                                                                                                                                                                     cTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGccacca                                                                                                                        >  1:354065/1‑212 (MQ=34)

                                                                                                                                                                                                                                          gTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGccc                                                                                                                                                                                             <  1:13770/139‑1 (MQ=34)

                                                                                                                                                                                                                                          gTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGccc                                                                                                                                                                                             >  2:13770/1‑139 (MQ=34)

                                                                                                                                                                                                                                          gTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCg                                                                                                                                                         >  2:280787/1‑174 (MQ=35)

                                                                                                                                                                                                                                                     cGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGc                                                                          >  2:62341/1‑241 (MQ=39)

                                                                                                                                                                                                                                                     cGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGc                                                                          <  1:62341/241‑1 (MQ=39)

                                                                                                                                                                                                                                                        cAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGCTTGGTGTGGAT‑ACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGt                       <  1:110918/288‑1 (MQ=39)

                                                                                                                                                                                                                                                                gTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGc                                                                          <  1:839436/230‑1 (MQ=40)

                                                                                                                                                                                                                                                                gTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGc                                                                          >  2:839436/1‑230 (MQ=40)

                                                                                                                                                                                                                                                                      ccAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTGATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGCTTGGTGTGGAT‑ACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAagag          >  2:471510/1‑287 (MQ=39)

                                                                                                                                                                                                                                                                           tCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGc                                                                          >  2:550881/1‑219 (MQ=39)

                                                                                                                                                                                                                                                                           tCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGc                                                                          <  1:550881/219‑1 (MQ=39)

                                                                                                                                                                                                                                                                            cGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGCTTGGTGTGGAT‑ACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTcc  <  2:103041/289‑1 (MQ=39)

TTTGCAATGAAGAGTGCATTCCGGTTTTTCAACAGCTGTTACAGTCATTTCATGAGTGCTCTGGATGAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCT‑TGAAG‑TGTCTGGTTTCAAGATTAGCCCC‑CGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATT‑GCGCCTCACCGGATAATTCCGGCTTGGTGTGGAT‑ACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCC  >  NC_000913/1269506‑1270062

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