BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_008463	3,670,543	+A	coding (543/1020 nt)	PA14_RS16675 →	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_008463	3,670,542	1	.	A	99.3%	1749.6 / 2.4	418	coding (542/1020 nt)	PA14_RS16675	ABC transporter permease
Reads supporting (aligned to +/- strand): ref base . (3/0); new base A (198/215); total (201/217)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 428 total aligned reads displayed.

ACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTTCCTGG‑AGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCCGCGG  >  NC_008463/3670392‑3670686

aCGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACACTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACGACGTCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTgg                                                                                                                                                   >  2:5538957/1‑150 (MQ=255)

   tGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGG‑aag                                                                                                                                               >  1:7157514/1‑148 (MQ=255)

    gCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGG‑aagt                                                                                                                                              >  2:9687450/1‑147 (MQ=255)

         cccGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGt                                                                                                                                          <  2:9479792/150‑1 (MQ=255)

            gggTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCgg                                                                                                                                       >  2:2905338/1‑150 (MQ=255)

                    tACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCg                                                                                                                               <  2:8427196/150‑1 (MQ=255)

                    tACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCg                                                                                                                               >  2:1161140/1‑150 (MQ=255)

                    tACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCg                                                                                                                               >  1:8001460/1‑150 (MQ=255)

                       ggCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTg                                                                                                                            >  1:4322405/1‑150 (MQ=255)

                        gCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGc                                                                                                                           >  2:10603559/1‑150 (MQ=255)

                        gCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGc                                                                                                                           >  1:5588049/1‑150 (MQ=255)

                            ccGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTa                                                                                                                       <  1:10665066/150‑1 (MQ=255)

                             cgcgTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGATCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGACGGGGCTGCCGGTGCTCTAc                                                                                                                      >  1:9861615/1‑150 (MQ=37)

                               cgTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACct                                                                                                                    <  1:7404170/150‑1 (MQ=255)

                                gtgtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctg                                                                                                                   >  1:6753265/1‑150 (MQ=255)

                                 tgtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgc                                                                                                                  <  2:1085807/150‑1 (MQ=255)

                                  gtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgct                                                                                                                 >  1:8612213/1‑150 (MQ=255)

                                  gtCGATCAACATCGCCCTGACCCTGACCGTGCTAAGATATCTGATCGGCGTCACCGTCGGCGACCAGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGACGGGGCTGCCTGAGCTCTACctgct                                                                                                                 >  2:8620740/1‑150 (MQ=255)

                                    cGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGa                                                                                                               <  1:7492169/150‑1 (MQ=255)

                                     gATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGat                                                                                                              <  1:9611677/150‑1 (MQ=255)

                                       tCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatca                                                                                                            >  1:7642786/1‑150 (MQ=255)

                                        ccTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatcat                                                                                                           <  2:8723488/150‑1 (MQ=255)

                                           ctTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCt                                                                                                        >  2:6927559/1‑150 (MQ=255)

                                                ccctgaccctgaCCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCa                                                                                                   >  1:4889720/1‑150 (MQ=255)

                                                ccctgaccctgaCCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGATGGGTCGACCTGATCGGCCCGCG‑TTTCCTGGAAGTCTCGTCGGGTCTGCCGGTGCTCTACCTGCTGATCATCATGGACa                                                                                                   >  2:4118905/1‑150 (MQ=255)

                                                   tgaccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCt                                                                                                <  1:1952296/150‑1 (MQ=255)

                                                      ccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCg                                                                                             >  1:9638679/1‑150 (MQ=255)

                                                        ctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTc                                                                                           >  1:3941795/1‑150 (MQ=255)

                                                         tgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTcc                                                                                          >  2:938216/1‑150 (MQ=255)

                                                           accGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAg                                                                                        >  1:1466608/1‑150 (MQ=255)

                                                             cGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGcc                                                                                      <  2:1584588/150‑1 (MQ=255)

                                                              gTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCg                                                                                     >  2:2403353/1‑150 (MQ=255)

                                                                gCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGaa                                                                                   >  1:8996045/1‑150 (MQ=255)

                                                                 ctgaGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAAc                                                                                  <  1:3056542/147‑1 (MQ=255)

                                                                     gTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCt                                                                              <  1:5880653/150‑1 (MQ=255)

                                                                      ttCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCtg                                                                             <  1:1656843/150‑1 (MQ=255)

                                                                        ctctGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCtggt                                                                           <  2:5236003/150‑1 (MQ=255)

                                                                           tGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgc                                                                        >  1:4005043/1‑150 (MQ=255)

                                                                              tCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgctgc                                                                     >  2:5715682/1‑150 (MQ=255)

                                                                                ggCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGctgctg                                                                   <  1:2500417/150‑1 (MQ=255)

                                                                                 gCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTgg                                                                  <  2:1848048/150‑1 (MQ=255)

                                                                                 gCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTgg                                                                  >  1:7022240/1‑150 (MQ=255)

                                                                                  cGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTggg                                                                 <  1:8220172/150‑1 (MQ=255)

                                                                                   gTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGc                                                                >  1:9443804/1‑150 (MQ=255)

                                                                                    tCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGca                                                               >  2:743799/1‑150 (MQ=255)

                                                                                     cACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcat                                                              >  1:6154141/1‑150 (MQ=255)

                                                                                       ccGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcatca                                                            >  1:5300442/1‑150 (MQ=255)

                                                                                       ccGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcatca                                                            <  2:6670994/150‑1 (MQ=255)

                                                                                         gTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtg                                                          >  2:10785148/1‑150 (MQ=255)

                                                                                         gTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtg                                                          >  1:6658797/1‑150 (MQ=255)

                                                                                          tCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgc                                                         <  1:10333073/150‑1 (MQ=255)

                                                                                          tCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgc                                                         <  2:340710/150‑1 (MQ=255)

                                                                                           cGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgct                                                        >  2:5636144/1‑150 (MQ=255)

                                                                                           cGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgct                                                        <  1:2680888/150‑1 (MQ=255)

                                                                                             gcgcCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgctgc                                                      <  2:2090331/150‑1 (MQ=255)

                                                                                                 ccTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGtt                                                  >  1:5684775/1‑150 (MQ=255)

                                                                                                  cTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGttt                                                 >  2:5852481/1‑150 (MQ=255)

                                                                                                   tGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGtttt                                                >  1:3750178/1‑150 (MQ=255)

                                                                                                   tGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGtttt                                                <  2:3695680/150‑1 (MQ=255)

                                                                                                   tGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGtttt                                                <  2:3071370/150‑1 (MQ=255)

                                                                                                   tGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGtttt                                                <  2:223468/150‑1 (MQ=255)

                                                                                                     ctaGGCTACGTCGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGTAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATGAGGCTGCTGTTTTcc                                              <  2:7284821/148‑1 (MQ=255)

                                                                                                     cAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTcc                                              <  2:5427781/150‑1 (MQ=255)

                                                                                                     cAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTcc                                              >  1:2319733/1‑150 (MQ=255)

                                                                                                      aaGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCt                                             <  2:7157514/150‑1 (MQ=255)

                                                                                                      aaGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCt                                             <  1:3231917/150‑1 (MQ=255)

                                                                                                      aaGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCt                                             >  2:2726320/1‑150 (MQ=255)

                                                                                                        ggCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTgg                                           >  2:4406184/1‑150 (MQ=255)

                                                                                                        ggCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTgg                                           >  2:10027252/1‑150 (MQ=255)

                                                                                                        ggCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTgg                                           <  1:1189366/150‑1 (MQ=255)

                                                                                                         gCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGa                                          <  2:164447/150‑1 (MQ=255)

                                                                                                           tactacGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATg                                        <  1:840467/150‑1 (MQ=255)

                                                                                                           tactacGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATg                                        <  2:10974564/150‑1 (MQ=255)

                                                                                                            actacGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGa                                       >  2:2540598/1‑150 (MQ=255)

                                                                                                            actacGGCGGCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGa                                       <  1:5048826/150‑1 (MQ=255)

                                                                                                                cggcggCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCt                                   <  2:1674469/150‑1 (MQ=255)

                                                                                                                cggcggCTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCt                                   >  1:10019833/1‑150 (MQ=255)

                                                                                                                 ggcggcTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCtg                                  >  2:8023973/1‑150 (MQ=255)

                                                                                                                  gcggcTGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCtgg                                 <  2:822589/150‑1 (MQ=255)

                                                                                                                     gcTGGGTCGACCTGCTCGGCCAGCG‑TGTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGGGGCTGCTGGGCCTCATGCTGCTGTTTTCCGGGATGAGCCgggtgg                              >  2:109010/1‑150 (MQ=255)

                                                                                                                       tGGGTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGAc                            <  1:346862/150‑1 (MQ=255)

                                                                                                                        gggTCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACg                           >  2:8798053/1‑150 (MQ=255)

                                                                                                                           tCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTgg                        >  1:8699415/1‑150 (MQ=255)

                                                                                                                           tCGACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTgg                        <  2:8978392/150‑1 (MQ=255)

                                                                                                                             gACCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTc                      <  1:33500/150‑1 (MQ=255)

                                                                                                                              aCCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTcc                     >  2:1034169/1‑150 (MQ=255)

                                                                                                                              aCCTGCTCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTcc                     >  2:1205424/1‑150 (MQ=255)

                                                                                                                                    tCGGCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCGGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCg               <  1:2905348/150‑1 (MQ=255)

                                                                                                                                       gCCAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGt            >  1:6978074/1‑150 (MQ=255)

                                                                                                                                        ccAGCG‑TTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGtt           <  1:6261977/150‑1 (MQ=255)

                                                                                                                                           agcgtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtcc         >  2:7852684/5‑150 (MQ=255)

                                                                                                                                           agcgtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtcc         >  2:4739111/5‑150 (MQ=255)

                                                                                                                                            gcgtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtcct        >  1:1387008/4‑150 (MQ=255)

                                                                                                                                            gcgtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtcct        <  1:64455/147‑1 (MQ=255)

                                                                                                                                             cgtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtcctc       >  1:4485081/3‑150 (MQ=255)

                                                                                                                                              gtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtcctcc      <  2:6335473/149‑1 (MQ=255)

                                                                                                                                              gtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtcctcc      >  2:5197031/2‑150 (MQ=255)

                                                                                                                                               tttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCcg     <  1:5417268/150‑1 (MQ=255)

                                                                                                                                                 tCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCcgcg   >  1:2818199/1‑150 (MQ=255)

                                                                                                                                                 tCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCcgcg   <  1:129369/150‑1 (MQ=255)

ACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTTCCTGG‑AGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCCGCGG  >  NC_008463/3670392‑3670686

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