BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,739,154	A→G	8.8%	G90G (GGT→GGC)	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,739,154	0	A	G	8.8%	881.0 / 50.4	294	G90G (GGT→GGC)	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base A (159/108); new base G (15/11); total (175/119)
Fisher's exact test for biased strand distribution p-value = 8.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 327 total aligned reads displayed.

ATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTG  >  NC_002947/5739005‑5739303

     gATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgac                                                                                                                                                  <  2:2970573/150‑1 (MQ=255)

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

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

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

      aTCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACaccgacc                                                                                                                                                 <  2:417419/150‑1 (MQ=255)

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

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

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

            cTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTg                                                                                                                                           <  2:6934009/150‑1 (MQ=255)

                 ggcggTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGa                                                                                                                                      <  2:1782366/150‑1 (MQ=255)

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

                             gCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACcagcaa                                                                                                                          >  2:1344117/1‑149 (MQ=255)

                             gCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTGGAATTGGCGATGAGAACGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACcagcaa                                                                                                                          >  2:6037348/1‑149 (MQ=255)

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

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

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

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

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

                                     tCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTtcgt                                                                                                                  >  2:5695434/1‑150 (MQ=255)

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

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

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

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

                                        ccAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTcgtcgt                                                                                                               >  2:3638758/1‑150 (MQ=255)

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

                                            gCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTc                                                                                                           <  2:6562099/150‑1 (MQ=255)

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

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

                                              ccTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGCTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCa                                                                                                         >  2:1464524/1‑150 (MQ=255)

                                               cTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAg                                                                                                        <  1:7389546/150‑1 (MQ=255)

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

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

                                                         aGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAgg                                                                                              <  2:5297187/150‑1 (MQ=255)

                                                          gTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAggg                                                                                             >  1:2031638/1‑150 (MQ=255)

                                                             cGTTGTAGTGGTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGttt                                                                                          <  2:7820081/150‑1 (MQ=255)

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

                                                                    gtgtTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGAt                                                                                   >  2:6427351/1‑150 (MQ=255)

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

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

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

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

                                                                          ccAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTc                                                                             <  2:4521663/150‑1 (MQ=255)

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

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

                                                                          ccAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTc                                                                             <  1:3027938/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                   ccAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGc                                                                    >  2:7386896/1‑150 (MQ=255)

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

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

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

                                                                                        gtgGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGt                                                               >  1:310328/1‑150 (MQ=255)

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

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

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

                                                                                                   ttGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCgcgcggcg                                                    >  1:9008320/1‑150 (MQ=255)

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

                                                                                                         ttGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCg                                              >  2:4618939/1‑150 (MQ=255)

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

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

                                                                                                             cAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATAc                                          <  1:1464524/150‑1 (MQ=255)

                                                                                                              gatgagaaCGCGCGGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGcc                                         >  2:7133865/8‑150 (MQ=255)

                                                                                                               aTCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATAccg                                        <  2:5447980/150‑1 (MQ=255)

                                                                                                                tCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATAccgc                                       >  2:468250/1‑150 (MQ=255)

                                                                                                                   gCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGt                                    >  2:4052976/1‑150 (MQ=255)

                                                                                                                    cACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGta                                   >  1:4968114/1‑150 (MQ=255)

                                                                                                                     aCGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGtat                                  >  2:3386278/1‑150 (MQ=255)

                                                                                                                      cgcgcgGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGtata                                 <  1:5478275/150‑1 (MQ=255)

                                                                                                                      cgcgcgGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCAtaga                                 <  1:7248109/150‑3 (MQ=255)

                                                                                                                      cgcgcgGTGCATCCTTGTGGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCAtaga                                 <  1:7716701/150‑3 (MQ=255)

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

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

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

                                                                                                                           ggTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATaccacc                            <  2:3013242/150‑1 (MQ=255)

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

                                                                                                                              gcatccttgtgGGTGCGGAACACGCCCACTGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGc                         <  2:8200825/143‑1 (MQ=255)

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

                                                                                                                                 tccttgtgGGTGCGGAACACGCCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGctca                      >  1:6196358/5‑149 (MQ=255)

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

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

                                                                                                                                      gtgGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCgg                 >  1:3374370/1‑150 (MQ=255)

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

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

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

                                                                                                                                             ttGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTcc          >  1:4563727/1‑150 (MQ=255)

                                                                                                                                             ttGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTcc          >  1:4931121/1‑150 (MQ=255)

                                                                                                                                             ttGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTcc          >  1:4361450/1‑150 (MQ=255)

                                                                                                                                             ttGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTcc          >  1:6111682/1‑150 (MQ=255)

                                                                                                                                               gAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCcgc        >  1:1592148/1‑150 (MQ=255)

                                                                                                                                               gAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCcgc        >  2:6304286/1‑150 (MQ=255)

                                                                                                                                                 acacgcCCACCGGCTTGCCCGACTGCACCAGCAAGGTTTCGTCTTCATTGAGCGTCTTGAGGGTCGCGACAATCTTGTCGAAATCATCCCATGTGCGCGCAGCGCGGCCAATGCCGCCATAGACGACCAGCTCATGCGGGCGTTCGGCGa      <  1:430571/150‑1 (MQ=255)

                                                                                                                                                   acacCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGAcc    <  1:1288751/150‑1 (MQ=255)

                                                                                                                                                     accgaccgGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCtg  >  2:6720666/1‑150 (MQ=255)

                                                                                                                                                     accgaccgGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCtg  >  1:4480140/1‑150 (MQ=255)

                                                                                                                                                     accgaccgGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCtg  >  1:2782993/1‑150 (MQ=255)

                                                                                                                                                     accgaccgGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCtg  >  2:4286111/1‑150 (MQ=255)

ATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTTGGCAATCAGCACGCGCGGTGCGTTGCTGTGGGTCTTGAACACACCGACCGGCTTGCCCGACTGCACCAGCAGGGTTTCGTCGTCTTCCAGGCGGGTCAGGGTTTCGACGATCTTGTCGTAGCAGGCCCAGTTGCGCGCGGCGCGGCCGATACCGCCGTATACCACCAGCTCTTTCGGGTTTTCCGCGACCTG  >  NC_002947/5739005‑5739303

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