BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,728,125	T→G	5.4%	Y416S (TAC→TCC)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,125	0	T	G	5.4%	1646.9 / 12.8	578	Y416S (TAC→TCC)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base T (164/364); new base G (43/5); total (207/371)
Fisher's exact test for biased strand distribution p-value = 1.00e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 590 total aligned reads displayed.

ATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727977‑3728274

aTAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGtc                                                                                                                                                      >  2:1478692/1‑149 (MQ=255)

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

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

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

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

        cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCGCCCCGTGCGGTGATCGGgcg                                                                                                                                              >  1:5494671/1‑149 (MQ=255)

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

          tAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCGCCCCCGGCCGGTAAGCGGgctgg                                                                                                                                            >  2:469114/1‑149 (MQ=255)

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

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

                      tCTTGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGata                                                                                                                                <  2:5481623/150‑1 (MQ=255)

                        tGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATc                                                                                                                              <  1:3774788/150‑1 (MQ=255)

                         ggTAGATGCGACTGTTCTGACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                             <  2:1652225/150‑1 (MQ=255)

                         ggTAGATGCGACTGTTCTCACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                             <  2:3589386/150‑1 (MQ=255)

                         ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                             <  1:651364/150‑1 (MQ=255)

                         ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                             <  1:2856812/150‑1 (MQ=255)

                         ggTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATcc                                                                                                                             <  1:1256/150‑1 (MQ=255)

                          gTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGGCTTCGGGCTGCTGGCCGCGGATGTCCg                                                                                                                            >  2:4932474/1‑150 (MQ=255)

                           tAGTTGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGc                                                                                                                           <  2:2711387/150‑1 (MQ=255)

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

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

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

                              aTGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAg                                                                                                                        <  1:1321531/150‑1 (MQ=255)

                                  gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGTAAGCGGGCTGCTGACCGCTGAGATCCGCCAGCGTc                                                                                                                    >  2:4980448/1‑150 (MQ=255)

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

                                       ttCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGACGGTAAGCGAGCTGCTGACCGCTGATATCCGCCAGCCTCTGGCg                                                                                                               >  2:2335196/1‑150 (MQ=255)

                                        tCTAAGGGTTGTTGCGGCAAAATTTCAATTTGCGATTTCTGTGTTTTCGCGGGTTTCTGCGCCAGCCTTGCTGCGCCCTGTGCTGGCCCCACATCCCCCCCCTTAGGTTAATCGAGGTGGTGACCGCTTATTTCCGCCCGCATCTGACGc                                                                                                              <  2:4401183/150‑1 (MQ=255)

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

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

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

                                         cTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGAGCTGCTGACCGCTGATATCCGCCAGCAGCTGGCGcc                                                                                                             >  2:4740949/1‑150 (MQ=255)

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

                                          tAACGGTAGGTGCGGCAACAATTCAATGAGCGCTTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                            <  2:5292404/150‑1 (MQ=255)

                                            aCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGTAATCGAGCTGCGGACCGCTGATAGCCGCCAGCGGCTGACGCCCGt                                                                                                          >  2:2034754/1‑150 (MQ=255)

                                             cGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTAc                                                                                                        >  2:4105962/1‑150 (MQ=255)

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

                                                  gTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTGATCGAGCTGCTGGCCGCTGAGCTCCGCCCGCGTCTGACGCCCGTGCTCAc                                                                                                    >  2:871088/1‑150 (MQ=255)

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

                                                     gCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAc                                                                                                 <  1:355986/150‑1 (MQ=255)

                                                          aacaaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGGAATCGAGCTGGTGACCGCTGGTATCCGCCCGCATCTGCCGCCAGTACGCGCTTGGcgcc                                                                                            >  2:3521741/1‑150 (MQ=255)

                                                             aaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGGGGGGGACTCGGGCTGCGGGGCGCGGGTGTCGGCCAGCATCGGGCGCGGGTGCGCACTCCGcgccccg                                                                                         >  1:4083364/1‑150 (MQ=255)

                                                                 cAATGAGCGATTTCTCTGGATTCGCCGCGCTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCgccccgcccc                                                                                     <  2:3759855/150‑1 (MQ=255)

                                                                  aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGccccc                                                                                    >  2:1108129/1‑150 (MQ=255)

                                                                  aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGTGCTGCTGGCCGCTGATATCCGCCAGCAGCTGACGCCAGTACTCACTACGCGCCCCGccccc                                                                                    >  2:1602161/1‑150 (MQ=255)

                                                                  aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGGTATCCGCCAGCAGCGGACGCCAGGACTCACGACGCGCCCCGccccc                                                                                    >  1:2742225/1‑150 (MQ=255)

                                                                    tGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCaa                                                                                  <  1:3243365/150‑1 (MQ=255)

                                                                     gagagaTTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt                                                                                 <  2:388985/150‑1 (MQ=255)

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

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

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

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

                                                                             tctcTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGAAGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTaa                                                                         <  1:5664229/150‑1 (MQ=255)

                                                                             tctcTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTACGGGAATCGAGCGGCTGACCGCTGCTATCCGCCAGCCTCTGACGCCAGGACTCACTACGCGCCCCGCCCCCACTCACCggga                                                                         >  2:4810871/1‑147 (MQ=255)

                                                                                ctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGTACTCGTGCTGCTGCCCGCGGTGGTCCGCCCGCCGGTGGCGCCAGTCCTCCCTGCGCGCCCCGCGCCCACTCAACggagccc                                                                      >  2:10363/1‑144 (MQ=255)

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

                                                                                      tCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGGATCGTGCTGCTGACCGCTGGTGTCCGCCCGCATCTGACGCCAGTTCTCACTACGCGCCCCGCCCCCAGTCGACGTGGCACGTTGCg                                                                >  2:502156/1‑150 (MQ=255)

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

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

                                                                                              aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTGCGGTAATCGCGCTGCTGACCGCGGGTATCCGCCAGCATCTGGCGCCAGTACTCACTACGCGCCCCGCCCCCCATCCACGTCACCCTTTGCGGTTTACTa                                                        >  2:543119/1‑150 (MQ=255)

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

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

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

                                                                                                     ccAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCaaa                                                 <  2:1260304/150‑3 (MQ=255)

                                                                                                       aGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCATCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCa                                               <  2:5159892/150‑1 (MQ=255)

                                                                                                          cTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGc                                            <  1:4732969/150‑1 (MQ=255)

                                                                                                          cTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGc                                            <  1:5152856/150‑1 (MQ=255)

                                                                                                           ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTCCGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCGGGCATGCa                                           >  2:4527666/1‑150 (MQ=255)

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

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

                                                                                                               tgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGAc                                       <  2:3131947/150‑1 (MQ=255)

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

                                                                                                                gcgcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCATCGTAACACTTTGCGGTTTTATACCGCAGGCATGCACGACa                                      <  2:613446/150‑1 (MQ=255)

                                                                                                                   cccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATcc                                   >  2:5502090/1‑150 (MQ=255)

                                                                                                                   cccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATcc                                   <  2:5603156/150‑1 (MQ=255)

                                                                                                                      cgTGCTGGCCCCACATCCCCCCCCGTCCGGTAAGCGAGCTGCTGCCCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCACCGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                                <  2:2409401/149‑1 (MQ=255)

                                                                                                                       gTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCCACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGc                               <  2:248338/150‑1 (MQ=255)

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

                                                                                                                           tGGCCCCACATCCCCCCCCGTTCGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGGACTCACTACGCGCCCCGCCCCCAATCAACGTACCACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATc                           >  2:534702/1‑150 (MQ=255)

                                                                                                                           tGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCCTCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTTATACCGCAGGCAGGCACGACATCCATGCCATc                           >  1:4779328/1‑150 (MQ=255)

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

                                                                                                                             gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTCACACATTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg                         >  1:2550163/1‑150 (MQ=255)

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

                                                                                                                             gCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAACAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTg                         >  1:2354914/1‑150 (MQ=255)

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

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

                                                                                                                              ccccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACAACGCGCCCCGCCCCCAATCAACGTACCACTTTGCGGTGTAATAACGCAGGCATGCACGACATCCATGCCATCAGc                        >  1:4305905/1‑150 (MQ=255)

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

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

                                                                                                                                ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCAGCCa                      >  2:2797708/1‑150 (MQ=255)

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

                                                                                                                                         cccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCCGCGCATAAc             >  2:3271876/1‑150 (MQ=255)

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

                                                                                                                                          ccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAcc            <  1:2700061/150‑1 (MQ=255)

                                                                                                                                          ccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAcc            <  1:5070717/150‑1 (MQ=255)

                                                                                                                                          ccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAcc            <  1:77831/150‑1 (MQ=255)

                                                                                                                                            ccGTACGGGACTCGAGCTGCTGACCGCTGATATCCGCCAGCATCGGACGCCAGTACTCACGACGCGCCCCGCCCCCAATCAACGTAACCCTTTGCGGTTTAAGACCGCAGGCATGCCCGACCTCCATGCCATCTGCCAGCGCATAAcccc          <  2:4858718/150‑2 (MQ=255)

                                                                                                                                              gTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCATCATCTGACGCCAGTACTCACTACGCGCCCCGCCACCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGTCATCCATGCCATCTGCCAGCGCATAACCCACg        <  2:807056/150‑1 (MQ=255)

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

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

                                                                                                                                                    tAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGTAGGCATGCACGACATCCATGCCATCTGCCAg                     >  1:836350/1‑131 (MQ=255)

ATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGCCTTCC  >  NC_000913/3727977‑3728274

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