BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,126	A→G	6.0%	Y416H (TAC→CAC) ‡	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,126	0	A	G	6.0%	596.5 / 6.9	201	Y416H (TAC→CAC) ‡	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base A (56/121); new base G (10/2); total (74/127)
Fisher's exact test for biased strand distribution p-value = 5.66e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 206 total aligned reads displayed.

AAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCA  >  NC_000913/3727979‑3728275

aaGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAAt                                                                                                                                                     <  1:2506063/150‑1 (MQ=255)

 aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATc                                                                                                                                                    <  1:882361/150‑1 (MQ=255)

 aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATc                                                                                                                                                    <  1:1917365/150‑1 (MQ=255)

 aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATc                                                                                                                                                    <  2:3272215/150‑1 (MQ=255)

 aGCGGCATAAAGCTGCGCTTCTGGTAGATGCTACTGTTGTAACGGTAGTTGCGGAAACAATTCAATGAGCGATTTCTCTGTATTCGCCGCTATCTGCTCATGCCTTGCTGCGCCCCGTGGTGGCCCCAAAT‑CCCCCCCCGTACGGTAATc                                                                                                                                                    <  2:185658/150‑1 (MQ=255)

  gCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCg                                                                                                                                                   <  2:1958693/150‑1 (MQ=255)

  gCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCg                                                                                                                                                   >  2:1840836/1‑150 (MQ=255)

  gCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCg                                                                                                                                                   <  2:1636514/150‑1 (MQ=255)

   cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGa                                                                                                                                                  <  1:383184/150‑1 (MQ=255)

    ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCNAg                                                                                                                                                 <  1:2593059/150‑1 (MQ=255)

    ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAg                                                                                                                                                 <  2:3032038/150‑1 (MQ=255)

    ggCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAg                                                                                                                                                 <  2:2645771/150‑1 (MQ=255)

      cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAgct                                                                                                                                               <  1:940919/150‑1 (MQ=255)

       aTAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCT‑CCCCCCCCGTTCGGGCATCGAgctg                                                                                                                                              >  2:3218989/1‑150 (MQ=255)

       aTAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCT‑CCCCCCCCGTACGGGACTCGAgctg                                                                                                                                              >  2:1545192/1‑150 (MQ=255)

       aTAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAGTCGAgctg                                                                                                                                              >  2:279641/1‑150 (MQ=255)

        tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAgctgc                                                                                                                                             <  1:24131/150‑1 (MQ=255)

        tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAgctgc                                                                                                                                             <  1:636875/150‑1 (MQ=255)

          aCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAG‑CCCCCCCCGTCCGGGACTCGAGctgctg                                                                                                                                           >  2:1115089/1‑150 (MQ=255)

           cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCT‑CCCCCCCCGTTCGGTCCGCGAgctgcggg                                                                                                                                          >  2:3405420/1‑147 (MQ=255)

             cTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGGAATCGAGCTGCTGAcc                                                                                                                                        >  2:755466/1‑150 (MQ=255)

              tGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCg                                                                                                                                       <  1:2404711/150‑1 (MQ=255)

              gtcgcATCTACCAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCg                                                                                                                                       >  1:3682943/3‑150 (MQ=255)

               gcgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGTCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGc                                                                                                                                      <  1:616492/150‑1 (MQ=255)

                cgcATCTGGTAGATGCGACTGTTCTATCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                                                                                                     <  1:1454556/150‑1 (MQ=255)

                       ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                              <  1:1099527/150‑1 (MQ=255)

                       ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAAGCGAGCTGCTGACCGCTGTTATcc                                                                                                                              >  2:365038/1‑150 (MQ=255)

                         tAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCATTGTGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCg                                                                                                                             <  2:25175/150‑1 (MQ=255)

                         tAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGc                                                                                                                            <  1:3071/150‑1 (MQ=255)

                         tAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGc                                                                                                                            >  2:3174779/1‑150 (MQ=255)

                         tAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGc                                                                                                                            <  2:3265817/150‑1 (MQ=255)

                         gaGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGc                                                                                                                            <  2:1590784/149‑1 (MQ=255)

                           gATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                                                          <  2:696626/150‑1 (MQ=255)

                              gCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGG‑CCGCCGCCGTGCGGGGATCGGGCTGCGGGCGGGTGGTATCCGCCGGcg                                                                                                                       >  1:1659893/1‑149 (MQ=255)

                                gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCCCCT‑CCGCCCCCGTGCGGGGAGCGGGCTGCTGGCGGCGGCGGTCCGCCAGcggc                                                                                                                     >  2:2590733/1‑147 (MQ=255)

                                 aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCt                                                                                                                    <  1:481856/150‑1 (MQ=255)

                                    gTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGAc                                                                                                                 <  1:1253311/150‑1 (MQ=255)

                                    gTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGAc                                                                                                                 <  1:50835/150‑1 (MQ=255)

                                      tCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                               <  2:805861/150‑1 (MQ=255)

                                      tCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGGACGGGCGGCGGGCTGCGGACCGCTGCTCTCCGCCAGCCTCTGCCGc                                                                                                               >  2:3338944/1‑150 (MQ=255)

                                        tAACGGTAGTTGCGGCAACAAGTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTCATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                             <  2:1593854/150‑1 (MQ=255)

                                             gTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCACCAGCATCTGACGCCAGTACt                                                                                                        <  2:3682943/150‑1 (MQ=255)

                                                gTTGCGGCAACACTTCAATGAGCTATTTCTGTGGATTCGCCGCTATCTGCGACAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACTTTAATCGATATGCTGACCGCTGATATCAGCCATCATCTGACGCCAGTACTCag                                                                                                     <  2:700438/150‑2 (MQ=255)

                                                gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAc                                                                                                     <  2:2100159/150‑1 (MQ=255)

                                                gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCAc                                                                                                     <  1:3370112/150‑1 (MQ=255)

                                                     ggCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGGGCGGGCCTCGAGCTGCTGACCGCTGATATCCGCCAGCGTCTGACGCCAGTGCTCGCTAcgc                                                                                                >  2:355941/1‑150 (MQ=255)

                                                       caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACgcgc                                                                                              <  2:462303/150‑1 (MQ=255)

                                                       caacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACgcgc                                                                                              <  2:2106322/150‑1 (MQ=255)

                                                            aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCT‑CCCCCCCCGTACGGTACTCGAGCTGCTGGCCGCTGATATCCGCCCGCAGCTGACGCCCGTACTCACTACGcgccccgc                                                                                         >  2:3407060/1‑150 (MQ=255)

                                                            aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCG‑CCCCCCCCGGGCGGTGAGCGAGGTGCTGACCGCGGGTGGCGGCCCGCGTCTGGCGGCCGGACTCGCTCCGcgcgccgc                                                                                         >  1:1768859/1‑150 (MQ=255)

                                                              tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgccc                                                                                       <  1:419534/150‑1 (MQ=255)

                                                              tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgccc                                                                                       <  1:1533236/150‑1 (MQ=255)

                                                                aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAG‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCAGCTGACGCCAGTACTCACTACGCGCCCCGccccc                                                                                     >  2:1419927/1‑150 (MQ=255)

                                                                  tGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCaa                                                                                   <  2:1067742/150‑1 (MQ=255)

                                                                  tGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCaa                                                                                   <  2:3727419/150‑1 (MQ=255)

                                                                    aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTACTCGAGCTGCTGACCGCGGATATCCGCCAGCATCTGACGCCAGTACTCACTCCGCGCCCCGccccccagc                                                                                 >  2:2832480/1‑146 (MQ=255)

                                                                    aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATc                                                                                 <  2:1349416/150‑1 (MQ=255)

                                                                      cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTTATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCaa                                                                               <  2:214923/150‑1 (MQ=255)

                                                                      cGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCaa                                                                               <  1:3745372/150‑1 (MQ=255)

                                                                       gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                                              <  2:2874564/150‑1 (MQ=255)

                                                                          ttCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGGACGGTAAGCGCGCTGCTGACCGCGGGTATCCGCCCGCATCTGACGCCGGTGCGCACTACGCGGCCCGCCCCCCGGCAACGTa                                                                           >  2:3285583/1‑150 (MQ=255)

                                                                           tctcTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTaa                                                                          <  1:1657408/150‑1 (MQ=255)

                                                                           tctcTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTaa                                                                          <  2:2423740/150‑1 (MQ=255)

                                                                           tctcTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGACAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTaa                                                                          <  2:2725196/150‑1 (MQ=255)

                                                                                       ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                              <  1:3337268/150‑1 (MQ=255)

                                                                                       ccgcGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGtt                                                              <  2:1969456/150‑1 (MQ=255)

                                                                                        cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                             <  2:3100302/150‑1 (MQ=255)

                                                                                         gcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑‑CCCCCCCGTACGGTCATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTaa                                                           <  1:700458/150‑1 (MQ=255)

                                                                                             tCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAc                                                        <  1:772761/150‑1 (MQ=255)

                                                                                              cTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTGATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAcc                                                       <  2:2902363/150‑1 (MQ=255)

                                                                                                 cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCGGGCCGCTGGTATCCGCCCGCCTCTGACGCCAGGTCTCACTACGCGCCCCGCCCCCCCGCTACGTAACACTTTGCGGGTTAATACCGCa                                                    >  2:3941134/1‑150 (MQ=255)

                                                                                                  gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCT‑CCCCCCCCGTACGGTCATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCCGTCCTCACTACGCGCCCCGCCCCCAATCAACGGACCACTTTGCGGGTGAATACCGccg                                                   >  2:1750209/1‑148 (MQ=255)

                                                                                                         ttGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCa                                            <  1:317185/150‑1 (MQ=255)

                                                                                                            ctgcGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGa                                         <  2:3034243/150‑1 (MQ=255)

                                                                                                               cgcCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACAt                                      >  2:171864/1‑150 (MQ=255)

                                                                                                                 cccAGTGCTGGCCCCACCT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGGTCTCCGCCCGCAGCTGGCGCCCGTCCTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATcc                                    >  1:1353742/1‑150 (MQ=255)

                                                                                                                  ccAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGTCCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCa                                   >  2:3182556/1‑150 (MQ=255)

                                                                                                                  ccAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCa                                   <  1:1315605/150‑1 (MQ=255)

                                                                                                                   cAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                                  <  1:2532042/150‑1 (MQ=255)

                                                                                                                   cAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                                  <  2:1996488/150‑1 (MQ=255)

                                                                                                                    aGTGCTGGCCCCACAT‑CCCCCCCCGTATGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                                 <  2:410143/150‑1 (MQ=255)

                                                                                                                    aGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                                 <  2:1956607/150‑1 (MQ=255)

                                                                                                                      tGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGcc                               <  1:1840883/150‑1 (MQ=255)

                                                                                                                        cTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCAGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTTATACCGCAGGCATGCACGACATCCATGCCaa                             >  1:2198261/1‑149 (MQ=255)

                                                                                                                          ggCCCCCCAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGAAGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCt                           <  2:1211703/150‑1 (MQ=255)

                                                                                                                            ccccACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGc                         >  2:2317173/1‑150 (MQ=255)

                                                                                                                             cccACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGcc                        >  2:2693632/1‑150 (MQ=255)

                                                                                                                             cccACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGcc                        >  2:1496695/1‑150 (MQ=255)

                                                                                                                             cccACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAACAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGcc                        >  1:3192884/1‑150 (MQ=255)

                                                                                                                              ccACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa                       <  1:2281564/150‑1 (MQ=255)

                                                                                                                              ccAAAT‑ACCCCCCCGTACGGTAATATAGCTGCACACCGCAGATATCCGCCAGCATCTGACGCCAGTACTAACAACGCGCCCCGCCCCCAAGCAACGTAACACTTTGAGGTTTAATAACGCACGAAGGCCCTACATCCATGCCATCTGccc                       >  2:3499665/1‑149 (MQ=255)

                                                                                                                                  aT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcgc                   <  1:1833286/150‑1 (MQ=255)

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

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

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

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

                                                                                                                                                 ggTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCtt     <  1:2149137/150‑1 (MQ=255)

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

                                                                                                                                                    aaTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCa  <  1:3108238/150‑1 (MQ=255)

                                                                                                                                                    aaTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCa  <  1:1196448/150‑1 (MQ=255)

AAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCCA  >  NC_000913/3727979‑3728275

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