BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,125	T→G	5.8%	Y416S (TAC→TCC)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,125	0	T	G	5.8%	1053.9 / 11.8	374	Y416S (TAC→TCC)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base T (119/223); new base G (25/6); total (145/229)
Fisher's exact test for biased strand distribution p-value = 8.83e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 379 total aligned reads displayed.

ATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727977‑3728274

 tAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTaa                                                                                                                                                     <  2:106564/150‑1 (MQ=255)

   aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATc                                                                                                                                                   <  1:4342571/150‑1 (MQ=255)

     cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                                                 <  1:5808157/150‑1 (MQ=255)

     cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                                                 <  1:3853543/150‑1 (MQ=255)

     cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                                                 <  1:3741930/150‑1 (MQ=255)

       gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGAGTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGTAATCGAgc                                                                                                                                               >  2:1397063/1‑150 (MQ=255)

        cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgct                                                                                                                                              <  1:2255559/150‑1 (MQ=255)

         aTAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTCGTCGGgctg                                                                                                                                             >  2:69139/1‑150 (MQ=255)

         aTAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTACGGTCAGCGGgctg                                                                                                                                             >  2:6277754/1‑150 (MQ=255)

           aaCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGGGCGGTCATCGGgctgcg                                                                                                                                           >  1:6139745/1‑149 (MQ=255)

               cTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGGCGTCGAgctgcggggc                                                                                                                                       >  2:2755965/1‑145 (MQ=255)

                  cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                                                                                                    <  2:469800/150‑1 (MQ=255)

                   gcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTg                                                                                                                                   <  2:6538883/150‑1 (MQ=255)

                    cATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGa                                                                                                                                  <  1:681866/150‑1 (MQ=255)

                     aTCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTGATCGTGCTGGTGGCCGCTGag                                                                                                                                 >  2:3616065/1‑149 (MQ=255)

                         ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                             <  2:5872106/150‑1 (MQ=255)

                            aGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGcc                                                                                                                          <  1:5579227/150‑1 (MQ=255)

                                 cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCAt                                                                                                                     <  2:5225088/150‑1 (MQ=255)

                                     tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTGCGGTGATCGAGCTGCTGACCGCTGATATCCGCCAGCGTCTgg                                                                                                                 >  2:5548888/1‑149 (MQ=255)

                                        tCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                              <  2:5447643/150‑1 (MQ=255)

                                          tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTTCGGTAGTCGAGCTGCTGACCGCTGATATCCGCCGGCATCTGACGccg                                                                                                            >  2:407464/1‑149 (MQ=255)

                                          tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                            <  1:6048365/150‑1 (MQ=255)

                                          tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                            <  1:5439029/150‑1 (MQ=255)

                                           aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAg                                                                                                           <  2:2649786/150‑1 (MQ=255)

                                           aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGCATCGCGGTGCTGGCCGCTGATCTCCGCCCGCCTCGGGCGccgg                                                                                                           >  2:204198/1‑148 (MQ=255)

                                            aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGAGCTGCTGACCGCGGAGATCCGCCAGCATCTGACGCCGGt                                                                                                          >  2:2437443/1‑150 (MQ=255)

                                             cGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                              >  2:6631979/1‑145 (MQ=255)

                                              ggTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGAGCTGCTGACCGCGGATCTCCGCCAGCAGCTGACGCCAGtgc                                                                                                        >  2:1856750/1‑148 (MQ=255)

                                                 aGTTGCGGCAACAATTCAATGAGCGCTTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGATGATATCCGCCAGCATCTGACGCCAGTACTCa                                                                                                     <  2:6589593/150‑1 (MQ=255)

                                                 aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCa                                                                                                     <  1:1839485/150‑1 (MQ=255)

                                                 aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCa                                                                                                     <  2:1838396/150‑1 (MQ=255)

                                                 aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCa                                                                                                     <  1:777416/150‑1 (MQ=255)

                                                   ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTACTCGGGCTGCGGACCGCGGTTATCCGCCCGCCGCGGACGCCAGTCCGCACt                                                                                                   >  2:4607956/1‑150 (MQ=255)

                                                     gCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAc                                                                                                 <  1:2306553/150‑1 (MQ=255)

                                                     gCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTAAGCGAGCGGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTGCTCACTAc                                                                                                 >  2:3177928/1‑150 (MQ=255)

                                                      cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCCTCGCGCTGCTGACCGCTGAGATCCGCCAGCATCTGACGCCAGTACTCACTCcg                                                                                                >  2:6193640/1‑150 (MQ=255)

                                                       ggCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAcgc                                                                                               <  1:4251124/150‑1 (MQ=255)

                                                       ggCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAcgc                                                                                               <  1:4999776/150‑1 (MQ=255)

                                                          aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTAGGGGTGTCGTGCTGCTGTCCGCTGATGTCGGGCGGCGTCTGTCGGCGGTGCTCTCTAGGggcc                                                                                            >  2:4934336/1‑150 (MQ=255)

                                                          aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGGCGGTCATCGAGCTGCTGACCGCGGGTGTCCGCCAGCAGCTGACGCCAGTACTCGCTACGcgcc                                                                                            >  2:2552482/1‑150 (MQ=255)

                                                            caaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGTGCGGGTGTCGGGCTGCTGACCGCGGGGATCCGCCGGCGTCTGACGCCGGGACTCTCTCCGcgcccc                                                                                          >  1:1252224/1‑150 (MQ=255)

                                                            caaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgcccc                                                                                          <  1:5259908/150‑1 (MQ=255)

                                                              aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGGGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCGCTACGcgccccgc                                                                                        >  2:1748130/1‑150 (MQ=255)

                                                              aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGTATCTGACGCCAGTACTCACTACGcgccccgc                                                                                        <  1:1067484/150‑1 (MQ=255)

                                                              aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgc                                                                                        <  2:5157512/150‑1 (MQ=255)

                                                              aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTGCGGGCGGCGTGCTGCTGGCCGCGGTGGGCCGCCGGCCTCTGGCGCGGGGGCTCGCGCCGcgccccgc                                                                                        >  1:3926565/1‑150 (MQ=255)

                                                               ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTAATCGGGCTGCTGACCGCTGAGATCCGCCCGCGGCGGGCGCCAGTGCTCGCGGCGcgccccgcc                                                                                       >  2:3942680/1‑150 (MQ=255)

                                                                tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgccc                                                                                      <  1:4607414/150‑1 (MQ=255)

                                                                  aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTGCGGTCGTCGAGCTGGTGGCCGCTGATATCCGCCGGCGGCTGGCGGCGGTGCGGACTGCGCGCGCCGccccc                                                                                    >  2:1317390/1‑150 (MQ=255)

                                                                      aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGGAGTCGGGCTGCTGGCCGCGGCTGTCCGCCCGCTTCTGACGCCGGTGCTCGCTGCGCGCCCCGCCCCCAATc                                                                                >  2:1723540/1‑150 (MQ=255)

                                                                      aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTAGGGGGATCGCGCTGCTGGCCGCTGATGTGCGCCCGCGGCTGCCGCCGGTTCTCCCTGCGCGCCCCGCCCCCACTc                                                                                >  2:4229225/1‑150 (MQ=255)

                                                                       gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGGTATCCGCCAGCAGCTGCCGCCAGTCCTCACTACGCGCCCCGCCCCCACTcc                                                                               >  2:3800077/1‑149 (MQ=255)

                                                                        cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                                             <  1:2152682/150‑1 (MQ=255)

                                                                         gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGACGGTAATCGAGCGGCTGACCGCTGAGCTCCGCCCGCATCTGACGCCAGTGCTCACTACGCGCCCCGCCCCCAATCAAc                                                                             >  2:2613534/1‑150 (MQ=255)

                                                                         gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                                             <  2:3350866/150‑1 (MQ=255)

                                                                         gATGTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                                             <  1:2701233/150‑1 (MQ=255)

                                                                                 tGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                                                                                                    >  1:2614897/1‑87 (MQ=255)

                                                                                  ggATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACtt                                                                    <  2:2696448/150‑1 (MQ=255)

                                                                                    aTTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTg                                                                  <  1:1560384/150‑1 (MQ=255)

                                                                                         ccgcGATCTGTGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                             <  2:6088874/150‑1 (MQ=255)

                                                                                         ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGACGGGGGTCGAGCTGCGGGCCGCGGGGGGCCGCCCGCCTCGGGCGCCGGTCCGCACTACGCGCCCCGCCCCCGGGCAGCGTCACACGTTGCgggg                                                             >  2:368634/1‑148 (MQ=255)

                                                                                            cgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTaa                                                          <  1:2936706/150‑1 (MQ=255)

                                                                                             gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAAt                                                         <  1:6606364/150‑1 (MQ=255)

                                                                                             gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAAt                                                         <  2:4503568/150‑1 (MQ=255)

                                                                                              aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATa                                                        <  2:3459626/150‑1 (MQ=255)

                                                                                                cTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAcc                                                      <  1:1694968/150‑1 (MQ=255)

                                                                                                   cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa                                                   <  2:580276/150‑1 (MQ=255)

                                                                                                   cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCa                                                   <  1:2473804/150‑1 (MQ=255)

                                                                                                    gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAg                                                  <  1:859062/150‑1 (MQ=255)

                                                                                                      cAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTGATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAGCAACGTAACACTTTGCGGGGTAATACCGCAGGc                                                >  2:3972616/1‑150 (MQ=255)

                                                                                                        gCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAt                                              <  1:6293044/150‑1 (MQ=255)

                                                                                                          cTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGc                                            <  1:5961339/150‑1 (MQ=255)

                                                                                                           ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCa                                           <  1:769066/150‑1 (MQ=255)

                                                                                                                gcgcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCCTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGTCGCCAGTAGTCACTCCGCGCCCCGCCCCCCATCAACGTCCCCCGTTGCGGTTGAATACCGCTGGCAGGCACGACa                                      >  2:2888770/1‑150 (MQ=255)

                                                                                                                gcgcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACa                                      <  1:4207079/150‑1 (MQ=255)

                                                                                                                 cgcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAACAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACAt                                     >  2:602067/1‑150 (MQ=255)

                                                                                                                   cccAGTGCTGGCCCCACCTCCCCCCCCGTGCGGTCCTCGAGCTGCTGACCGCTGATAGCCGCCCGCATCTGGCGCCCGTGCTCACGACGCGCCCCGCCCCCCCTCAACGGAACACGTTGCGGGTTGATACCGCAGGCATGCACGAGATcc                                   >  1:3836963/1‑150 (MQ=255)

                                                                                                                   cccAGTGCTGGCCCCACATCCCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTACCACTTTGCGGTTTAATACCGCAGGCATGCACGACATcc                                   >  2:759502/1‑150 (MQ=255)

                                                                                                                   cccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATcc                                   <  1:3577790/150‑1 (MQ=255)

                                                                                                                    ccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGCCGCCAGTACTCACTACGCGCCCCGCCCCCGATCAACGGACCACTTTGCGGTTGAATAACGCAGGCATGCACgcccgcca                                  >  1:2446004/1‑143 (MQ=255)

                                                                                                                     cAGTGCTGGCCCCACATCCCCCCCCGTACGGGAATCGAGCTGCGGACCGCGGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACGTACCCCGTTGCGGTGTGATACCGCAGGCATGCACGccatccat                                 >  1:3482264/1‑150 (MQ=255)

                                                                                                                      aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCCGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCCATCCACGTACCACTTTGCGGGTTAATACCGCAGGCATGCACGACATCCATg                                >  2:1538583/1‑150 (MQ=255)

                                                                                                                      aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                                <  2:787007/150‑1 (MQ=255)

                                                                                                                      aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                                <  2:4154172/150‑1 (MQ=255)

                                                                                                                      aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                                <  2:6317348/150‑1 (MQ=255)

                                                                                                                           tGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCGGATATCCGCCAGCATCTGACGCCAGTACGCACTACGCGCCCCGCCCCCAATCAACGGACCACTTTGCGGTGTAATACCGCAGGCATGCACGACATCCCTGCCATc                           >  2:1180568/1‑150 (MQ=255)

                                                                                                                                ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa                      >  2:567796/1‑150 (MQ=255)

                                                                                                                                ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa                      >  1:1023992/1‑150 (MQ=255)

                                                                                                                                ccACATCCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACGTAACACTTTTCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa                      <  2:196653/150‑1 (MQ=255)

                                                                                                                                ccACATCCCCCCCCGAAAGGAAAACGAGCTGCTGACCGCTGATAACAGCCACCATCTGACGCCAGAACTCACAACGCGCCCCGCCCCCAATCAACCAAACACTTTGCGGTATAAAACCGCAGGCATGCACGACAACCATGCCAACAGCCa                      >  1:3277566/1‑150 (MQ=255)

                                                                                                                                  cactCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGTCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgc                    <  2:508742/147‑1 (MQ=255)

                                                                                                                                        ccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAc                                                       <  1:5857818/109‑1 (MQ=255)

                                                                                                                                        ccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATaa              <  1:5040618/150‑1 (MQ=255)

                                                                                                                                          ccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCa                                           >  2:3092930/1‑119 (MQ=255)

                                                                                                                                          ccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCa                                           <  1:3092930/119‑1 (MQ=255)

                                                                                                                                          ccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAcc            <  1:4792613/150‑1 (MQ=255)

                                                                                                                                             cGTACGGTAAGCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAc         <  2:2298899/150‑1 (MQ=255)

                                                                                                                                               tACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCCACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGc       <  2:5721457/150‑1 (MQ=255)

                                                                                                                                               tACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGc       <  1:3137480/150‑1 (MQ=255)

                                                                                                                                               tACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGc       <  1:6576133/150‑1 (MQ=255)

ATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727977‑3728274

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