BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,738,962	G→C	28.8%	V154V (GTC→GTG)	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,738,962	0	G	C	28.8%	299.9 / 152.7	177	V154V (GTC→GTG)	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base G (56/70); new base C (27/24); total (83/94)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 186 total aligned reads displayed.

TCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTT  >  NC_002947/5738813‑5739111

  gATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTcgac                                                                                                                                                     >  1:335780/1‑150 (MQ=255)

  gATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTcgac                                                                                                                                                     >  1:4703797/1‑150 (MQ=255)

   aTACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTcgacg                                                                                                                                                    >  2:4392013/1‑150 (MQ=255)

    tACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTcgacga                                                                                                                                                   >  1:3300016/1‑150 (MQ=255)

     aCGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGaa                                                                                                                                                  <  1:3593322/150‑1 (MQ=255)

     aCGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGaa                                                                                                                                                  <  2:3728485/150‑1 (MQ=255)

      cGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAg                                                                                                                                                 <  2:938329/150‑1 (MQ=255)

      cGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAg                                                                                                                                                 <  1:2996710/150‑1 (MQ=255)

        gCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGt                                                                                                                                               >  2:5034326/1‑150 (MQ=255)

        gCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGt                                                                                                                                               >  2:5034513/1‑150 (MQ=255)

            tgctgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCa                                                                                                                                           <  1:636986/150‑1 (MQ=255)

            tgctgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCa                                                                                                                                           <  2:720866/150‑1 (MQ=255)

              ctgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATa                                                                                                                                         <  1:4074827/150‑1 (MQ=255)

               tgGCATTCGTTGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAg                                                                                                                                        <  1:2279824/150‑1 (MQ=255)

                    ttCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGcc                                                                                                                                   >  2:4521018/1‑150 (MQ=255)

                    ttCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGGTTCATAGGTGcc                                                                                                                                   >  2:2584226/1‑150 (MQ=255)

                      cGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCt                                                                                                                                 <  1:1870469/150‑1 (MQ=255)

                           ttCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACg                                                                                                                            >  2:1963768/1‑150 (MQ=255)

                            tCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACACACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGa                                                                                                                           <  2:1716899/150‑1 (MQ=255)

                             cAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGAt                                                                                                                          >  1:2040838/1‑150 (MQ=255)

                              aGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATg                                                                                                                         >  1:781624/1‑150 (MQ=255)

                              aGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATg                                                                                                                         >  1:519685/1‑150 (MQ=255)

                                gcagcaagcaCCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGcc                                                                                                                       >  1:3777254/3‑150 (MQ=255)

                                 cagcaagcaCCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGccc                                                                                                                      >  2:579571/2‑150 (MQ=255)

                                 cagcaagcaCCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGccc                                                                                                                      <  2:3023356/149‑1 (MQ=255)

                                  agcaagcaCCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCt                                                                                                                     >  2:492633/1‑150 (MQ=255)

                                  aGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGACGGTTTCATAGGTGCCCTGGACGATGCCCt                                                                                                                     >  2:82573/1‑150 (MQ=255)

                                  aGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCt                                                                                                                     >  2:2478229/1‑150 (MQ=255)

                                    cAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTgg                                                                                                                   <  1:717014/150‑1 (MQ=255)

                                    cAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTgg                                                                                                                   >  2:3992768/1‑150 (MQ=255)

                                     aaGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGc                                                                                                                  <  1:3462261/150‑1 (MQ=255)

                                     aaGCACCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTgtg                                                                                                                  >  2:604542/1‑148 (MQ=255)

                                        ccccAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGc                                                                                                               <  1:3123604/150‑1 (MQ=255)

                                        ccccAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGc                                                                                                               >  2:933880/1‑150 (MQ=255)

                                         cccAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGcc                                                                                                              <  1:93835/150‑1 (MQ=255)

                                          ccagccagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCAGAGGTGCCCTGGACGATGCCCTGGCTGCCg                                                                                                             >  1:1649348/1‑150 (MQ=255)

                                             gccagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATg                                                                                                          >  1:3455139/1‑150 (MQ=255)

                                             gccagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATg                                                                                                          >  1:1126304/1‑150 (MQ=255)

                                               cagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTa                                                                                                        <  2:1361872/150‑1 (MQ=255)

                                               cagGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTa                                                                                                        >  1:2742989/1‑150 (MQ=255)

                                                agGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAg                                                                                                       <  1:19551/150‑1 (MQ=255)

                                                   gTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATc                                                                                                    >  1:1207184/1‑150 (MQ=255)

                                                    tCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATcc                                                                                                   <  1:4589827/150‑1 (MQ=255)

                                                        ggCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCt                                                                                               >  1:4225014/1‑150 (MQ=255)

                                                         gCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTg                                                                                              >  2:2658813/1‑150 (MQ=255)

                                                           cAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGcc                                                                                            >  1:1843757/1‑150 (MQ=255)

                                                           cAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGcc                                                                                            >  1:4037380/1‑150 (MQ=255)

                                                                     gcgcCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATc                                                                                  <  1:3237069/150‑1 (MQ=255)

                                                                     gcgcCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATc                                                                                  <  1:2199744/150‑1 (MQ=255)

                                                                     gcgcCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATc                                                                                  <  1:1513381/150‑1 (MQ=255)

                                                                      ctccgccCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCt                                                                                 <  2:4062053/148‑1 (MQ=255)

                                                                       tccgccCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCtc                                                                                >  1:2476657/2‑149 (MQ=255)

                                                                       tccgccCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTg                                                                                <  2:2050041/149‑1 (MQ=255)

                                                                          gccCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGcc                                                                             >  2:2966482/1‑150 (MQ=255)

                                                                          gccCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGcc                                                                             <  1:433059/150‑1 (MQ=255)

                                                                            ccATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCAt                                                                           >  2:3172904/1‑150 (MQ=255)

                                                                             cATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATa                                                                          >  2:4228754/1‑150 (MQ=255)

                                                                                 ccgccAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACAtg                                                                      <  1:4123231/150‑2 (MQ=255)

                                                                                  cgccAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACAtgg                                                                     <  2:4579160/150‑3 (MQ=255)

                                                                                   gccCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGc                                                                    >  1:3921583/1‑150 (MQ=255)

                                                                                   gccAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGc                                                                    <  2:2105876/150‑1 (MQ=255)

                                                                                    ccaaGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGcg                                                                   >  1:3328161/4‑149 (MQ=255)

                                                                                     ccAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCa                                                                  <  1:1963768/150‑1 (MQ=255)

                                                                                       aGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAga                                                                <  1:564185/150‑2 (MQ=255)

                                                                                        ggCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAgac                                                               >  1:4353013/1‑148 (MQ=255)

                                                                                         gCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGcc                                                              >  2:1253962/1‑150 (MQ=255)

                                                                                           caccGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCt                                                            >  1:1533269/4‑150 (MQ=255)

                                                                                           caccGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCt                                                            <  1:358130/147‑1 (MQ=255)

                                                                                           cAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCt                                                            >  2:3549306/1‑150 (MQ=255)

                                                                                           cAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCt                                                            >  1:4385838/1‑150 (MQ=255)

                                                                                             gCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTg                                                          <  1:481940/150‑1 (MQ=255)

                                                                                                  gAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTc                                                     <  1:2524943/150‑1 (MQ=255)

                                                                                                          aCTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGt                                             >  2:3537798/1‑150 (MQ=255)

                                                                                                              gCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGaa                                         <  2:3312586/150‑1 (MQ=255)

                                                                                                              gCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGaa                                         <  2:1905393/150‑1 (MQ=255)

                                                                                                              gCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGaa                                         <  1:5007513/150‑1 (MQ=255)

                                                                                                              gCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGaa                                         <  1:4157945/150‑1 (MQ=255)

                                                                                                                 tttCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAgtg                                      >  1:2027558/1‑150 (MQ=255)

                                                                                                                 tttCAGGCTGCCGCCGTAGTGCTGGAGACCGGCTACGACGAAGGTTTCATAGGTGCAATGGACGATGCCCTGGCTGCCGATGAAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGGTCGTTGAAgtg                                      >  2:4900597/1‑150 (MQ=255)

                                                                                                                 cttgaGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGcga                                                                  <  2:2183028/118‑3 (MQ=255)

                                                                                                                 cttgaGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGcga                                                                  >  1:2183028/5‑120 (MQ=255)

                                                                                                                 cttgaGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAatg                                      >  1:1802923/5‑150 (MQ=255)

                                                                                                                        cTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGTCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAg                               >  2:4821373/1‑150 (MQ=255)

                                                                                                                         tGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGt                              >  2:4463674/1‑150 (MQ=255)

                                                                                                                          gccgtcATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAggt                             >  2:4139796/1‑148 (MQ=255)

                                                                                                                           ccgccgTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTg                            <  1:1403869/150‑1 (MQ=255)

                                                                                                                              ccgTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGcc                         <  2:4053989/150‑1 (MQ=255)

                                                                                                                               cgTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGccc                        <  1:2915949/150‑1 (MQ=255)

                                                                                                                                   atGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAatg                    <  1:1302783/149‑1 (MQ=255)

                                                                                                                                    tGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAatgc                   <  1:2777579/150‑2 (MQ=255)

                                                                                                                                     gCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGtgtg                  <  1:5034513/150‑1 (MQ=255)

                                                                                                                                     gCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGtgtg                  <  1:4521018/150‑1 (MQ=255)

                                                                                                                                             ccGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAga          <  1:3748581/150‑2 (MQ=255)

                                                                                                                                               ggCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGtt        <  2:121270/150‑1 (MQ=255)

                                                                                                                                               ggCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGtt        <  2:1712995/150‑1 (MQ=255)

                                                                                                                                                gCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTg       <  2:3721391/150‑1 (MQ=255)

                                                                                                                                                gCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGAGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTg       <  2:4224977/150‑1 (MQ=255)

                                                                                                                                                gCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCAGGTCGCCCAATGCGGCACGAGATTg       <  2:4579731/150‑1 (MQ=255)

                                                                                                                                                   tCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAg    <  2:4484000/150‑1 (MQ=255)

                                                                                                                                                     gacgaAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCAAGTTGGCCCAGTGTGGCACCAGGTTGGAGtt  >  1:4051489/1‑150 (MQ=255)

TCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCAGTGTGGCACCAGGTTGGAGTT  >  NC_002947/5738813‑5739111

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