BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,738,026	A→T	5.0%	A466A (GCT→GCA)	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,738,026	0	A	T	5.0%	884.4 / 22.9	258	A466A (GCT→GCA)	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base A (138/107); new base T (9/4); total (147/111)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 261 total aligned reads displayed.

TCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCAGCTTGGCGCGAAG  >  NC_002947/5737877‑5738175

tCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACa                                                                                                                                                        >  1:4575131/1‑150 (MQ=255)

   tcggcATCTTCCGAACCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCg                                                                                                                                                     >  2:2892307/2‑150 (MQ=255)

         tCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGaa                                                                                                                                               <  1:5346000/150‑1 (MQ=255)

            tCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCa                                                                                                                                            >  1:6596038/1‑150 (MQ=255)

               gaaCCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTct                                                                                                                                         >  1:8831925/3‑149 (MQ=255)

                    gTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCAt                                                                                                                                    <  2:1799913/150‑1 (MQ=255)

                     tCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATg                                                                                                                                   <  1:3588422/150‑1 (MQ=255)

                        cAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGcc                                                                                                                                >  2:5320522/1‑150 (MQ=255)

                                  ccATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgc                                                                                                                      >  2:1973632/1‑150 (MQ=255)

                                  ccATGCCAGAGCGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGAGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTGTcgc                                                                                                                      >  2:3880671/1‑150 (MQ=255)

                                   cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg                                                                                                                     >  1:8575649/1‑150 (MQ=255)

                                    aTGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCgg                                                                                                                    <  2:8817816/150‑1 (MQ=255)

                                        ccGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑g                                                                                                               <  2:2540543/150‑1 (MQ=255)

                                          gAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑ggg                                                                                                             >  1:1327361/1‑150 (MQ=255)

                                          gAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑ggg                                                                                                             <  1:4903663/150‑1 (MQ=255)

                                          gAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑ggg                                                                                                             <  2:3532745/150‑1 (MQ=255)

                                            gTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCt                                                                                                           <  1:8630200/150‑1 (MQ=255)

                                             tgctgcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTg                                                                                                          <  1:3363576/150‑1 (MQ=255)

                                             tgctgcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTg                                                                                                          <  2:4292017/150‑1 (MQ=255)

                                              gctgcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTgg                                                                                                         <  1:370538/150‑1 (MQ=255)

                                               ctgcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGa                                                                                                        <  2:6321134/150‑1 (MQ=255)

                                                tgcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGAt                                                                                                       >  1:2649747/1‑150 (MQ=255)

                                                tgcGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGAt                                                                                                       >  1:1575028/1‑150 (MQ=255)

                                                  cGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATAc                                                                                                     >  2:3340452/1‑150 (MQ=255)

                                                  cGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATAc                                                                                                     <  1:7950319/150‑1 (MQ=255)

                                                  cGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATAc                                                                                                     <  1:2039034/150‑1 (MQ=255)

                                                  cGAGAAGCCAATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATAc                                                                                                     <  1:2139568/150‑1 (MQ=255)

                                                  cGAGAAGCCAATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATAc                                                                                                     <  1:5479970/150‑1 (MQ=255)

                                                   gagaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATAcc                                                                                                    >  1:2690737/1‑150 (MQ=255)

                                                    agaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCg                                                                                                   >  1:2986664/1‑150 (MQ=255)

                                                     gaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGTGGTTC‑GGGCTGGATACCGa                                                                                                  >  2:5624239/1‑150 (MQ=255)

                                                      aaGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAg                                                                                                 <  2:1559181/150‑1 (MQ=255)

                                                         cccATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCg                                                                                              >  2:8298542/1‑150 (MQ=255)

                                                          ccATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCgg                                                                                             >  1:5279505/1‑150 (MQ=255)

                                                           cATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGa                                                                                            >  2:9046513/1‑150 (MQ=255)

                                                            aTGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAg                                                                                           <  2:2711990/150‑1 (MQ=255)

                                                            aTGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAg                                                                                           >  1:7801196/1‑150 (MQ=255)

                                                            aTGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAg                                                                                           >  2:6133012/1‑150 (MQ=255)

                                                            aTGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAg                                                                                           >  2:1168466/1‑150 (MQ=255)

                                                                ccACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGa                                                                                       >  1:6405057/1‑150 (MQ=255)

                                                                 cACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAg                                                                                      <  2:866182/150‑1 (MQ=255)

                                                                 cACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAg                                                                                      <  2:4915556/150‑1 (MQ=255)

                                                                   cGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGAGGACAGTCGGAAACATCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAggt                                                                                    <  1:5787249/150‑1 (MQ=255)

                                                                     ccACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAggtgg                                                                                  >  1:7008654/1‑150 (MQ=255)

                                                                     ccACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAggtgg                                                                                  <  1:7557249/150‑1 (MQ=255)

                                                                      cgccgccgTGATGGAGTGACACCCACGTTGCGCCCGAAGCGGTGTTGAGCAGAGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCATCCTTCATGGCTTCCGTTTCGCGGTTC‑GGCGAAGCCACCGAGCCGCTATCGAggtgat                                                                                 >  2:6880071/3‑148 (MQ=255)

                                                                           ccATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACgg                                                                            >  1:2190692/1‑150 (MQ=255)

                                                                            cATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGc                                                                           >  1:4167685/1‑150 (MQ=255)

                                                                              tGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCg                                                                         <  2:70567/150‑1 (MQ=255)

                                                                                gTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGAt                                                                       <  1:907215/150‑1 (MQ=255)

                                                                                    aGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATcacc                                                                   >  1:8218306/1‑150 (MQ=255)

                                                                                     gCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATcacca                                                                  >  1:5663353/1‑150 (MQ=255)

                                                                                            ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTc                                                                                <  2:7553899/129‑1 (MQ=255)

                                                                                            ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGc                                                           <  1:2779211/150‑1 (MQ=255)

                                                                                            ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGc                                                           >  1:3995982/1‑150 (MQ=255)

                                                                                            ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGc                                                           >  2:4116355/1‑150 (MQ=255)

                                                                                            ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGc                                                           >  1:1130957/1‑150 (MQ=255)

                                                                                            ccAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGc                                                           >  1:6193735/1‑150 (MQ=255)

                                                                                             cAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGcc                                                          <  2:6207454/150‑1 (MQ=255)

                                                                                              aGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCg                                                         >  1:1556159/1‑150 (MQ=255)

                                                                                                 tGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACa                                                      <  1:5885824/150‑1 (MQ=255)

                                                                                                     gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc                                                  >  2:4458257/1‑150 (MQ=255)

                                                                                                     gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc                                                  >  2:7504058/1‑150 (MQ=255)

                                                                                                     gccgccTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTc                                                  <  2:3969308/150‑1 (MQ=255)

                                                                                                           tGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCt                                            >  1:3775525/1‑150 (MQ=255)

                                                                                                            gCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCCTCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTcgcccggg                                           <  2:5246002/150‑5 (MQ=255)

                                                                                                              ggTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCg                                         >  2:6539748/1‑150 (MQ=255)

                                                                                                                tgtTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGa                                       <  2:4575131/150‑1 (MQ=255)

                                                                                                                  ttCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGAcc                                     >  1:7345953/1‑150 (MQ=255)

                                                                                                                   tCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCa                                    >  2:6021683/1‑150 (MQ=255)

                                                                                                                   tCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCa                                    >  1:2526400/1‑150 (MQ=255)

                                                                                                                     agcagGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCAtt                                  >  1:1928086/1‑150 (MQ=255)

                                                                                                                     agcagGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCAtt                                  >  1:4129637/1‑150 (MQ=255)

                                                                                                                      gcagGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCAttt                                 >  2:3403960/1‑150 (MQ=255)

                                                                                                                      gcagGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCAttt                                 >  2:1407020/1‑150 (MQ=255)

                                                                                                                         gggCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATcaccac                                                                 <  1:4285175/115‑1 (MQ=255)

                                                                                                                         gggCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCagggcgt                              >  2:7409455/1‑144 (MQ=255)

                                                                                                                          ggCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGtt                             <  2:3039777/150‑1 (MQ=255)

                                                                                                                          ggCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGtt                             <  2:7801196/150‑1 (MQ=255)

                                                                                                                          ggCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGtt                             <  2:8215390/150‑1 (MQ=255)

                                                                                                                           gCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTg                            >  2:8088487/1‑150 (MQ=255)

                                                                                                                           gCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTg                            <  2:6709994/150‑1 (MQ=255)

                                                                                                                            cGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGa                           >  2:3120812/1‑150 (MQ=255)

                                                                                                                               tGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCAGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGaaa                         <  2:3796765/150‑1 (MQ=255)

                                                                                                                               tGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAg                        >  1:6935057/1‑150 (MQ=255)

                                                                                                                                  gcagcGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCAGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCa                     >  2:2126740/1‑150 (MQ=255)

                                                                                                                                    agcGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAgg                   <  2:7657320/150‑1 (MQ=255)

                                                                                                                                        gCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCa               >  2:434495/1‑150 (MQ=255)

                                                                                                                                        gCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAAGctgg               <  2:3775525/150‑4 (MQ=255)

                                                                                                                                           aGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCAGCt            >  1:3090075/1‑150 (MQ=255)

                                                                                                                                             tCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCAGCTTg          >  2:6969498/1‑150 (MQ=255)

                                                                                                                                             tCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCAGCTTg          <  1:5701651/150‑1 (MQ=255)

                                                                                                                                               ggAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCAGCTTGgc        <  1:1100300/150‑1 (MQ=255)

                                                                                                                                                gAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCAGCTTGgcg       >  2:4807672/1‑150 (MQ=255)

                                                                                                                                                 aaaCAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCAGCTTGgcgc      <  2:5092987/150‑1 (MQ=255)

                                                                                                                                                 aaaCAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCAGCGGACCATTTCGTTGAAAGCCAGGCCCAGCTTGgcgc      >  2:8953810/1‑150 (MQ=255)

                                                                                                                                                  aaCAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAG‑CTGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCAGCTTGGCGCGa    >  2:2342274/1‑150 (MQ=255)

                                                                                                                                                   aCAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCAGCTTGGCGCGa    >  2:4854628/1‑150 (MQ=255)

                                                                                                                                                    cAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCAGCTTGGCGCGaa   <  2:8974880/150‑1 (MQ=255)

                                                                                                                                                     aGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCAGCTTGGCGCGAAg  >  2:5699049/1‑150 (MQ=255)

TCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTC‑GGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACGGCCGATCACCACCGGTGCCGACAGCTCGCCGCTGCGGACCATTTCGTTGAAAGCCAGGCCCAGCTTGGCGCGAAG  >  NC_002947/5737877‑5738175

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