BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,125	T→G	6.2%	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	6.2%	1886.8 / 21.9	642	Y416S (TAC→TCC)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base T (187/397); new base G (50/7); total (238/404)
Fisher's exact test for biased strand distribution p-value = 1.30e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 647 total aligned reads displayed.

TAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727978‑3728274

tAAGCGGCATAACGCTGCGCATCTGGTAGTTGCGACTGTTCTAACGGTTGTTGCGGCAACAATTCAATGAGCGATTTCTCTGTATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTaa                                                                                                                                                     <  2:3650566/150‑1 (MQ=255)

 aaGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGGCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAAt                                                                                                                                                    <  2:3771509/150‑1 (MQ=255)

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

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

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

        aTAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGGCATCGAgctg                                                                                                                                             >  2:6600302/1‑150 (MQ=255)

        aTAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGGGCGGTAAtggggcgg                                                                                                                                             >  1:5965687/1‑143 (MQ=255)

         tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAG‑CCCCCCCCGTACGGGAAGCGAgctgc                                                                                                                                            >  2:1536315/1‑150 (MQ=255)

          aaCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAgctgct                                                                                                                                           <  2:1729511/150‑1 (MQ=255)

             gCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTCCGGGCGTCGGGctgctggc                                                                                                                                        >  2:2895020/1‑148 (MQ=255)

              cTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCACCAATTCAATGATATCTTGCTCTGGAGGCGCCTCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCT‑CCCCCCCCGTACGGTAGTCGAGCTGCGGAcc                                                                                                                                       >  2:5763945/1‑150 (MQ=255)

              cTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGAcc                                                                                                                                       >  1:1815032/1‑150 (MQ=255)

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

                 cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCT‑CCCCCCCCGTACGGGAATCGAGCTGCTGACCGCt                                                                                                                                    <  2:5797548/150‑1 (MQ=255)

                   cATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGa                                                                                                                                  <  2:264858/150‑1 (MQ=255)

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

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

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

                            gATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAAGGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCa                                                                                                                         <  2:4209107/150‑1 (MQ=255)

                                 gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATc                                                                                                                    <  1:5671825/150‑1 (MQ=255)

                                  aCTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCT‑CCCCCCCCGGGCGGGGGGCGAGCTGGTGGCCGCTGATCGCCGCCCGCGTCt                                                                                                                   >  2:2939519/1‑150 (MQ=255)

                                    tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                                                                 <  1:3746096/150‑1 (MQ=255)

                                    tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                                                                 <  1:5136559/150‑1 (MQ=255)

                                    tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                                                                 <  1:6242577/150‑1 (MQ=255)

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

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

                                         tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CGGGCCCCGTAGGGTAATGGAGCTGGTGAGGGGTGATTTCCGCGAGCGTCTGACggcg                                                                                                            >  2:4643701/1‑147 (MQ=255)

                                          aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGGACGGGAGTCGAGCTGCTGGCCGCGGGTATCCGCCCGCCTCTGACGCCAg                                                                                                           >  2:5383624/1‑150 (MQ=255)

                                                     cGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGGACGGTACTCGAGCTGCTGGCCGCTGAGATCCGCCCGCATCTGTCGCCCGTCCTCagtccg                                                                                                >  2:641477/1‑145 (MQ=255)

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

                                                           caaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGT‑CCCGCCCCGGGCGGGAGTCGGGGTGGTGAGCGGTGTTGGGCGCGGGGGTCTGACGCCAGTGCGCTGTACGcgcgcc                                                                                          >  1:5513951/1‑150 (MQ=255)

                                                            aaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCG‑CCCCCCCCGTCCGGGACTCGCGCGGCGGACCGCGGCGCTCCGCCCGCCGCGGCCGCCCGTCCGCACGCCGcgccccg                                                                                         >  2:890162/1‑150 (MQ=255)

                                                             aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGG‑CCCCCCCCGGGCGGGAATCGGGCTGGGGACCGCTGTTGTCCGCCGGCCTCTGACGCCCGGACGCTCTAGGcgccccgc                                                                                        >  1:4407984/1‑150 (MQ=255)

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

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

                                                                 aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGccccc                                                                                    <  1:2841900/150‑1 (MQ=255)

                                                                   tgtgCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACTT‑CCCGCCCCGTTCGTTATTCGTGCTGCTGACCGCTGATATCCGCGAGCATCTGACGCCAGTACTCACTACGCGGCCCGCCCCCaa                                                                                  <  2:402248/147‑1 (MQ=255)

                                                                    gAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACGT‑CCCCCCCCGGACGGGCGTGGAGCTGCTGTCGGCTGGTAGCCGCCCGCATCGGACGGCCGTCCGCACTACGCGCCCCGcccccggg                                                                                 >  1:3904200/1‑147 (MQ=255)

                                                                     aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAG‑CCCCCCCCGGCCGGGCATCGAGCTGCTGCCCGCTGCTCTCCGCCCGCATCTGCCGCCCGTACTCGCTACGCGCCCCGCCCCCaagc                                                                                >  2:4809500/1‑148 (MQ=255)

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

                                                                        gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCG‑CCCCCCCCGGACGGTAATCGAGCTGCTGGCCGCTGTTCTCCGCCCGCATCTGACGGCGGTCCTCACTAGGCGCCCCGcgcccgctgcgc                                                                             >  2:5264852/1‑143 (MQ=255)

                                                                        gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAACAGCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAAc                                                                             >  2:3417954/1‑150 (MQ=255)

                                                                              tctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTACTCGCGCTGCTGGCCGCTGATCTCCGCCCGCATCTGACGCCCGTACTCACTACGCGCCCCGCCCCCAATCAACGGAacc                                                                       >  2:3924994/1‑149 (MQ=255)

                                                                                  gATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACttt                                                                   <  2:786561/150‑1 (MQ=255)

                                                                                     tCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCg                                                                <  2:190526/150‑1 (MQ=255)

                                                                                     tCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCg                                                                <  2:2998318/150‑1 (MQ=255)

                                                                                     tCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCg                                                                <  1:6380951/150‑1 (MQ=255)

                                                                                        ttccgATCT‑CGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                            <  1:4847694/147‑1 (MQ=255)

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

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

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

                                                                                         cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                            <  2:4402612/150‑1 (MQ=255)

                                                                                          gcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                            <  1:2527911/150‑1 (MQ=255)

                                                                                          gcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTCCGGTAAGCGGGCGGCGGACCGCGGGGAGCCGCCAGCAGCTGAGGCCCGGACTCACGGCGCGCCCCGCCCCCACGCACCGTACCACTTGGCgggtga                                                           >  2:4879605/1‑146 (MQ=255)

                                                                                            gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAAt                                                         <  1:5331967/150‑1 (MQ=255)

                                                                                             aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATa                                                        <  1:19150/150‑1 (MQ=255)

                                                                                             aTCTGCGCCAGCCTTGCGGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATa                                                        <  2:5200831/150‑1 (MQ=255)

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

                                                                                               cTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAcc                                                      <  1:700553/150‑1 (MQ=255)

                                                                                                tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CGCCCCCCGTTCGGTAATCGTGCTGCTGACCGCTGATATCCGCCCGCATCTGACGCCAGGACTCACTACGCGCCCCGCCGCCACTCACCGTACCACGTTGCGGTTTAAGACCg                                                     >  2:1194590/1‑150 (MQ=255)

                                                                                                tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTCCGGTCCTCGGGGTGCTGACCGCTGGTCTCCGCCAGCGGCTGCCGCCGGTGCGCGCGCCGCGCCCCGCCCCCGAGCGGCGTGACGCTTTGCGGGTGAATAGCg                                                     >  2:2398145/1‑150 (MQ=255)

                                                                                                tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCGTCTGACGCCCGTACTCACTACGCGCCCCGCCCCCACTCAACGTACCAGTTTGCGGTGTAATACCg                                                     >  2:5596778/1‑150 (MQ=255)

                                                                                                tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg                                                     <  1:992915/150‑1 (MQ=255)

                                                                                                tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg                                                     <  2:2485387/150‑1 (MQ=255)

                                                                                                tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg                                                     <  2:6281218/150‑1 (MQ=255)

                                                                                                 gcgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTGCGGTAATCGAGCTGCTGGCCGCTGTTATCCGCCCGCATCTGACGCGAGGACTCACTACGCGCCCCGCCCCCTAGCCACGTAACGCTTTGCGGTTTATTACCGc                                                    >  2:1986539/1‑150 (MQ=255)

                                                                                                  cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGTTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa                                                   <  1:2385858/150‑1 (MQ=255)

                                                                                                  cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACTGCa                                                   <  1:6149247/150‑1 (MQ=255)

                                                                                                  cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa                                                   <  2:4294709/150‑1 (MQ=255)

                                                                                                      aGCCTTGCTGCGCCCAGTGCTGGCCCCACCT‑CCCCCCCCGTACGGTAGTCGGGCTGCTGGCCGCTGATATCCGCCAGCATCTGGCGCCTGTACTCACTACGCGCCCCGCCCCCAATCGCCGTGACGCTTTGGGGTTTTATGCCGCCGGcg                                               >  2:1246569/1‑149 (MQ=255)

                                                                                                         cTTGCTGCGCCCAGTGCTGGCCCCACCT‑CCCCCCCCGTCCGGTCAGCGCGCTGCGGTCCGCGGATAGCCGCCCGCAGCGGGCGCCGGGACTCCCTTCGCGCCCCGCCCCCCAGCGACGTACCCGGGTGCGGGGTTATACCGCAGGcgggc                                            >  2:3346527/1‑146 (MQ=255)

                                                                                                            gctgcGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACg                                         <  2:4135649/150‑1 (MQ=255)

                                                                                                               gcgcCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTTACCGCTTATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACa                                      <  1:130982/150‑1 (MQ=255)

                                                                                                                 gcCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATc                                    <  1:1604729/150‑1 (MQ=255)

                                                                                                                  cccAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAAGCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCCGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATcc                                   >  1:5835961/1‑150 (MQ=255)

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

                                                                                                                   ccAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCGGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACGTTGCGGTTTAATACCGCAGGCATGCACGACATCCa                                  >  2:4573431/1‑150 (MQ=255)

                                                                                                                    cAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGGCCGCTGATTTCCGCCAGCATCGGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTGACACGTTGCGGTTTAATACCGCAGGCATGCACGccatccat                                 >  1:475781/1‑150 (MQ=255)

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

                                                                                                                    cAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAt                                 >  2:1379394/1‑150 (MQ=255)

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

                                                                                                                     aGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAGCCACGTACCACTTTGCGGTTTAATACCGCCGGCATGCCCGACCTCCATg                                >  2:4732546/1‑150 (MQ=255)

                                                                                                                      gTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGc                               <  2:3071544/150‑1 (MQ=255)

                                                                                                                         cTGGCCCCACCG‑CCCCCCCCGGACGGTACTCGGGCTGCTGACCGCTGATAGCCGCCCGCATCTGGCGCCCGTACTCACTACGCGCCCCGCCCCCCCGCAACGGAACCCGTTGCGGGGGACTACCGCAGGCATGCACGACAGCCAGGcccg                            >  2:2171811/1‑148 (MQ=255)

                                                                                                                         cTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCAt                            <  1:6061135/150‑1 (MQ=255)

                                                                                                                            gCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg                         >  2:3158908/1‑150 (MQ=255)

                                                                                                                             ccccACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCTCTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGc                        <  2:3684685/150‑1 (MQ=255)

                                                                                                                              cccACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAGCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCAAGCCATCTGcc                       >  2:6071332/1‑150 (MQ=255)

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

                                                                                                                                  caT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgcg                   <  1:3148961/150‑1 (MQ=255)

                                                                                                                                   cagccccccccGTACGGTAAGCGAGCTGCGGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACCCTTTGCGGTTTAATACCGCAGGCCTGCACGACATCCATGCCATCTGCCAgcg                   <  2:2098814/147‑1 (MQ=255)

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

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

                                                                                                                                         cccccGTACGGTACTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAc             <  2:4439381/150‑1 (MQ=255)

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

                                                                                                                                            ccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCa          >  1:228798/1‑150 (MQ=255)

                                                                                                                                              gTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACg        <  2:6352570/150‑1 (MQ=255)

                                                                                                                                              gTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCACTCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACg        <  1:6493877/150‑1 (MQ=255)

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

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

TAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727978‑3728274

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