BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,120	T→G	5.6%	T418P (ACG→CCG)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,120	0	T	G	5.6%	753.3 / 3.7	252	T418P (ACG→CCG)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base T (68/163); new base G (19/2); total (87/165)
Fisher's exact test for biased strand distribution p-value = 4.41e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

GGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACG  >  NC_000913/3727972‑3728268

 gCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTAc                                                                                                                                                    >  2:4389689/1‑150 (MQ=255)

   tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCGCCCGTGCgg                                                                                                                                                  >  2:2789235/1‑150 (MQ=255)

   tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCgg                                                                                                                                                  >  2:3758891/1‑150 (MQ=255)

       aaGGGGAATAAGGCTGCGCATCTGGTTGATGCGACTGTTATAACGGTAGTTGGGGCAACAATTCAAGGAGCGTTTTCTATGGATTCGCGGCGTTCTGCGCCAGCCTTGCTGGGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAAt                                                                                                                                              <  1:3730286/150‑1 (MQ=255)

       aaGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAAt                                                                                                                                              <  2:4453708/150‑1 (MQ=255)

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

        aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGATCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATc                                                                                                                                             <  2:3716216/150‑1 (MQ=255)

            gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGTCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgc                                                                                                                                         <  2:1461695/150‑1 (MQ=255)

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

                 aCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGctgctg                                                                                                                                    <  1:1312413/150‑1 (MQ=255)

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

                      gcgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGc                                                                                                                               <  2:3467911/150‑1 (MQ=255)

                      gcgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGc                                                                                                                               <  1:9598/150‑1 (MQ=255)

                      gcgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGc                                                                                                                               <  2:523647/150‑1 (MQ=255)

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

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

                         cATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTGGTCGGGCTGCTGGCCGCTGa                                                                                                                            >  2:3431031/1‑150 (MQ=255)

                             tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGTTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATc                                                                                                                        <  2:91779/150‑1 (MQ=255)

                             tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGGGCCCCCCCGTACGGTCGGCGGGCGGCTGGCCGGTGagagg                                                                                                                        >  1:3756595/1‑146 (MQ=255)

                              ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTACTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                       <  1:2975794/150‑1 (MQ=255)

                              ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGAGTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                       <  1:2669124/150‑1 (MQ=255)

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

                                       gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATc                                                                                                              <  2:4149854/150‑1 (MQ=255)

                                       gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATc                                                                                                              <  2:786279/150‑1 (MQ=255)

                                          tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGACGGTAATCGGGCTGCTGACCGCTGATAGCCGCCCGCagcgga                                                                                                           >  2:3911725/1‑145 (MQ=255)

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

                                            ttCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGc                                                                                                                               >  1:3038631/1‑128 (MQ=255)

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

                                              cTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCGCGCCCGTTGGGGGGGGGGGCTGCTGACCGCGGGTCGCCGCCCGCGTCTGTCGcc                                                                                                       >  1:2874001/1‑150 (MQ=255)

                                               tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCATCCCCCCCCGTGCGGGACGCGAGCTGCTGGCCGCGGATGGCCGCCGGCAGCTGgcgcgc                                                                                                      >  2:960613/1‑148 (MQ=255)

                                                aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCATCCCCCCCCGTACGGTGAGCGTGCTGCTGACCGCTGCTATCCGCCGGCAGCTGACGCCAg                                                                                                     >  2:4238689/1‑150 (MQ=255)

                                                aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTCCGGGCGTCGCGCTGCTGCCCGCTGGTGGCCGCCAGCAGCGGGCGCCAg                                                                                                     >  2:1279439/1‑150 (MQ=255)

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

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

                                                      aGTTGCGGCAACAATTCAATGAGCGAGTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCa                                                                                               <  2:5070428/150‑1 (MQ=255)

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

                                                        ttGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTACGGTCCTCGAGCTGCTGGCCGCTGATATCCGCCGGCATCTGACGCCAGTACTCCCt                                                                                             >  2:3768146/1‑150 (MQ=255)

                                                         tGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTa                                                                                            <  2:2902050/150‑1 (MQ=255)

                                                           cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGGTGTCCGCCCGCATCTGCCGCCCGTGCTCACTAcg                                                                                          >  2:2470344/1‑150 (MQ=255)

                                                           cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTCCGGGCCGCGGGCTGCTGCCCGCGGAGAGCCGCCGGCGGCTGACGCCGGTGCTCCCTCcg                                                                                          >  2:5012661/1‑150 (MQ=255)

                                                             gCAACAATTCAATGAGCGATTTCTCTGGAGTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTGGGCGGGCTGCTGGCCGGTGGTATCCGCCGGCGTCTGGCGCCGGGAGTCAGTGcgcg                                                                                        >  2:4738455/1‑150 (MQ=255)

                                                              caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCACACATCCCCCCCCGTACGGGAATCGAGGTGCTGACCGCTGATATCCGCCAGAGTCTGACGCGAGTGCGCACTgcgcgc                                                                                       >  2:734112/1‑150 (MQ=255)

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

                                                                    ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGGGGGGAAGCGGGCTGGTGGCCGCGGGTCTGCGCCGGCCTCTGAGGCCAGGAGTCCGTGCGcggcccgcc                                                                                 >  1:333203/1‑150 (MQ=255)

                                                                     tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGTAAGCGAGCGGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgccc                                                                                >  2:1897482/1‑150 (MQ=255)

                                                                      cAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCACCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCgccacgcccc                                                                               <  2:3011264/150‑1 (MQ=255)

                                                                       aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTAATCGGGCTGCTGACCGCTGGTGTCCGCCAGCATCTGTCGCCCGTCCGCACTACGCGCCCCGccccc                                                                              >  2:324035/1‑150 (MQ=255)

                                                                       aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCAGCGCGCTGCGGCCCGCGGTGCTCCGCCCGCGTCTGACGCCCGTACTCACTCCGCGCCCCGccccc                                                                              >  1:4210337/1‑150 (MQ=255)

                                                                             cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGAGGTGCTGACCGCTGAGGTCCGCCCGCGTCTGACGCCCGTGCTCGCTGCGCGCCCCGCCCCCaggcaa                                                                        >  2:3386856/1‑145 (MQ=255)

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

                                                                              gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCATCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                                       <  2:166642/150‑1 (MQ=255)

                                                                                      tGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATc                                                                                                                                             <  2:2011644/72‑1 (MQ=255)

                                                                                      tGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATc                                                                                                                                             >  1:2011644/1‑72 (MQ=255)

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

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

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

                                                                                             gccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGt                                                        <  1:3424409/150‑1 (MQ=255)

                                                                                              ccgcGATCTGCGCCAGCCTTGCTGCGCCCCGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                       <  2:4641015/150‑1 (MQ=255)

                                                                                              ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                       <  2:174176/150‑1 (MQ=255)

                                                                                              ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGGCTCGAGCTGCGGGCCGGTGAGATCCGCCAGCCGCTGACGCCCGGCGTCACTCCGCGCGCCGCCCCCAAGCCACGGGCCACGTTGCgggt                                                       >  2:242803/1‑148 (MQ=255)

                                                                                               cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGGGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCCACGTACCACTTTGCGGtgt                                                      >  2:1528693/1‑148 (MQ=255)

                                                                                                    tCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAc                                                 <  2:4816002/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                          ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTCCTCGGGCTGCTGGCCGCTGATCTCCGCCGGCGTCTGACGCCGGTGCTCACTCCGCGCCCCGCCCCCGAGCCACGGACCGCGTTGCGGGGTAAGACCGCAgg                                           >  2:2867611/1‑150 (MQ=255)

                                                                                                           cAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAAGCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCaaag                                          <  1:1798989/150‑4 (MQ=255)

                                                                                                             tccTTGCTGCGACCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCTCTGATATCCTCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACTTAACACTTTTCGGTTTAATACCGCAGGCAt                                        <  1:320173/149‑1 (MQ=255)

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

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

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

                                                                                                                    tgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGAc                                 <  1:3350869/150‑1 (MQ=255)

                                                                                                                         ccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGATCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGGACCCCGTTGCGGTTTAATACCGCCGGCATGCACGACATCCa                            >  2:1839867/1‑150 (MQ=255)

                                                                                                                          cAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGGCCGCTGATCTCCGCCAGCAGCTGACGCCAGTACTCACGACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                           >  1:152907/1‑150 (MQ=255)

                                                                                                                           aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGGCCGCTGATATCCGCCAGCCTCTGACGCCAGTTCTCACGACGCGCCCCGCCCCCAATCAACGTAACACTGTGCGGGTTAATACCGCAGGCATGCACGACAGCCATg                          >  1:1227524/1‑150 (MQ=255)

                                                                                                                             tGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGcc                        <  1:2944645/150‑1 (MQ=255)

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

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

                                                                                                                               cTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCCGGCAGGCACGACCTCCAAGCCAt                      >  2:1117202/1‑150 (MQ=255)

                                                                                                                               cTGGCCCCACATCCCCCCCCGTACGGTAATAGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCCCGACATCCAGGCCAt                      >  1:3839373/1‑150 (MQ=255)

                                                                                                                                tGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGTTATCCGCCCGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATc                     >  1:2216157/1‑150 (MQ=255)

                                                                                                                                 ggCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCt                    <  1:3433754/150‑1 (MQ=255)

                                                                                                                                 ggCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCAACt                    >  2:4691219/1‑150 (MQ=255)

                                                                                                                                  gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg                   <  2:1872478/150‑1 (MQ=255)

                                                                                                                                   ccccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGc                  <  2:1020541/150‑1 (MQ=255)

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

                                                                                                                                      acatcccccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAg               <  1:2064196/145‑1 (MQ=255)

                                                                                                                                       acaTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgc              >  2:4905793/1‑150 (MQ=255)

                                                                                                                                        caTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcg             <  2:4000199/150‑1 (MQ=255)

                                                                                                                                        caTCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcg             <  1:4733473/150‑1 (MQ=255)

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

                                                                                                                                           ccccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCAt          >  2:4055155/1‑150 (MQ=255)

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

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

                                                                                                                                                  cGTACGGTATTCGAGCTGCTGACCTCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAc   <  2:3691019/150‑1 (MQ=255)

                                                                                                                                                  cGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAc   <  2:722420/150‑1 (MQ=255)

                                                                                                                                                  cGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAc   <  1:4977947/150‑1 (MQ=255)

                                                                                                                                                   gTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACg  >  2:1309633/1‑150 (MQ=255)

GGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACG  >  NC_000913/3727972‑3728268

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