BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,125	T→G	7.9%	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	7.9%	513.9 / 2.4	191	Y416S (TAC→TCC)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base T (62/108); new base G (19/2); total (81/110)
Fisher's exact test for biased strand distribution p-value = 1.93e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 193 total aligned reads displayed.

AGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727980‑3728274

  cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGTCCCCCCCCGTGCGGtggtcgg                                                                                                                                                 >  2:1309527/1‑144 (MQ=255)

  cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                                                 >  1:1318903/1‑150 (MQ=255)

   ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAg                                                                                                                                                <  1:939038/150‑1 (MQ=255)

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

     cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGGCGGTAAGCGAgcg                                                                                                                                              >  2:2933389/1‑149 (MQ=255)

         aCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCCTCCCCCCCCGTACGGGACTCGCGctgctg                                                                                                                                          >  2:438527/1‑150 (MQ=255)

         aCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGctgctg                                                                                                                                          >  2:523701/1‑150 (MQ=255)

           gCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGAc                                                                                                                                        <  2:444304/150‑1 (MQ=255)

           gCTGAGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGATGCTGAc                                                                                                                                        <  1:291903/150‑1 (MQ=255)

             tGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCg                                                                                                                                      <  2:2810277/150‑1 (MQ=255)

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

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

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

               cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGGACGGGCTGCGAGctgctgggcgcg                                                                                                                                    >  1:3475215/1‑144 (MQ=255)

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

                   tCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCGCCCGGGCGGTGCTGGAGCGGCGGGCCGcgggtg                                                                                                                                >  1:1354308/1‑145 (MQ=255)

                     tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTAATCGAGCTGCTGACCGCTGATATc                                                                                                                              >  2:1765027/1‑150 (MQ=255)

                        tAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGGCCGGGAAGCGAGGTGCGGGCCGCGGGGATCCGc                                                                                                                           >  2:3045544/1‑150 (MQ=255)

                        tAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGc                                                                                                                           <  2:1056274/150‑1 (MQ=255)

                         aGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGGCCCCCCCCGGACGGTCCTCGGGCTGCTGACCGCTGAGAgccgcc                                                                                                                          >  2:141651/1‑150 (MQ=255)

                             gCGACTGTTCTAACGGTAGGTGCGGCAACACTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCa                                                                                                                      <  2:1879952/150‑1 (MQ=255)

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

                                aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                                                   <  2:3443254/150‑1 (MQ=255)

                                 cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTACGGGGATCGAGCTGCGGACCGCTGATCTCCGCCAGCGTCTg                                                                                                                  >  2:3779092/1‑150 (MQ=255)

                                  tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTTCGGGAATCGGGCGGCTGACCGCTGGGAGCCGCCAGCCTCTgg                                                                                                                 >  2:137764/1‑149 (MQ=255)

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

                                      cTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTGCGCGGGCTGCTGGCCGCGGATCTCCGCCAGCATCTgccgcc                                                                                                             >  2:513784/1‑150 (MQ=255)

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

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

                                           ggTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTGCGGGGCTCGAGCTGCTGTCCGCTGAGCTCCGCCAGCGGCTGGCGccggggc                                                                                                        >  2:1639871/1‑145 (MQ=255)

                                           ggTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTAATCGAGCTGCTGACCGCTGGTATCCGCCCGCGGCTGACGCCAGTAc                                                                                                        >  2:776699/1‑150 (MQ=255)

                                               gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAc                                                                                                    <  2:2443472/150‑1 (MQ=255)

                                               gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAc                                                                                                    <  2:2953995/150‑1 (MQ=255)

                                               gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAc                                                                                                    <  1:1729536/150‑1 (MQ=255)

                                                  gCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCGCCGTGCGGTCATCGAGCGGCTGACCGCTGGTCTCCGCCAGCGTCTGGCGCCAGTCCTcgcgac                                                                                                 >  2:1063288/1‑145 (MQ=255)

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

                                                      caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACgcgc                                                                                             <  1:192616/150‑1 (MQ=255)

                                                        acaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccc                                                                                           <  2:1894362/150‑1 (MQ=255)

                                                           aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTCCGGGGGTCGGGCTGCTGGCCGCTGGTAGGCGCCGGGTTCTGTCGCCGGTGCGCTCGACGcgccccgc                                                                                        >  1:200995/1‑150 (MQ=255)

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

                                                            ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCCGCGAGCTGCGGACCGCGGATAGCCGCCAGCATCTGGCGCCAGGACTCGCTACGcgccccgcc                                                                                       >  2:3347227/1‑150 (MQ=255)

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

                                                                   aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGccccccctc                                                                                >  2:3842558/1‑146 (MQ=255)

                                                                     cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGTACGGTAATCGAGCTGCTGACCGGTGATATCCGCCAGCCTCTGACGCCAGTACTCCCTACGCGCCCCGCCCCCACTccc                                                                              >  2:3150319/1‑148 (MQ=255)

                                                                          tctcTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTaa                                                                         <  2:145195/150‑1 (MQ=255)

                                                                            tctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGGACTCGAGCTGCTGGCCGCTGATATCCGCCGGCATCTGCCGCCCGTACTCCCTGCGCGCCCCGCCCCCGATCAGCGTatcc                                                                       >  2:2624687/1‑147 (MQ=255)

                                                                                gATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACttt                                                                   <  1:637443/150‑1 (MQ=255)

                                                                                 cttCGCCGCGATCTGCGCCAGCCTTGGTGCGCCGTGTGCTGGCCCCACATCCCCCGCCTTACGGTATTCGAGCTTCTGACCGCTGATATCCGCCAGCTTCTGTCGCCAGTACTCACTACTCGCCCCGCCCCCAATCAACGTAACACTTTg                                                                  <  2:3096772/149‑1 (MQ=255)

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

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

                                                                                  ttCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGGCCGCTGGTATCCGCCAGCATCTGACGCCCGTACTCACTACGCGCCCCGCCCCCACTCACCGTACCACGTTGc                                                                 >  2:3930807/1‑150 (MQ=255)

                                                                                   tCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCg                                                                <  2:2788062/150‑1 (MQ=255)

                                                                                    cgccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCgg                                                               <  2:2486160/150‑1 (MQ=255)

                                                                                      ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                             <  2:2409782/150‑1 (MQ=255)

                                                                                       cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                            <  1:3534068/150‑1 (MQ=255)

                                                                                       cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                            <  1:2303967/150‑1 (MQ=255)

                                                                                       cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCGGGTCTCCGCCCGCAGCTGACGCCCGGGCGCCCTACGCGCCCCGCCCCCAATCCACGTAACACTTTGCGGtgg                                                            >  2:1802397/1‑148 (MQ=255)

                                                                                       cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGGCGGTAATCGGGCTGCGGACCGCTGAGATCCGCCCGCATCTGACGCCCGGGCTCACGACGCGCCCCGCCCCCCCTCCACGTCTCGCTGTGCGGgtt                                                            >  2:2392288/1‑150 (MQ=255)

                                                                                       cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTACGGTGCTCGAGCGGCTGACCGCGGATATCCGCCAGCATCTGCCGCCAGTACTCACTCCGCGCCCCGCCCCCCAGCGACGGAGCACGTTGCGGttg                                                            >  2:1399257/1‑149 (MQ=255)

                                                                                        gcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTa                                                           <  1:1309527/150‑1 (MQ=255)

                                                                                        gcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTa                                                           <  1:1309510/150‑1 (MQ=255)

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

                                                                                              tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg                                                     <  2:1175606/150‑1 (MQ=255)

                                                                                              tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCGGTACTCACTACGCGCCCCGCCCCCCATCAACGTAACACTTTGCGGTTTAATCCCg                                                     >  2:795837/1‑150 (MQ=255)

                                                                                                cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCCGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa                                                   >  2:3999750/1‑150 (MQ=255)

                                                                                                 gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTTCGGTAATCGAGCTGCTGCCCGCGGATATGCGCCAGCGTCTGACGCCAGTACTCACTCGGCGCCCCGCCCCCGATCAACGTAACGCGTTGCGGTGGAATACCGCAg                                                  >  2:2101790/1‑150 (MQ=255)

                                                                                                 gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCGGACCGCGGAGGTCCGCCAGCGGCTGGCGCCGGTACGCACTACGCGCCCCGCCCCCGCTCAACGTCACACGGTGCGGGTTCATCCCGCAg                                                  >  1:3135923/1‑150 (MQ=255)

                                                                                                  ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCGCCCCCCGTACGGTATTCGAGCTGCTGACCGGTGTTATCCGCCAGGTTCTGTGGGCAGTAGTGTCTACGGGCGGCGCCCCCAATCAACGTAACACTTTGCGGTTTTATAGCGCTgg                                                 >  2:2696596/1‑150 (MQ=255)

                                                                                                   cAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCaaag                                                <  2:3272408/150‑4 (MQ=255)

                                                                                                      ccTTGCTGCGCCCAGTGCTGGCCCCACATGCCCCCCCGTAGGGTAGTCGAGCTGGTGAGGGCTGATATCCGCCAGGATCTGACGGGAGTAGTCAGTTGGCGCGCCGGGCCCCATGCGCGTTAGAGTTTGGGGTTTATTAGGGCAGGCTTg                                             >  2:3674683/1‑150 (MQ=255)

                                                                                                      ccTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATg                                             <  2:2105211/150‑1 (MQ=255)

                                                                                                        ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTACGCGCGCTGCTGACCGCTGATATCCGCCAGCAGCTGACGCCAGTCCTCACTACGCGCCCCGCCCCCACTCAACGTAACACGTTGCGGTTTACTACCGCAGGCATGcc                                           >  2:3604196/1‑149 (MQ=255)

                                                                                                            tgcGCCCAGTGCTGGCCCCACATCCCCGCCCGTACGGTTATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACTCCAGTACTCACTACGCGCCCCGCCCCCAATAAACTTAACACTTTGCGGTTTAATACCGCAGGCATTCACGAc                                       <  1:3238481/150‑1 (MQ=255)

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

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

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

                                                                                                              cgcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACAt                                     <  1:431422/150‑1 (MQ=255)

                                                                                                                cccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCCTCTGACGCCAGTGCTCACTACGCGCCCCGCCCCCCCTCAACGGAACACTGTGCGGGGGCCTCCCGCAGGCATGCGCGACCTcc                                   >  1:2945638/1‑150 (MQ=255)

                                                                                                                 ccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCCTCTGACGCCAGTACTCTCTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTATTACCGCAGGCATGCACGACATCCa                                  >  1:2796857/1‑150 (MQ=255)

                                                                                                                 ccAGTGCTGGCCCCACATCCCCCCCCGGACGGGAATCGAGCTGCTGACCGCTGATATCCGCCCGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCCACGTCACACTTTGCGGGGTAATACCGCCGGCCTGCCCGACAGCCa                                  >  1:3019835/1‑150 (MQ=255)

                                                                                                                  cAGTGCTGGCCCCACATCCCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGCCATGCACGACATCCAt                                 <  2:2799553/150‑1 (MQ=255)

                                                                                                                      gCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCa                             <  1:3866001/150‑1 (MQ=255)

                                                                                                                        tGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCATCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATc                           <  2:3801201/150‑1 (MQ=255)

                                                                                                                        tGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATc                           <  2:3622385/150‑1 (MQ=255)

                                                                                                                        tGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACGTTGCGGTGTAATACCGCAGGCATGCACGACATCCATGCCATc                           >  1:1090388/1‑150 (MQ=255)

                                                                                                                           ccccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCACCATCTGACGCCAGTCCTCACTACCCGCCCCGCCCCCAATCAACGTAACAATTTGCGGTGTAATAACGCAGGCATGCACGACATCCATGCCAACTGc                        >  1:2442455/1‑150 (MQ=255)

                                                                                                                            cccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGcc                       >  2:219412/1‑150 (MQ=255)

                                                                                                                             ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa                      <  1:278391/150‑1 (MQ=255)

                                                                                                                                    cccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATa               >  1:3711689/1‑150 (MQ=255)

                                                                                                                                     ccccccGTCCGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATaa              <  2:3143927/150‑1 (MQ=255)

                                                                                                                                     ccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATaa              >  1:977214/1‑150 (MQ=255)

                                                                                                                                     acccaagTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGACAGCGCATaa              >  2:2568669/7‑150 (MQ=255)

                                                                                                                                      cccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAc             >  1:3569263/1‑150 (MQ=255)

                                                                                                                                      cccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAc             >  2:1473343/1‑150 (MQ=255)

                                                                                                                                      cccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAc             <  1:18323/150‑1 (MQ=255)

                                                                                                                                       ccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAcc            >  1:3045259/1‑150 (MQ=255)

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

                                                                                                                                                gTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTc   <  1:4064652/150‑1 (MQ=255)

                                                                                                                                                gTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTc   <  1:3860837/150‑1 (MQ=255)

                                                                                                                                                 tAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTcc  <  1:603950/150‑1 (MQ=255)

AGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727980‑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: `‑`