BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,739,145	C→G	8.7%	K93N (AAG→AAC) ‡	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,739,145	0	C	G	8.7%	979.6 / 63.0	299	K93N (AAG→AAC) ‡	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base C (156/117); new base G (15/11); total (171/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 301 total aligned reads displayed.

TGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTT  >  NC_002947/5738996‑5739293

tGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCAAACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCCGATGGCCCAGTGTGGCACCAGGTTGGCGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTc                                                                                                                                                      >  2:1465469/1‑150 (MQ=255)

  gCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCtt                                                                                                                                                    <  1:7847212/150‑1 (MQ=255)

     gCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGaa                                                                                                                                                 <  1:5050390/150‑1 (MQ=255)

     gCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGaa                                                                                                                                                 <  1:4547579/150‑1 (MQ=255)

       cGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAaca                                                                                                                                               <  1:1541597/150‑1 (MQ=255)

           gTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACacc                                                                                                                                           >  1:8838064/1‑150 (MQ=255)

              gATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgac                                                                                                                                        >  2:2394436/1‑150 (MQ=255)

              gATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgac                                                                                                                                        >  2:1580631/1‑150 (MQ=255)

                tCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgaccg                                                                                                                                      >  1:510359/1‑150 (MQ=255)

                   acgaaccGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCt                                                                                                                                   >  1:4758170/6‑150 (MQ=255)

                   acgaaccGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCt                                                                                                                                   >  1:4030484/6‑150 (MQ=255)

                   aGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACAC‑CCGACCGGCtt                                                                                                                                  >  2:7842037/1‑150 (MQ=255)

                        ccGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGccc                                                                                                                              >  2:8177835/1‑150 (MQ=255)

                             ggTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTg                                                                                                                         <  1:6425485/150‑1 (MQ=255)

                                 aTCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCAcc                                                                                                                     >  2:5832542/1‑150 (MQ=255)

                                  tCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACca                                                                                                                    <  1:2016347/150‑1 (MQ=255)

                                    tGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGTCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACcagc                                                                                                                  >  1:7316822/1‑150 (MQ=255)

                                    tGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACcagc                                                                                                                  >  1:7475676/1‑150 (MQ=255)

                                      gCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACcagcag                                                                                                                >  2:660624/1‑150 (MQ=255)

                                      gCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACaagaag                                                                                                                >  2:650912/1‑150 (MQ=255)

                                       ccATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACcagcaag                                                                                                               <  1:6900669/150‑3 (MQ=255)

                                       ccATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                                                                                               <  2:8717010/150‑1 (MQ=255)

                                       ccATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAgg                                                                                                               <  2:3590198/150‑1 (MQ=255)

                                         aTACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGt                                                                                                             >  2:5750761/1‑150 (MQ=255)

                                             aTCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcg                                                                                                         >  1:3776200/1‑150 (MQ=255)

                                             aTCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcg                                                                                                         >  1:2062489/1‑150 (MQ=255)

                                               cGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcgtc                                                                                                       <  2:7078174/150‑1 (MQ=255)

                                               cGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcgtc                                                                                                       <  1:944221/150‑1 (MQ=255)

                                               cGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcgtc                                                                                                       <  1:3214978/150‑1 (MQ=255)

                                               cGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcgtc                                                                                                       <  1:8112879/150‑1 (MQ=255)

                                                  cAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCgtcgtc                                                                                                    >  1:4790562/1‑150 (MQ=255)

                                                   aGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCt                                                                                                   >  1:5432984/1‑150 (MQ=255)

                                                     gCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTc                                                                                                 >  2:5234754/1‑150 (MQ=255)

                                                       ccTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCa                                                                                               >  1:6158326/1‑150 (MQ=255)

                                                       ccTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCa                                                                                               >  1:1425997/1‑150 (MQ=255)

                                                       ccTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCa                                                                                               <  1:5084589/150‑1 (MQ=255)

                                                       ccTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCa                                                                                               >  2:1857053/1‑150 (MQ=255)

                                                              gTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGt                                                                                        <  1:5284035/150‑1 (MQ=255)

                                                              gTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGt                                                                                        >  1:3669992/1‑150 (MQ=255)

                                                               tCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTc                                                                                       >  1:3071611/1‑150 (MQ=255)

                                                                         tGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcga                                                                             >  2:3671997/1‑150 (MQ=255)

                                                                         tGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcga                                                                             >  1:8206755/1‑150 (MQ=255)

                                                                         tGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcga                                                                             <  2:4771129/150‑1 (MQ=255)

                                                                          gAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGTCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGAACGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCTGGGTTTcgac                                                                            <  2:3239642/150‑1 (MQ=255)

                                                                          gAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTcgac                                                                            >  1:4324528/1‑150 (MQ=255)

                                                                                ttCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCtt                                                                      >  1:2618232/1‑150 (MQ=255)

                                                                                ttCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCtt                                                                      >  2:4198581/1‑150 (MQ=255)

                                                                                 tCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTg                                                                     >  1:794098/1‑150 (MQ=255)

                                                                                  cccAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGt                                                                    >  2:7057220/1‑150 (MQ=255)

                                                                                   ccAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTc                                                                   >  1:669668/1‑150 (MQ=255)

                                                                                   ccAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTc                                                                   >  1:2255405/1‑150 (MQ=255)

                                                                                   ccAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTc                                                                   >  2:2457749/1‑150 (MQ=255)

                                                                                    cAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCg                                                                  >  1:2457919/1‑150 (MQ=255)

                                                                                    cAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCg                                                                  >  2:4658500/1‑150 (MQ=255)

                                                                                      gTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTa                                                                <  2:6172835/150‑1 (MQ=255)

                                                                                       ttGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCCAAGGTCATGtagtag                                                               >  2:2666040/1‑150 (MQ=255)

                                                                                          gCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAg                                                            >  2:9047960/1‑150 (MQ=255)

                                                                                          gCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAg                                                            >  1:1675040/1‑150 (MQ=255)

                                                                                           cccAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAgg                                                           <  2:6923283/150‑1 (MQ=255)

                                                                                           cccAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAgg                                                           >  1:3411579/1‑150 (MQ=255)

                                                                                           cccAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGGTTCGACGATCTGGTCGTCGCAgg                                                           >  2:8513345/1‑150 (MQ=255)

                                                                                             cAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGcc                                                         <  2:8672557/150‑1 (MQ=255)

                                                                                              aGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGccc                                                        <  2:1375345/150‑1 (MQ=255)

                                                                                               gtgtGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCa                                                       >  1:8341299/1‑150 (MQ=255)

                                                                                                    gCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTgc                                                  >  2:6378064/1‑150 (MQ=255)

                                                                                                    gCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTgc                                                  >  2:1170062/1‑150 (MQ=255)

                                                                                                    gCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTgc                                                  >  2:1778432/1‑150 (MQ=255)

                                                                                                      aCCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTgcgc                                                >  2:7368509/1‑150 (MQ=255)

                                                                                                       ccAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTgcgcg                                               <  1:7975206/150‑1 (MQ=255)

                                                                                                           gTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCgcgcggc                                           >  1:1331199/1‑150 (MQ=255)

                                                                                                           gTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCgcgcggc                                           >  2:8397183/1‑150 (MQ=255)

                                                                                                               gAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGcgcggcgcgg                                       <  1:4690068/150‑1 (MQ=255)

                                                                                                                   tGGCGATGAGAACGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGccaa                                   <  1:5216039/150‑3 (MQ=255)

                                                                                                                     gCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATa                                 >  1:423532/1‑150 (MQ=255)

                                                                                                                     gCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATa                                 >  1:8219676/1‑150 (MQ=255)

                                                                                                                     gCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATa                                 >  1:4936047/1‑150 (MQ=255)

                                                                                                                     gCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATa                                 >  2:5398305/1‑150 (MQ=255)

                                                                                                                     gCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATa                                 <  1:7057220/150‑1 (MQ=255)

                                                                                                                     gCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCCTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATa                                 >  1:9055272/1‑150 (MQ=255)

                                                                                                                        atgagaaCGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGccg                              <  2:7572815/144‑1 (MQ=255)

                                                                                                                        atgagaaCGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAACAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGccg                              >  1:5496448/7‑150 (MQ=255)

                                                                                                                                gcgcgGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATac                      <  1:1168644/150‑1 (MQ=255)

                                                                                                                                    ggTGCGTTGCTGTGGGTCGTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACAGCCGTATaccacc                  <  2:3012924/150‑1 (MQ=255)

                                                                                                                                     gTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATAccacca                 <  1:2324482/150‑1 (MQ=255)

                                                                                                                                     gTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATAccacca                 <  1:775026/150‑1 (MQ=255)

                                                                                                                                     gTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATAccacca                 <  1:770844/150‑1 (MQ=255)

                                                                                                                                       gcatccttgtgGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGc               >  2:4551120/8‑150 (MQ=255)

                                                                                                                                        catccttgtgGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGct              >  2:5157741/7‑150 (MQ=255)

                                                                                                                                          ttGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGctct            >  2:2029339/1‑150 (MQ=255)

                                                                                                                                           ccttgtgGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGctcat           <  2:2964733/147‑3 (MQ=255)

                                                                                                                                             cTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTc         >  2:414831/1‑150 (MQ=255)

                                                                                                                                             cTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTc         >  1:3007672/1‑150 (MQ=255)

                                                                                                                                             cTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTc         <  2:2233695/150‑1 (MQ=255)

                                                                                                                                              tgtgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCg        <  2:4466720/150‑1 (MQ=255)

                                                                                                                                              tgtgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCg        >  1:2999124/1‑150 (MQ=255)

                                                                                                                                               gtgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCgg       <  2:3933214/150‑1 (MQ=255)

                                                                                                                                                tgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCggg      >  1:6197956/1‑150 (MQ=255)

                                                                                                                                                tgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCggg      <  2:2330489/150‑1 (MQ=255)

                                                                                                                                                   gTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGttt   >  1:763872/1‑150 (MQ=255)

                                                                                                                                                    tCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGtttt  >  2:5634232/1‑150 (MQ=255)

TGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTT  >  NC_002947/5738996‑5739293

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