BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,125	T→G	5.7%	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.7%	1774.6 / 19.4	607	Y416S (TAC→TCC)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base T (175/379); new base G (42/9); total (219/388)
Fisher's exact test for biased strand distribution p-value = 1.39e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 610 total aligned reads displayed.

TAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727978‑3728274

tAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCATTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTaa                                                                                                                                                     <  2:2016575/150‑1 (MQ=255)

 aaGCGGCATAACGCTGTGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAAt                                                                                                                                                    >  1:10167761/1‑150 (MQ=255)

 aaGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCACTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAAt                                                                                                                                                    <  2:1191013/150‑1 (MQ=255)

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

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

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

      gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgc                                                                                                                                               <  1:4726692/150‑1 (MQ=255)

       cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgct                                                                                                                                              >  1:5217849/1‑150 (MQ=255)

         tAACGCTGTGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgctgc                                                                                                                                            <  2:10167761/150‑1 (MQ=255)

            cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGGCCCCGCCCGTGCGGTCGGCGAGctggtgg                                                                                                                                         >  1:4846804/1‑149 (MQ=255)

            cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGa                                                                                                                                         <  1:1123217/150‑1 (MQ=255)

              cTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGAcc                                                                                                                                       <  2:7100305/150‑1 (MQ=255)

               tGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGGGCTGCTGGCCg                                                                                                                                      >  2:10008272/1‑150 (MQ=255)

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

                gcgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGTCCGc                                                                                                                                     >  2:7982476/1‑150 (MQ=255)

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

                 cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTGCGGTCATCGAGCTGCTGACCGcgg                                                                                                                                   >  1:7462123/1‑148 (MQ=255)

                   cATCTGGTAGATGCGACTGTTCTTACGGTAGTTGCGGCAACAATTCAATGAGCGCTTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCTGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGa                                                                                                                                  <  2:9663208/150‑1 (MQ=255)

                   cATCTGGTACATGCGACTGTTCTAACGGGAGTTGCGCCAACAATTCAATCAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTAGCTCCGCCCAGTGCTGGCCCCCCAGCCCCCCCCGTCCGGCCAGCGCACAGCTGACCGcggc                                                                                                                                  >  2:3681501/1‑147 (MQ=255)

                     tCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGata                                                                                                                                >  2:3360211/1‑150 (MQ=255)

                       tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGATATc                                                                                                                              >  2:9116694/1‑150 (MQ=255)

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

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

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

                            gATGCGACTGTTCTAACTGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTTACCGCTGATATCCGCCa                                                                                                                         <  1:2885464/150‑1 (MQ=255)

                            gATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                                                         <  1:3019763/150‑1 (MQ=255)

                                    tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                                                                 <  1:7696434/150‑1 (MQ=255)

                                        ctcaCGCTAGTCGCCGCAACAATTCACTGAGCGATTTCTCGGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGATATCCTCCAGCATCTGACGcc                                                                                                             <  2:9488194/147‑1 (MQ=255)

                                         tAACGGTAGTTGCGACAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                            <  1:3278938/150‑1 (MQ=255)

                                           aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACg                                                                                                               >  1:7137203/1‑145 (MQ=255)

                                           aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTa                                                                                                         <  2:1474426/150‑1 (MQ=255)

                                            cGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTa                                                                                                         <  1:451087/150‑1 (MQ=255)

                                            cGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCACTGATATCCGCCAGCATCTGACGCCAGTa                                                                                                         <  1:5797512/150‑1 (MQ=255)

                                             ggTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTAc                                                                                                        <  2:1537624/150‑1 (MQ=255)

                                              gTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGGACGGGCCGCGGGCTGCTGACCGCTGCTGGCCGCCCGCCGCTGGCGCCCGTACt                                                                                                       >  2:5210653/1‑150 (MQ=255)

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

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

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

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

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

                                                            aaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGTCCCCCCCCGGGCGGGAGTCGGGCTGCTGGCCGGGGATTTCCGCCGGCAGGTGGCGCCAGGACTCACTACGcgcgccg                                                                                         >  1:2796599/1‑150 (MQ=255)

                                                             aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGGTCGGGGTGCGGGCTGGTGACCGCGGGGGGCCGCCAGCAGCTGACGCCGGGACTCACTCCGcgccccgc                                                                                        >  1:4172990/1‑150 (MQ=255)

                                                              ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgcc                                                                                       <  1:2118874/150‑1 (MQ=255)

                                                               tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCAGCCCCCCCCGGACGGGCCTCGCGCGGCTGACCGCTGCTCGCCGCCAGCCGCTGACGCCCGGCCTCACTACGcgccccgccc                                                                                      >  2:1826422/1‑150 (MQ=255)

                                                                ccaTGAGCGAGTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGCGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCgccccgcccc                                                                                     <  2:2508631/148‑1 (MQ=255)

                                                                cAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCgccccgcccc                                                                                     <  2:6172068/150‑1 (MQ=255)

                                                                     aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATc                                                                                <  1:8113728/150‑1 (MQ=255)

                                                                       cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCaa                                                                              <  1:8547232/150‑1 (MQ=255)

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

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

                                                                        gATTTCTCTGGAGTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCATTCAAc                                                                             <  2:8162140/150‑1 (MQ=255)

                                                                          tttCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGTAATCGGGCGGCTGACCGCGGAGATCCGCCAGCAGCTGGCGCCCGGACTCACTACGCGCCCCGCCCCCCATCACCGt                                                                           >  2:2186705/1‑150 (MQ=255)

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

                                                                               ctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGTCCCCCCCCGGACGGGACTCGAGCTGCTGGCCGCTGCTGTCCGCCAGCATCTGACGCCGGTACTCACTACGCGCCCCGCCCCCAATCACCGGAacac                                                                      >  2:5851288/1‑150 (MQ=255)

                                                                                tGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTCAGCGGGCTGCTGTGCGCTGGTGTCCGCCAGCATCTGACGCCAGTCCTCGCTCCGCGCCCCGCCCCCAGGCAACGTAACACt                                                                     >  2:297506/1‑150 (MQ=255)

                                                                                     tCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAGTGGGGCTGCTGACCGCTGATGTCCGCCAGCGTCTGGCGCCCGTGCTCTCTGCGCGCCCCGCCCCCAATCAACGGGACGCGTTGCg                                                                >  2:5366852/1‑150 (MQ=255)

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

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

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

                                                                                        ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGACGGTACTCGGGCTGCTGACCGCGGGTGTCCGCCGGCCTCTGGCGGCAGTGCGCGCGGCGCGCCCCGCCCCCGAGCGGCGTGACACTGTGCGGtt                                                             >  2:3483859/1‑150 (MQ=255)

                                                                                        ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                             >  2:7238114/1‑86 (MQ=255)

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

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

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

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

                                                                                                  cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCACCGTAACACTTTGCGGTTTAATACCGCa                                                   <  2:8771091/150‑1 (MQ=255)

                                                                                                   gcCAGCCTTGCTGCGCCCAGTGCTGGTCCCACATGCCCCCCCGTTCGGTAATCGGGCTGCTGACCGCTGGTAGCCGCCCGCCGGTGACGCCGGTACTCACTACGCGCCCCGCCCCCAATCAGGGGGACCCGTTGCGGGGTAGTACCGCAg                                                  >  2:4985921/1‑150 (MQ=255)

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

                                                                                                    ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTACGGTCAGCGTGCTGCTGGGCGCGGATGTCCGCCCGCGTCTGACGCCCGGGCTCACGACGCGCCCCGCCCCCAATCGACGGAACACTGTGCGGTGTCGTACCGCAgg                                                 >  1:3602572/1‑150 (MQ=255)

                                                                                                    ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGGACGGGGGTCGAGCGGCTGCCCGCGGGTATGCGCCCGCGGCGGCCGCCAGGACTCCCTACGCGCCCCGCCCCCCCTCGCCGGGACGCGGTGCGGTTTAATACCGCCgg                                                 >  2:10258772/1‑150 (MQ=255)

                                                                                                      aGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTGCGGTTATCGGGCTGCTGACCGCTGATATCCGCCCGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGGAACACTTTGCGGGTGAGTACCGCAGGCa                                               >  2:417729/1‑150 (MQ=255)

                                                                                                      aGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCa                                               <  1:4002771/150‑1 (MQ=255)

                                                                                                           tgctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc                                          <  1:3427592/150‑1 (MQ=255)

                                                                                                           tgctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc                                          <  1:10008272/150‑1 (MQ=255)

                                                                                                            gctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACg                                         <  1:1710936/150‑1 (MQ=255)

                                                                                                            gctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACg                                         <  1:9427685/150‑1 (MQ=255)

                                                                                                             ctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGCTAGCCGCCAGCATCTGCCGCCGGTGCGCCCTGCGCGCCCCGCCCCCGAGCCACGGAACACTTTGCGGGGTAATACCGCAGGCATGCACgc                                        >  2:896377/1‑149 (MQ=255)

                                                                                                             ctgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGa                                        <  2:932065/150‑1 (MQ=255)

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

                                                                                                                   ccAGTGCTGGCCCCACATCCCCCCCCGGGCGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACGTTGCGGTTTAATACCGCAGGCATGCACGACAGccc                                  >  2:437587/1‑149 (MQ=255)

                                                                                                                    cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                                 <  2:2776628/150‑1 (MQ=255)

                                                                                                                      gTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGc                               <  2:2679008/150‑1 (MQ=255)

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

                                                                                                                         cTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCGGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCCACGTAACCCTGTGCGGGGTAATACCGCAGGCAGGCACGACAGCCATGcccg                            >  2:9863880/1‑148 (MQ=255)

                                                                                                                            gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAGCAACGGAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCCTCTg                         >  1:8373657/1‑150 (MQ=255)

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

                                                                                                                                  caTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCTgcg                   <  2:3750524/150‑1 (MQ=255)

                                                                                                                                   ctCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcgc                  <  2:5968972/149‑1 (MQ=255)

                                                                                                                                   aTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcgc                  <  1:2422696/150‑1 (MQ=255)

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

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

                                                                                                                                            cGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTACTACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAc         <  2:7521617/150‑1 (MQ=255)

                                                                                                                                             gTACTGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACg        <  2:3539714/150‑1 (MQ=255)

                                                                                                                                               aCGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGcc      >  1:5637060/1‑150 (MQ=255)

                                                                                                                                                cGGTAATCGCGCGGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTTACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCt     <  2:3302271/150‑1 (MQ=255)

                                                                                                                                                 ggTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCtt    >  2:766499/1‑150 (MQ=255)

                                                                                                                                                 ggTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCtt    >  1:5414670/1‑150 (MQ=255)

                                                                                                                                                 ggTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCtt    >  2:5571678/1‑150 (MQ=255)

                                                                                                                                                 ggTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATAAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACCCCtt    >  2:83113/1‑150 (MQ=255)

                                                                                                                                                  gTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTc   >  2:2398865/1‑150 (MQ=255)

TAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727978‑3728274

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