BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	BS168_reference_genome	2,581,732	(T)_6→7	intergenic (‑115/+39)	pstS ← / ← pbpA	phosphate ABC transporter (phosphate binding lipoprotein)/transpeptidase (penicillin‑binding protein 2A)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	2,581,726	1	.	T	93.7%	724.9 / 39.3	207	intergenic (‑109/+45)	pstS/pbpA	phosphate ABC transporter (phosphate binding lipoprotein)/transpeptidase (penicillin‑binding protein 2A)
Reads supporting (aligned to +/- strand): ref base . (6/7); new base T (75/118); total (82/125)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

GCAGCCATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑G‑‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAA  >  BS168_reference_genome/2581577‑2581870

gcagcCATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAGTTTGTCCTACACGTATAGAGTGACGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑g                                                                                                                                                       <  1:4121734/150‑1 (MQ=255)

   gcCATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑G‑‑t‑tt                                                                                                                                                 >  1:1346084/1‑150 (MQ=255)

     cATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑G‑‑t‑ttt‑t                                                                                                                                              <  1:3884110/150‑1 (MQ=255)

               atagcccccAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑aaaa                                                                                                                                    <  1:3466488/143‑1 (MQ=255)

               aTAAGCACCAATTTGTTTGTTTTCATTTCCTGAATTCCCCCTGCGTTTAATGTCGAATTTGACCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGATTTGAACCGAATGATAAGGTTTTGTAAAACAGTGCGG‑‑T‑TTT‑TTGTAT‑aaaa                                                                                                                                    >  1:2222446/1‑150 (MQ=255)

                tAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑aaaaa                                                                                                                                   <  1:2924203/150‑1 (MQ=255)

                tAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑aaaaa                                                                                                                                   >  1:700930/1‑150 (MQ=255)

                   gCACCATTTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑aaaaaaaa                                                                                                                                <  1:250519/150‑1 (MQ=255)

                     aCCAATTTGTTTTTTTTCATTTCCTGAATCCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAAct                                                                                                                              <  1:3156934/150‑2 (MQ=255)

                      ccAATTTGTTTTTTCTCAGTTCCTGAATCCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGc                                                                                                                             <  2:2945018/150‑1 (MQ=255)

                        aaTTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTAAAATGTGTAATTTGTCCTTCACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAATTCGccc                                                                                                                           <  2:799008/150‑1 (MQ=255)

                         aTTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCca                                                                                                                          >  1:2712767/1‑150 (MQ=255)

                            tGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCAcat                                                                                                                       <  1:3388578/150‑1 (MQ=255)

                             gTTTTTTTTCATTTCCTGAATTCCCCTTGCGTATATTGTGTAATTTGTCCTACACGTATAGAGTAAAGGTCCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAACCAGTGC‑GT‑T‑TTT‑TTGTAT‑‑AAAAAAACGCCCAcatca                                                                                                                     <  2:2184309/150‑1 (MQ=255)

                             gTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCAcatc                                                                                                                      >  1:14280/1‑150 (MQ=255)

                                 tttttCATTTCCTGAATTCCCCTTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCAtgt                                                                                                                  <  2:3789290/150‑1 (MQ=255)

                                 tttttCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCAtgg                                                                                                                  >  1:311191/1‑149 (MQ=255)

                                    ttCATTTCCTGAACTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAGGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAg                                                                                                               <  1:857509/150‑1 (MQ=255)

                                         ttCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGGGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAGACGCCCACATCATGTGAGCGttt                                                                                                          <  1:1851336/150‑1 (MQ=255)

                                         ttCCTGAATTCCCCCTGCGTATAATGTGTAATTTGCCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAACCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGttt                                                                                                          <  1:1134836/150‑1 (MQ=255)

                                          tCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTCCACGTATGGAGTGAAGGTCCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGtttt                                                                                                         <  2:3007620/150‑1 (MQ=255)

                                          tCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGtttt                                                                                                         <  2:3742700/150‑1 (MQ=255)

                                          tCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACAGCGTGTGcgcggtct                                                                                                         >  1:3464174/1‑146 (MQ=255)

                                           ccTGAATTCCCCTTGCGTTTAATGTGTAATTTGTCCTCCACGTATAGGGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGttttt                                                                                                        <  1:2849392/150‑1 (MQ=255)

                                           ccTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGttttt                                                                                                        <  2:690592/150‑1 (MQ=255)

                                            cTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCGG‑‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACCTCATGGGAGCGgttttt                                                                                                       >  2:3348004/1‑150 (MQ=255)

                                             tGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGttttttt                                                                                                      <  1:4505673/150‑1 (MQ=255)

                                                  tCCCCCTGCGTATAATGTGTATTTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAATACGCCCACATCATGTGAGCGTTTTTTTTTAGt                                                                                                 <  1:4782009/150‑1 (MQ=255)

                                                   cccccTGCGTATAATGTGTAATTTGTCCTCCACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGtt                                                                                                <  2:2372232/150‑1 (MQ=255)

                                                     cccTGCGTATATTGTGTAATTTGTCCTCCACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTAt                                                                                              <  1:3772543/150‑1 (MQ=255)

                                                      ccTGCGTCTAATGTGTCATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAGGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAACCGCCCACATCATGTGAGCGTTTTTTTTTAGTTATc                                                                                             <  2:4569934/150‑1 (MQ=255)

                                                       cTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTGTTTTTTTAGTTATCa                                                                                            >  2:3461671/1‑150 (MQ=255)

                                                         gCGTATAATGTGTAATTTGCCCTACATGTATAGAGTAAAGGTTCTATGTAAAAGGGGGTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACTCCCACATCATGTGAGCGTTTTTTTTTAGTTATCaga                                                                                          <  1:3275790/150‑1 (MQ=255)

                                                          cGTATAATGTGAAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGACCCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAATACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCagaa                                                                                         <  2:2456301/150‑1 (MQ=255)

                                                           gTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTATTTTTTTAGTTATCagaag                                                                                        >  1:1479556/1‑150 (MQ=255)

                                                             ataATGTGTAATTTGTCCTAAACGTATAGAGTAAAGGTTCTATGTAAAAAAGGTTTGAACCGAATGGTAAGGTTTTGTAAATCAGTGC‑G‑‑T‑TTTGTTGTAT‑AAAAAAAACGCCCACATCATGAGAGCGAATTGTTTTAGTTATCAGAAGAc                                                                                      >  1:1089467/1‑150 (MQ=255)

                                                              taATGTGTATTTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAGGGTTTTGTAAACCAGTGC‑GT‑T‑TTT‑TTGCAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACg                                                                                     <  1:2053747/150‑1 (MQ=255)

                                                               aaTGTGTAGTTTTCCCTACACGTATAGTGAAAGGGTTTTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTTTAAACCAGTGC‑GT‑T‑TTT‑TTTTAT‑AAAAAATACCCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGt                                                                                    <  1:445174/150‑1 (MQ=255)

                                                                 tgtgTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCACGTGAGCGCAGTTCTTTAGTTATCAGAGGACGTtg                                                                                  >  2:2759118/1‑150 (MQ=255)

                                                                   tgtAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTGTTTTTATAGTTATCAGAAGACGTtgtg                                                                                >  2:4555388/1‑150 (MQ=255)

                                                                    gttaTTTGTCCAGCTGGTAAAGGGTTAGGGTTCATTGTACAATGTGTTTGAACCGAATGTTTAGTTTTTTTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAATTCGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTgtgt                                                                               <  2:284005/147‑1 (MQ=255)

                                                                    gtACTTTGACCAACCCGTATAGAGTAATGGTTCTATGTAAAATGGGTTTGACCTGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTgtgt                                                                               <  2:1660256/150‑1 (MQ=255)

                                                                     tattttGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAGAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACCCGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGtt                                                                              <  1:1605835/147‑1 (MQ=255)

                                                                     tAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGtt                                                                              >  1:3545759/1‑150 (MQ=255)

                                                                        tttGTGCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGTACCGATTGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACTCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTc                                                                           <  1:876027/150‑1 (MQ=255)

                                                                          tGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTg                                                                         <  1:732931/150‑1 (MQ=255)

                                                                          ggTCCTACCCGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCTT‑‑T‑TTT‑TTGTAT‑AAAAAACACACCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTg                                                                         <  1:3733634/149‑1 (MQ=255)

                                                                           gTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTATTTTTTTAGTTATCAAAAGACGTTGGGGTTTTTgg                                                                        >  2:639855/1‑149 (MQ=255)

                                                                             ccTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTATTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCAt                                                                      <  2:1520015/150‑1 (MQ=255)

                                                                              cTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCAtt                                                                     <  2:4650623/150‑1 (MQ=255)

                                                                               tACACGTATAGAGAAAAGGTTGCATGTAAAAGGGGTTTGAAGCGAAGGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCAttt                                                                    <  1:4804146/150‑1 (MQ=255)

                                                                                acacGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAACCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCAtttt                                                                   <  2:2822264/150‑1 (MQ=255)

                                                                                 cacGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTttc                                                                  <  2:1346084/150‑1 (MQ=255)

                                                                                 cacGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTttc                                                                  <  1:660787/150‑1 (MQ=255)

                                                                                  acGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTttct                                                                 <  2:4586932/150‑1 (MQ=255)

                                                                                   cGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATGAGTGC‑GT‑T‑TTT‑TTGTATAAAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTttct                                                                 <  1:449438/150‑1 (MQ=255)

                                                                                     tataGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCa                                                              >  2:774434/1‑150 (MQ=255)

                                                                                       taGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGAGTCTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAt                                                            >  1:1193434/1‑150 (MQ=255)

                                                                                        agagTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGACCGTCTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAtt                                                           >  1:4425172/1‑150 (MQ=255)

                                                                                         gagTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAttt                                                          <  2:2008722/150‑1 (MQ=255)

                                                                                          agaaaaGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAtttt                                                         <  2:3713838/147‑1 (MQ=255)

                                                                                          agTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGT‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAAtttt                                                         <  1:1521767/150‑1 (MQ=255)

                                                                                            tAAGGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGCTT‑‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGt                                                       <  1:3991352/150‑1 (MQ=255)

                                                                                            tAAAGGTCCTATGTAAAAAGGGTTTGAACCGAATGTAAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAA‑AAAAAAGACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGt                                                       <  2:1249552/150‑1 (MQ=255)

                                                                                              aaGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGAGTATTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATCTTCTTCAATTTTGTCt                                                     >  1:1958437/1‑150 (MQ=255)

                                                                                              aaGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAACCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCt                                                     <  2:86845/150‑1 (MQ=255)

                                                                                                    ttATGTAAAAAGGGTTTGAACGGATTGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTc                                               <  2:1076700/149‑1 (MQ=255)

                                                                                                    cTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTc                                               <  2:3943588/150‑1 (MQ=255)

                                                                                                     tATGTAAAAGGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCTTTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACCCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCt                                              <  1:567747/150‑1 (MQ=255)

                                                                                                      ttGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGGGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCtt                                             <  2:2895197/149‑1 (MQ=255)

                                                                                                      aTGTCAAAGGGGTTTGACCCGATTGTTAAGGTTTTGTAAACCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAACAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCtt                                             <  2:3079735/150‑1 (MQ=255)

                                                                                                       tGTAAAAGGGGTTTGAACCGAATGTTAAGGTTTTGTAAACCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCttt                                            <  2:953957/150‑1 (MQ=255)

                                                                                                          aaaaaGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAATCCCCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGc                                         <  1:2429070/150‑1 (MQ=255)

                                                                                                            aaaGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCGTTTTCTTCAATTTTGTCTTTATTCTTTTGCTg                                       >  2:4036228/1‑150 (MQ=255)

                                                                                                              aGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAGAAAAACGCCCACATCATGTGAGCGTTATTTTTTAGTTATCAGAAGACGTGGGGTTTTCTGCATTTTCTTCAATTTTTACTTTATTCTTTTCCTGAg                                     >  2:4233943/1‑150 (MQ=255)

                                                                                                              aGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTCTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTgtg                                     >  2:1282089/1‑148 (MQ=255)

                                                                                                                ggTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAATACGCCCACATCATGTGAGCGTTATTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTa                                   <  1:46401/150‑1 (MQ=255)

                                                                                                                   ttGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GTTT‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGGTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGTTGAGTAcc                                 >  2:4393430/1‑150 (MQ=255)

                                                                                                                     gAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTATTTTTTTAGTTATCAGAAGACGTTGTGTTTTTAGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGTGTACcttc                              >  2:1796793/1‑150 (MQ=255)

                                                                                                                       aCCGAATGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGGAGGCGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGGGTACcttctt                            >  2:351496/1‑150 (MQ=255)

                                                                                                                          gAATGTTAAGGTTTTGTAAATCAGTGCTT‑‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGt                         <  1:809464/150‑1 (MQ=255)

                                                                                                                           aaTGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTc                        <  2:1129588/150‑1 (MQ=255)

                                                                                                                           aaTGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAACCGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTc                        <  2:14280/150‑1 (MQ=255)

                                                                                                                            ttGTTAGGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACCCCCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTcc                       <  1:1629512/149‑1 (MQ=255)

                                                                                                                             tGTTAAGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAACACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCt                      <  2:2830944/150‑1 (MQ=255)

                                                                                                                               ttAGGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAATACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGc                    <  1:4709664/150‑1 (MQ=255)

                                                                                                                                  aGGTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTc                 <  1:1845304/150‑1 (MQ=255)

                                                                                                                                  aGGTTTTGTAAATCAGTGC‑G‑‑TATTT‑TTGTAT‑AAACACCACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGTAGACGTTGTGTTTTCTGCATTTTCTTCAATCTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTc                 <  2:535323/150‑1 (MQ=255)

                                                                                                                                   ggTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAATACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCg                <  1:94334/150‑1 (MQ=255)

                                                                                                                                   ggTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCg                >  2:3750341/1‑150 (MQ=255)

                                                                                                                                   ggTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGGTTTCTGCATTTTCTTCAATTTTGGCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCg                >  2:904304/1‑150 (MQ=255)

                                                                                                                                    gTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCTCATCATGTGAGCGTATTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGa               >  1:1833293/1‑150 (MQ=255)

                                                                                                                                    gTTTTGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGa               >  2:3547020/1‑150 (MQ=255)

                                                                                                                                     ttttGTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGAt              >  2:1871248/1‑150 (MQ=255)

                                                                                                                                         gTAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGGCGTTGTGGTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCa          >  2:2092024/1‑150 (MQ=255)

                                                                                                                                          tAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAACCGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAg         <  2:1556397/150‑1 (MQ=255)

                                                                                                                                          tAAATCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTTCCTTCTTTGTCCTGCTTCGATTTCAg         >  1:3320787/1‑150 (MQ=255)

                                                                                                                                          tAAACCAGTGCTT‑‑T‑TTT‑TTGTAT‑AAAAAACACGCCCCCATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAg         <  2:3221567/150‑1 (MQ=255)

                                                                                                                                           aaaTCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAATACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGc        <  1:2676060/150‑1 (MQ=255)

                                                                                                                                           aaaTCAGTGC‑GT‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTGTTTTTTTAGTTATCAGAAGGCGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGGCCTGCTTCGATTTCAGc        >  2:1186935/1‑150 (MQ=255)

GCAGCCATCAGAAGCATAAGCACCAATTTGTTTTTTTTCATTTCCTGAATTCCCCCTGCGTATAATGTGTAATTTGTCCTACACGTATAGAGTAAAGGTTCTATGTAAAAAGGGTTTGAACCGAATGTTAAGGTTTTGTAAATCAGTGC‑G‑‑T‑TTT‑TTGTAT‑AAAAAAAACGCCCACATCATGTGAGCGTTTTTTTTTAGTTATCAGAAGACGTTGTGTTTTCTGCATTTTCTTCAATTTTGTCTTTATTCTTTTGCTGAGTACCTTCTTTGTCCTGCTTCGATTTCAGCTCAAAA  >  BS168_reference_genome/2581577‑2581870

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