BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,054,888	A→T	9.6%	V186E (GTA→GAA)	ygfB ←	UPF0149 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,054,888	0	A	T	9.6%	1935.2 / 28.7	617	V186E (GTA→GAA)	ygfB	UPF0149 family protein
Reads supporting (aligned to +/- strand): ref base A (160/384); new base T (65/3); total (230/387)
Fisher's exact test for biased strand distribution p-value = 1.19e-27
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 625 total aligned reads displayed.

TACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGA  >  NC_000913/3054739‑3055037

tACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTa                                                                                                                                                       <  2:3310526/150‑1 (MQ=255)

tACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTa                                                                                                                                                       >  1:11046474/1‑150 (MQ=255)

 aCTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAc                                                                                                                                                      <  1:1438312/150‑1 (MQ=255)

      tGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTg                                                                                                                                                 <  1:5154505/150‑1 (MQ=255)

       ggTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTgg                                                                                                                                                <  1:6886245/150‑1 (MQ=255)

         tGCAACAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTGTTTGTTCTTCTGgcg                                                                                                                                              >  1:8379087/1‑150 (MQ=255)

           cagcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGcgcg                                                                                                                                            <  1:6815711/150‑1 (MQ=255)

           cagcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGcgcg                                                                                                                                            <  1:10675789/150‑1 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTCTGGCGcggtc                                                                                                                                         >  2:8191769/1‑150 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggtc                                                                                                                                         <  2:5904592/150‑1 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggtc                                                                                                                                         >  1:9248496/1‑150 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggtc                                                                                                                                         >  2:511550/1‑150 (MQ=255)

               aaaaaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggtcg                                                                                                                                        >  1:4689538/1‑150 (MQ=255)

               aaaaaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTTCTTCTGGggcgggcg                                                                                                                                        >  1:3782260/1‑150 (MQ=255)

                aaaaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggtcgg                                                                                                                                       >  2:7460544/1‑150 (MQ=255)

                   aTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTg                                                                                                                                    <  2:4689319/150‑1 (MQ=255)

                    tCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTCTTATTGCGCGGTGGGGTGc                                                                                                                                   >  1:9315959/1‑150 (MQ=255)

                    tCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGc                                                                                                                                   <  1:3974336/150‑1 (MQ=255)

                      aGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCgg                                                                                                                                 <  2:11125952/150‑1 (MQ=255)

                      aGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGGCGGTTTCgg                                                                                                                                 >  2:2791442/1‑150 (MQ=255)

                        cgcgGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGAt                                                                                                                               <  1:10205019/150‑1 (MQ=255)

                          cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                             <  2:3048616/150‑1 (MQ=255)

                          cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                             <  2:1729059/150‑1 (MQ=255)

                          cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                             <  1:9494187/150‑1 (MQ=255)

                          cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                             <  1:10914510/150‑1 (MQ=255)

                          cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                             <  1:7666233/150‑1 (MQ=255)

                               cTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTaaa                                                                                                                        <  1:2376111/150‑1 (MQ=255)

                                tGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAg                                                                                                                       <  2:5402633/150‑1 (MQ=255)

                                tGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAg                                                                                                                       <  2:6635776/150‑1 (MQ=255)

                                tGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAg                                                                                                                       <  2:9476089/150‑1 (MQ=255)

                                tGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAg                                                                                                                       <  1:10361558/150‑1 (MQ=255)

                                        tGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTg                                                                                                               <  1:1578297/150‑1 (MQ=255)

                                         gCATTTGCTCGTCCAGGGCATGTCGGCGACGCTGAATCTCTTGGCGTGATATCTCACTCATAACACTCTCCTTACTTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTgg                                                                                                              <  1:4855/150‑1 (MQ=255)

                                         gCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTgg                                                                                                              <  2:5599827/150‑1 (MQ=255)

                                         gCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTTCTTCTGGCGCGGTCGGGTGCGGGTGAGTAAAGGTGTCGTgg                                                                                                              >  2:8739968/1‑150 (MQ=255)

                                          cATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGc                                                                                                             <  2:4571031/150‑1 (MQ=255)

                                          cATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGc                                                                                                             <  2:9906236/150‑1 (MQ=255)

                                          cATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTATGGGGCGGTCGGGTGCGGGTGAGTAAAGGTGTCGTGGc                                                                                                             >  1:7334403/1‑150 (MQ=255)

                                           aTTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCa                                                                                                            <  2:3860172/150‑1 (MQ=255)

                                           aTTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCa                                                                                                            <  1:4491640/150‑1 (MQ=255)

                                           aTTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGGTGCGGATGCGTTAAGGTGTCGTGGCa                                                                                                            >  2:1582335/1‑150 (MQ=255)

                                            tttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTATTGTGGCGCGGTGGGGTGGGGTTGAGGGAAGGTGTGGTGGGAt                                                                                                           >  2:349785/1‑150 (MQ=255)

                                            tttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTATTGCGCGGTCGGGTGCGGGTGAGTAAAGGTGTCGTGGCAt                                                                                                           >  1:10019937/1‑150 (MQ=255)

                                            tttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAAATCTGGCGCGGGCGGGTGCGGATGAGTAAAGGTGTCGTGGCAt                                                                                                           >  1:3710394/1‑150 (MQ=255)

                                               gCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAAc                                                                                                        <  1:6174525/150‑1 (MQ=255)

                                                cTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGGCGGGTGCGGGTGAGGAAAGGTGTAGTGGCATaaaa                                                                                                       >  1:8200486/1‑148 (MQ=255)

                                                   caccaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTGACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTATTCTGGCGCGGTCGGTTGCGGGTGAGGAAAGGTGTCGTGGGATAAAAgcg                                                                                                    >  1:7975762/1‑150 (MQ=255)

                                                      caGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCgg                                                                                                 <  2:10400471/150‑1 (MQ=255)

                                                       aGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGc                                                                                                <  2:8075528/150‑1 (MQ=255)

                                                         ggCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCaa                                                                                              <  2:6353526/150‑1 (MQ=255)

                                                         ggCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCaa                                                                                              <  1:5390176/150‑1 (MQ=255)

                                                             tGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTATGGCGCGGGCGGTTGCGGATGGGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACg                                                                                          >  1:1486950/1‑150 (MQ=255)

                                                                cGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAAc                                                                                       <  2:10521517/150‑1 (MQ=255)

                                                                   cGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTAATATGGCGCGGTGGGGTTCGGATTAGTAAAGGTGTCGTTGCAAAAAAAAGCCGCAACACCaaacta                                                                                    >  2:1009313/1‑146 (MQ=255)

                                                                      cGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTc                                                                                 <  2:8696803/150‑1 (MQ=255)

                                                                      cGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTc                                                                                 <  1:5740503/150‑1 (MQ=255)

                                                                        cTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACATATTCga                                                                               <  2:4883199/150‑1 (MQ=255)

                                                                         tGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGGTGCGGGTGAGTAAAGGTGTCGTTGCATAAAAGCGCGGCAACACCAACGGATTCgat                                                                              >  2:7865888/1‑150 (MQ=255)

                                                                          gAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACGTCTGGGGCGGGCGGTTGCGGGTGGGGAAAGGTGTCGTGGCAAAACAGCGCGGCAACACCAACGTATTCgatg                                                                             >  2:4414120/1‑150 (MQ=255)

                                                                           aaaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACATCAGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgatga                                                                            >  2:931943/1‑150 (MQ=255)

                                                                           aaaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgatga                                                                            >  2:7272837/1‑150 (MQ=255)

                                                                            aaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCAGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgatgat                                                                           <  1:5497762/150‑1 (MQ=255)

                                                                             aCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTATTGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAAAAGCGCGGCAACACGAACGGATTCgatgata                                                                          >  2:1842682/1‑149 (MQ=255)

                                                                              ctctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGAtct                                                                         <  1:8878063/150‑1 (MQ=255)

                                                                               tctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTATGGCGCCGTGGGGTGCGGAGGAGTAAAAGTGTAGTGGCATAAAAGCCCCGCAACAACAACGGATTAGAAGTtctc                                                                        >  2:8470290/1‑150 (MQ=255)

                                                                               tctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGG‑GTTTGTTCTTATTGCGCGGTCGGGTGCGGATGAGGGAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTAGATTAtctat                                                                       >  1:3552243/1‑148 (MQ=255)

                                                                                 ttGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCtt                                                                      <  2:6508939/150‑1 (MQ=255)

                                                                                    ccGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTTCTTGAGGCGCGGGGGGGTGCGGGTGAGTAAAGGTGTCGTGGCAATAAAGCGCGGCAAAAACAACGGATTAGGTGAtatcataaa                                                                   >  1:2728563/1‑145 (MQ=255)

                                                                                      gggATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGc                                                                 <  1:5041827/150‑1 (MQ=255)

                                                                                            tctcACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTCTGGGGCGGTCGGTTGCGGGTGTGTAGAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTTTTCGGTGATGTCTTCATGCGGCAtt                                                           >  2:7166965/1‑150 (MQ=255)

                                                                                             ctcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACATCTGGCGCGGGCGGGTTCGGATGAGGAAAGGTGTCGTTGCATAACAGAGCGGCAACACGAAAGAATTCGAAGATATCTTCAAGCGACAttt                                                          >  2:3996624/1‑150 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTATTATTGCGCGGTCGGTTTCGGATGAGTAAAGGTGTTGTGGCATAAAAGCGCGGCAACAAGAACGTATTAGATGATATATTCAAGCGAAAttta                                                         >  2:2517437/1‑149 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTc                                                         <  1:1996631/150‑1 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGGTGCGGATGTGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACCAACGTATTCGATGATCTCTTCAAGCGACATTTc                                                         >  2:8035549/1‑150 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGGCGGTTGCGGCTGAGTAAAGGGGGTGTGGCATAACAGCGCCGCCACACGCACGTATTCGATGTTCTCTTCAAGCGACATTTc                                                         >  2:9122670/1‑150 (MQ=255)

                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCaa                                                       <  1:10517576/150‑1 (MQ=255)

                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCaa                                                       <  1:222645/150‑1 (MQ=255)

                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTAATTGTTGGGCGGTCGGTTTCGGTTGTTTTAAGGTGTGGTGGCAAAACAGCCCCGCAACAACAAACTAATCGATTATATATTAAAGAGACAATTaaa                                                       >  1:5912243/1‑150 (MQ=255)

                                                                                                 ctcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAg                                                      <  2:1345922/150‑1 (MQ=255)

                                                                                                    aTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTAATATTGCGCGGTCGGTTTCGGGTGTGTTAAGGTTTGGTGGGATAAACGCGCGGCAACACGAACGGATTAGATGATGTCATCAAGCGACATATAAAGcac                                                   >  2:10316331/1‑148 (MQ=255)

                                                                                                      aaCACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCtt                                                 <  1:9958997/150‑1 (MQ=255)

                                                                                                       acacTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTc                                                <  1:4178820/150‑1 (MQ=255)

                                                                                                                ttACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCCAAAAAGCGCGGCAACACGAACGGATTCGCTGGTCTCTTCAAGCGCCATTTCAAGCTCTTCCGGAtcttc                                       >  2:678566/1‑150 (MQ=255)

                                                                                                                      tttttGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttcgtcttc                                 <  2:10695339/150‑1 (MQ=255)

                                                                                                                          tGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTtcgtcg                             <  1:8665657/150‑1 (MQ=255)

                                                                                                                            tttttAGTGTAGAGTCGGTTTTTTTAATTCTGGCGCGGTCGGTTTCGGATTAGTAAAGGTTTCGTGGCATAAAAGCCCGGCAACACCAAACTATTAGATGATCTAATCAAGCGACAATTCAAGCTCATACAGATAATAGTCATAGTCGTa                           >  2:2134240/1‑150 (MQ=255)

                                                                                                                             ttttAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTaa                          <  1:7460544/150‑1 (MQ=255)

                                                                                                                              tttAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAAc                         <  1:3087966/150‑1 (MQ=255)

                                                                                                                                  gtgtAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAg                     <  1:3246070/150‑1 (MQ=255)

                                                                                                                                    gtAGAGTCGGTTTTTTTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGtt                   <  1:8521668/150‑1 (MQ=255)

                                                                                                                                    gtAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGtt                   <  2:6120255/150‑1 (MQ=255)

                                                                                                                                    gtAGAGTCGGTTTTTGTACGTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACCGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGtt                   <  2:7006613/150‑1 (MQ=255)

                                                                                                                                      agagTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTgc                 <  1:9588789/150‑1 (MQ=255)

                                                                                                                                          tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCaa             <  2:2070197/150‑1 (MQ=255)

                                                                                                                                           cGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAAt            >  1:5998555/1‑150 (MQ=255)

                                                                                                                                           cGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAAt            <  1:459454/150‑1 (MQ=255)

                                                                                                                                            ggTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATg           >  2:9615461/1‑150 (MQ=255)

                                                                                                                                             gTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGt          >  2:2225325/1‑150 (MQ=255)

                                                                                                                                               ttttGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTa        >  1:1686084/1‑150 (MQ=255)

                                                                                                                                                 ttGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACg      <  1:9085081/150‑1 (MQ=255)

TACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGA  >  NC_000913/3054739‑3055037

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: `‑`