BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,412,858	G→T	5.6%	A158S (GCT→TCT)	yjfM →	DUF1190 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,412,858	0	G	T	5.6%	1923.2 / 26.4	549	A158S (GCT→TCT)	yjfM	DUF1190 family protein
Reads supporting (aligned to +/- strand): ref base G (191/314); new base T (42/2); total (233/316)
Fisher's exact test for biased strand distribution p-value = 8.57e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 555 total aligned reads displayed.

AGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGT  >  NC_000913/4412709‑4413007

  cagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCt                                                                                                                                                     <  2:5452257/150‑1 (MQ=255)

      aCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcgc                                                                                                                                                 <  1:11889332/150‑1 (MQ=255)

      aCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcgc                                                                                                                                                 <  1:4735445/150‑1 (MQ=255)

        tGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCgcgc                                                                                                                                               <  2:1429239/150‑1 (MQ=255)

         gTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGcc                                                                                                                                              <  1:2021448/150‑1 (MQ=255)

              tCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTc                                                                                                                                         <  1:8824454/150‑1 (MQ=255)

              tCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTc                                                                                                                                         <  1:1726064/150‑1 (MQ=255)

              tCTAAGTACGAAAATTGCGACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTc                                                                                                                                         <  2:3206705/150‑1 (MQ=255)

                    tACGAAAAGTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgcgc                                                                                                                                   <  2:4157953/150‑1 (MQ=255)

                     aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca                                                                                                                                  <  2:7632130/150‑1 (MQ=255)

                     aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTGCGGGATTTTTGTTAATCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTACTCTTTTGCTTCGCGCCCAGTCTGGCgca                                                                                                                                  >  2:10014331/1‑150 (MQ=255)

                      cGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcag                                                                                                                                 <  1:565570/150‑1 (MQ=255)

                              gCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGGTATTAGGAAAGATCGCGATGAGCCGGTTGTTTATAACAGCGGCGGTTCCACTTTTGCTTAGCGCCCAGTCTGGCGCCGCACTTCTg                                                                                                                         >  2:951053/1‑150 (MQ=255)

                                tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                                                                                                       <  2:9876274/150‑1 (MQ=255)

                                tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                                                                                                       <  2:11402455/150‑1 (MQ=255)

                                tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                                                                                                       <  2:5011061/150‑1 (MQ=255)

                                tactacGACAATGTTGAACAGAGGTGGATCCCGGTCGTGTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                                                                                                       <  2:1717154/150‑1 (MQ=255)

                                  ctacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGa                                                                                                                     <  1:11676102/150‑1 (MQ=255)

                                   tacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGAt                                                                                                                    <  1:1327520/150‑1 (MQ=255)

                                      gACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGAATTTTGTTAATCCGGGTTAATCGTAAATATAGCGATGATCCGTTTTTTTATAAAAGCGGCGGTTCCTATTTTTCTTAGCGCCCAGTCTGGCGCAGCACATCTTGTGATTaa                                                                                                                 >  2:5078922/1‑149 (MQ=255)

                                      gACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGGTTTTTTTTAAGCCCGGTTAATACTAAAAATAGCCATGATCCCGTTTTTTATAAAACCCGCCGTTCCTCTTTTTCTTAGCGCCCCCTCAGGCGCACCACATCGGGTGATTaa                                                                                                                 >  2:9357144/1‑149 (MQ=255)

                                         aaTGTTGAACAGAGTTGGATCCCGGTCGTTTGTTGATTTTTTTTAATCCCGGTTATTATTAAAAATAGCGAATAGCCCGTTTTTTTTAAAACCGGCGGGTACTCTTTTTGTTCGCGCCCCCGCGGGCCCACCACATATGGTGGTTaaaaa                                                                                                              >  2:10390301/1‑146 (MQ=255)

                                           tGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTTTTTATAAAAGCGGCGGTTACTCTTTTTCTTAGCGCCCAGTCAGGCGCAGCACGTCTGGTGATTACTccag                                                                                                            >  2:2875093/1‑148 (MQ=255)

                                           tGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTg                                                                                                            <  2:7331055/150‑1 (MQ=255)

                                             ttGAACAGAGTTGGATCCCGGTCGTTTCTGGAATTTTGTTAAGACGGGTTATTCGTAAATATAGCGAATATGCGATTGTTTATAAAAACCGCGGTTACTCAAATTCTTCGCGCCCAGTATGGCGCAGCACATCAGGTGATTAATCATGgc                                                                                                          >  2:4033193/1‑150 (MQ=255)

                                               gAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGgcgc                                                                                                        <  1:2873235/150‑1 (MQ=255)

                                                 aCAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTc                                                                                                      <  1:7041656/150‑1 (MQ=255)

                                                   agagTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCg                                                                                                    <  1:2803647/150‑1 (MQ=255)

                                                        tGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCt                                                                                               <  2:11867596/150‑1 (MQ=255)

                                                        tGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCt                                                                                               <  1:11745667/150‑1 (MQ=255)

                                                         ggATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                                                                                              <  2:4806847/150‑1 (MQ=255)

                                                          gATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTgg                                                                                             <  1:6135675/150‑1 (MQ=255)

                                                           aTCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGc                                                                                            <  1:2164226/150‑1 (MQ=255)

                                                               cGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATAGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTACTCTTTTTCTTCGCGCCCAGTAAGGCGCAGCACATCTGGTGATTACAACAGGCGCTCCGGTTCTGGCaaaa                                                                                        >  2:12163576/1‑150 (MQ=255)

                                                                ggTCGTTTCAGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGGTTGTTTATAACAGCGGCGGGTCCTCTTTTGCGTCCCGCCCACTCTGGCGCACCACATCTGGTGATAACACCAGGCGCACCCGTTCTGGCcaaaa                                                                                       >  2:5628219/1‑150 (MQ=255)

                                                                   cGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga                                                                                    <  1:713429/150‑1 (MQ=255)

                                                                   cGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga                                                                                    <  1:4604513/150‑1 (MQ=255)

                                                                     tttCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCTCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGt                                                                                  >  2:7793363/1‑150 (MQ=255)

                                                                       tCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCt                                                                                <  1:3342999/150‑1 (MQ=255)

                                                                         tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTa                                                                              >  2:8392272/1‑150 (MQ=255)

                                                                         tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTa                                                                              >  1:5048610/1‑150 (MQ=255)

                                                                          ggATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTAc                                                                             <  1:7392894/150‑1 (MQ=255)

                                                                            aTTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActc                                                                           >  2:11574190/1‑150 (MQ=255)

                                                                             tttttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActct                                                                          >  1:10221409/1‑150 (MQ=255)

                                                                             tttttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActct                                                                          >  2:3924018/1‑150 (MQ=255)

                                                                               tttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTc                                                                        <  2:11189256/150‑1 (MQ=255)

                                                                                     aaGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcgg                                                                  <  1:956887/150‑1 (MQ=255)

                                                                                     aaGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcgg                                                                  >  2:8056219/1‑150 (MQ=255)

                                                                                     aaGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcgg                                                                  >  2:6157056/1‑150 (MQ=255)

                                                                                      aGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTTTTTATAACAGCGGCGGGTGCTCGTTTTATTAGCCCCAAGGAGGGCGCAGCACATCTGGTGGGTTATACTGGCGCGCCGGTTCTGGCAAAAAAGATTAATAATCTTAGgggggc                                                                 >  1:5411847/1‑150 (MQ=255)

                                                                                      aGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcggc                                                                 <  1:773423/150‑1 (MQ=255)

                                                                                      aGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcggc                                                                 <  2:450054/150‑1 (MQ=255)

                                                                                      aGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcggc                                                                 <  2:6419429/150‑1 (MQ=255)

                                                                                        ccGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCtt                                                               >  2:5425273/1‑150 (MQ=255)

                                                                                        ccGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCtt                                                               >  2:8874269/1‑150 (MQ=255)

                                                                                        ccGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCtt                                                               >  2:2119810/1‑150 (MQ=255)

                                                                                        ccGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCAGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCtt                                                               >  1:7996216/1‑150 (MQ=255)

                                                                                        ccGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCACTTTTGCTTCGCGCCCCGTCAGGCGCAGCACATCTGGTGATTACTCCTGGCGCTCCCGTTCTGGCCCAAAAGAGTCTTACTCTTCGGGCGGCtt                                                               >  1:2506414/1‑150 (MQ=255)

                                                                                         cGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTATGGTGATTAATCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCttt                                                              >  1:8132352/1‑150 (MQ=255)

                                                                                         cGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCGCTTTTTCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCGCCGGTTCTGGCAAAAAAGAGTCTTAATATTCGGGCGGCttt                                                              >  1:3996590/1‑150 (MQ=255)

                                                                                            gTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTacc                                                           <  2:1896836/150‑1 (MQ=255)

                                                                                             ttATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTacca                                                          <  1:6769536/150‑1 (MQ=255)

                                                                                              tATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTaccac                                                         <  1:8905498/150‑1 (MQ=255)

                                                                                              tATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAAAAGCGGCGGTTCCGCTGTTGCTTCGCGCCCAGTCTGGCGCAGCACTTATGGTGATTACTCCTGGCGCGCCGGGTCTGGCAAAAAAGATTCTTACTCTTCGGGCGGCTTTaacac                                                         >  1:11389085/1‑150 (MQ=255)

                                                                                               aTTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCGTTaccacc                                                        >  2:371218/1‑150 (MQ=255)

                                                                                                   gTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCaaaa                                                    <  1:3427499/150‑1 (MQ=255)

                                                                                                    tAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCaaaaa                                                   >  2:4715548/1‑150 (MQ=255)

                                                                                                    tAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCaaaaa                                                   >  2:6945485/1‑150 (MQ=255)

                                                                                                     aaaGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCaaaaaa                                                  <  2:11729772/150‑1 (MQ=255)

                                                                                                     aaaGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCaaaaaa                                                  <  1:5846854/150‑1 (MQ=255)

                                                                                                         agcgcgATGAGCCGTGTGCTTTTAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTCCTCCTGGCGCTCCGGTTCTGGCAAAAATGAGTCTTACTCTTCGGGCGGCTTTACCACCACAAAAGCGt                                              <  2:6581948/148‑1 (MQ=255)

                                                                                                            gcgATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTa                                           >  2:9164267/1‑150 (MQ=255)

                                                                                                             cgATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTAc                                          <  2:2455874/150‑1 (MQ=255)

                                                                                                                  aGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGttt                                     >  2:8730560/1‑150 (MQ=255)

                                                                                                                    ccGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCCAAAAATAGTCTTACTCTTCGGGCGGCGTTACCACCAAAAAAGCGTCTACCGTTtct                                   >  1:1870907/1‑150 (MQ=255)

                                                                                                                       tttgtttATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTcgc                                >  2:7152435/1‑150 (MQ=255)

                                                                                                                       tttgtttATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTcgc                                >  2:10183634/1‑150 (MQ=255)

                                                                                                                           tttATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCgcggcg                            >  2:8454893/1‑150 (MQ=255)

                                                                                                                           tttATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCgcggcg                            >  2:10892101/1‑150 (MQ=255)

                                                                                                                           tttATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCgcggcg                            >  2:6877538/1‑150 (MQ=255)

                                                                                                                           tttATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCgcggcg                            >  1:10723869/1‑150 (MQ=255)

                                                                                                                           tttATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACGTCTGGGGATTACTCCTGGCGCGCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTGTCTCgcggcg                            >  1:3697687/1‑150 (MQ=255)

                                                                                                                            ttATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGcggcgg                           <  1:8499978/150‑1 (MQ=255)

                                                                                                                             tataACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCggcggc                          >  2:8226510/1‑150 (MQ=255)

                                                                                                                                    gcggcgGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTc                   <  1:987459/150‑1 (MQ=255)

                                                                                                                                         ggTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCGAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGttct              >  1:6386415/1‑150 (MQ=255)

                                                                                                                                           ttCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGttcttc            >  1:8606662/1‑150 (MQ=255)

                                                                                                                                           ttCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGttcttc            >  2:9663773/1‑150 (MQ=255)

                                                                                                                                            tCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTcc           >  2:1164359/1‑150 (MQ=255)

                                                                                                                                            tCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTcc           >  1:10531850/1‑150 (MQ=255)

                                                                                                                                             ccTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCa          <  1:10020633/150‑1 (MQ=255)

                                                                                                                                               tctTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAgc        >  1:1501312/1‑150 (MQ=255)

                                                                                                                                               tctTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAgc        >  2:10597068/1‑150 (MQ=255)

                                                                                                                                                 ttttGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAgcgc      >  2:6006956/1‑150 (MQ=255)

                                                                                                                                                 ttttGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAgcgc      <  1:4390444/150‑1 (MQ=255)

                                                                                                                                                 ttttGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAgcgc      >  2:3401144/1‑150 (MQ=255)

                                                                                                                                                  tttGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGcc     >  1:923409/1‑150 (MQ=255)

                                                                                                                                                   ttGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGGTCTTACAGCGccc    >  2:7921180/1‑150 (MQ=255)

                                                                                                                                                    tGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCg   <  2:11003943/150‑1 (MQ=255)

                                                                                                                                                     gCTTCGCGCCCAGTCTGGCGCAGCACTTCCGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGt  <  2:201990/150‑1 (MQ=255)

AGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTCCAGCGCCCGT  >  NC_000913/4412709‑4413007

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