BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,623,143	A→T	70.7%	L170* (TTG→TAG)	wzy ←	O‑antigen polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,623,143	0	A	T	70.7%	567.8 / 302.0	358	L170* (TTG→TAG)	wzy	O‑antigen polymerase
Reads supporting (aligned to +/- strand): ref base A (58/47); new base T (149/104); total (207/151)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 360 total aligned reads displayed.

ACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAG  >  NC_002947/5622994‑5623291

acaACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCa                                                                                                                                                      >  2:3470803/1‑150 (MQ=255)

 caACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCaa                                                                                                                                                     <  1:7457977/150‑1 (MQ=255)

  aaCGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAg                                                                                                                                                    >  1:2830191/1‑150 (MQ=255)

    cGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGc                                                                                                                                                  >  1:850334/1‑150 (MQ=255)

    cGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGc                                                                                                                                                  <  1:670565/150‑1 (MQ=255)

        gcCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGaa                                                                                                                                              >  1:1885176/1‑150 (MQ=255)

        gcCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGaa                                                                                                                                              >  1:6204058/1‑150 (MQ=255)

            tCACCACGGTAATCACCAGTGCGAGAACCGCACCACGGATCTGGCTGAGCTTGTCTTATGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCaa                                                                                                                                          <  2:656712/150‑1 (MQ=255)

            tCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCaa                                                                                                                                          <  2:5310744/150‑1 (MQ=255)

               ccacGGTGATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGACCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGa                                                                                                                                       <  1:5563784/150‑1 (MQ=255)

               ccacGGTAATAACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGTCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCAGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGa                                                                                                                                       <  1:2437535/150‑1 (MQ=255)

                cacGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGAc                                                                                                                                      >  1:5608868/1‑150 (MQ=255)

                    gTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATaa                                                                                                                                  >  1:2639011/1‑150 (MQ=255)

                     tAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAg                                                                                                                                 >  1:6374744/1‑150 (MQ=255)

                        tCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGccc                                                                                                                              >  2:5829717/1‑150 (MQ=255)

                        tCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGACCCAATGACATCAGccc                                                                                                                              >  1:4614294/1‑150 (MQ=255)

                         cACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCGGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGcccc                                                                                                                             >  1:5283221/1‑150 (MQ=255)

                           ccAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCa                                                                                                                           >  2:3820425/1‑150 (MQ=255)

                            cAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAg                                                                                                                          >  1:5128654/1‑150 (MQ=255)

                              gTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCTAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGa                                                                                                                        <  1:5993573/150‑1 (MQ=255)

                              gTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAACCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGa                                                                                                                        >  1:4004757/1‑150 (MQ=255)

                               tGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGAt                                                                                                                       <  2:3333055/150‑1 (MQ=255)

                               tGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGAt                                                                                                                       <  2:5660348/150‑1 (MQ=255)

                                gCCAGAACCGCACCACGGCTCTGTCTGAGCATGGCGAAGGCCCACAGGCAAGCCAGTGCCTCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGTATCGTATAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGAtt                                                                                                                      <  2:5425566/150‑2 (MQ=255)

                                 ccAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATgg                                                                                                                     >  1:1005640/1‑150 (MQ=255)

                                 ccAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATgg                                                                                                                     >  1:1543501/1‑150 (MQ=255)

                                  cAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATggg                                                                                                                    >  2:7705439/1‑150 (MQ=255)

                                  cAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATggg                                                                                                                    <  1:215761/150‑1 (MQ=255)

                                   aGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATgggg                                                                                                                   >  2:223279/1‑150 (MQ=255)

                                      aCCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTgg                                                                                                                <  1:4448490/150‑1 (MQ=255)

                                      aCCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTgg                                                                                                                <  2:3840018/150‑1 (MQ=255)

                                        cGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGa                                                                                                              >  2:4945961/1‑150 (MQ=255)

                                        cGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGa                                                                                                              >  2:5965654/1‑150 (MQ=255)

                                         gCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGat                                                                                                             >  1:5689521/1‑150 (MQ=255)

                                         gCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGat                                                                                                             <  2:7163537/150‑1 (MQ=255)

                                         gCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGat                                                                                                             >  1:7890086/1‑150 (MQ=255)

                                             cacGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATAtct                                                                                                         >  1:6468023/1‑150 (MQ=255)

                                                ggCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCGCGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGc                                                                                                      <  2:4346571/150‑1 (MQ=255)

                                                ggCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGc                                                                                                      <  1:6338523/150‑1 (MQ=255)

                                                     tGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATa                                                                                                 <  1:4838706/150‑1 (MQ=255)

                                                       gCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATAcc                                                                                               >  1:939473/1‑150 (MQ=255)

                                                         tGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCg                                                                                             >  1:2457762/1‑150 (MQ=255)

                                                         tGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCg                                                                                             >  1:3365238/1‑150 (MQ=255)

                                                         tGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCg                                                                                             >  2:8668257/1‑150 (MQ=255)

                                                                  cGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGa                                                                                    <  1:477253/150‑1 (MQ=255)

                                                                   gAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTCGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGCTACCCGCCAGCCTGac                                                                                   >  2:4409427/1‑149 (MQ=255)

                                                                      ggCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAg                                                                                <  2:485022/150‑1 (MQ=255)

                                                                       gCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGc                                                                               >  2:8753767/1‑150 (MQ=255)

                                                                         ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaa                                                                             >  2:2966627/1‑150 (MQ=255)

                                                                          cccAGGCAAGCCAGTGTCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaaa                                                                            >  2:7484982/1‑150 (MQ=255)

                                                                          cccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCCCGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTCGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaaa                                                                            >  2:6889834/1‑150 (MQ=255)

                                                                          cccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaaa                                                                            >  2:3927267/1‑150 (MQ=255)

                                                                          cccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaaa                                                                            >  1:2977207/1‑150 (MQ=255)

                                                                          cccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaaa                                                                            <  2:1906993/150‑1 (MQ=255)

                                                                          cccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaaa                                                                            <  2:8553803/150‑1 (MQ=255)

                                                                          cccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaaa                                                                            >  2:6451222/1‑150 (MQ=255)

                                                                          cccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaaa                                                                            >  1:4122223/1‑150 (MQ=255)

                                                                           ccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAg                                                                           >  2:7609267/1‑150 (MQ=255)

                                                                           ccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAg                                                                           <  2:80844/150‑1 (MQ=255)

                                                                             aGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGgc                                                                         <  1:8668257/150‑1 (MQ=255)

                                                                                 aaGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCAc                                                                     >  1:3794361/1‑150 (MQ=255)

                                                                                        tGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCAc                                                              <  1:773272/150‑1 (MQ=255)

                                                                                          ccgccAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCAccc                                                            <  1:4328510/150‑1 (MQ=255)

                                                                                            gccagccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCa                                                          >  1:1253864/1‑150 (MQ=255)

                                                                                            gccagccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCa                                                          <  1:3185284/150‑1 (MQ=255)

                                                                                              cagccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATa                                                        <  2:5197906/150‑1 (MQ=255)

                                                                                                gccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATACCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGa                                                      >  2:1749063/1‑150 (MQ=255)

                                                                                                gccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATACCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGa                                                      >  2:6172892/1‑150 (MQ=255)

                                                                                                 ccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGaa                                                     <  1:8784252/150‑1 (MQ=255)

                                                                                                  cacaACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAAc                                                    >  2:1448241/1‑150 (MQ=255)

                                                                                                  cacaACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAAc                                                    >  2:2605433/1‑150 (MQ=255)

                                                                                                   acaacaACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACt                                                   >  1:7098178/1‑150 (MQ=255)

                                                                                                   acaacaACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACt                                                   >  1:8234598/1‑150 (MQ=255)

                                                                                                    caacaaCTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACtt                                                  <  2:4226690/150‑1 (MQ=255)

                                                                                                     aacaacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCTATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCCTCGAACTtg                                                 >  2:6949783/1‑150 (MQ=255)

                                                                                                      acaacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtga                                                <  2:315080/150‑1 (MQ=255)

                                                                                                         tctacgaGCCACCTCGGCGACGTGGCTTGCAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtgatga                                             >  1:5810448/7‑150 (MQ=255)

                                                                                                          cTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACa                                                                                                                                     >  2:229628/1‑61 (MQ=255)

                                                                                                           tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGatgatg                                           <  1:7589730/150‑1 (MQ=255)

                                                                                                           tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGatgatg                                           <  2:8364922/150‑1 (MQ=255)

                                                                                                             cAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGga                                         <  1:5250799/150‑1 (MQ=255)

                                                                                                              gaGCCACCTCGGCGACGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGgag                                        >  2:2940633/2‑150 (MQ=255)

                                                                                                               aGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCATGCTGCACCCCATAGAACTTGATGATGgaga                                       <  2:2566181/150‑1 (MQ=255)

                                                                                                                gCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGAc                                      <  1:3688193/150‑1 (MQ=255)

                                                                                                                gCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGAc                                      >  2:5412450/1‑150 (MQ=255)

                                                                                                                 ccACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCAATAGAACTTGATGATGGAGAcc                                     >  1:2022044/1‑150 (MQ=255)

                                                                                                                           gAGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGa                           >  2:3665098/1‑150 (MQ=255)

                                                                                                                            aGGTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGAt                          >  1:1649010/1‑150 (MQ=255)

                                                                                                                             ggTGGCTCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATg                         <  2:7745685/150‑1 (MQ=255)

                                                                                                                                   tCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCcagc                   >  1:8309419/1‑150 (MQ=255)

                                                                                                                                   tCGTAGAGCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCcagc                   >  1:714813/1‑150 (MQ=255)

                                                                                                                                         agCATCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAgcgcg             <  1:8067675/150‑1 (MQ=255)

                                                                                                                                           cATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCt           <  1:7282975/150‑1 (MQ=255)

                                                                                                                                            aTCAAGAGCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTg          <  2:4641657/150‑1 (MQ=255)

                                                                                                                                             tCAAGAGCTAGGCCGAACCGATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGc         >  1:3254474/1‑150 (MQ=255)

                                                                                                                                             tCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGc         >  1:965487/1‑150 (MQ=255)

                                                                                                                                                agagCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCaa      >  2:7314373/1‑150 (MQ=255)

                                                                                                                                                agagCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCaa      <  2:665228/150‑1 (MQ=255)

                                                                                                                                                 gagCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAag     <  2:6204058/150‑1 (MQ=255)

                                                                                                                                                  agCTAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAagc    <  2:125777/150‑1 (MQ=255)

ACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAG  >  NC_002947/5622994‑5623291

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