BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,738,941	G→A	14.8%	Y161Y (TAC→TAT)	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,738,941	0	G	A	14.8%	821.5 / 111.2	325	Y161Y (TAC→TAT)	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base G (116/161); new base A (13/35); total (129/196)
Fisher's exact test for biased strand distribution p-value = 5.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 335 total aligned reads displayed.

TAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCA  >  NC_002947/5738792‑5739090

tAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGccgccg                                                                                                                                                       >  1:7770814/1‑150 (MQ=255)

  gCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTa                                                                                                                                                     >  1:4241286/1‑150 (MQ=255)

  gCGGGGGTCCAGGCGGCCGTCGATCCGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTa                                                                                                                                                     <  1:4581667/150‑1 (MQ=255)

   cGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAg                                                                                                                                                    >  2:7427484/1‑150 (MQ=255)

    gggTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGt                                                                                                                                                   >  2:3851752/1‑150 (MQ=255)

    gggTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGt                                                                                                                                                   >  1:6823034/1‑150 (MQ=255)

     ggTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTg                                                                                                                                                  <  1:1873557/150‑1 (MQ=255)

     ggTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTg                                                                                                                                                  >  1:4853409/1‑150 (MQ=255)

            aGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGa                                                                                                                                           >  1:7981863/1‑150 (MQ=255)

             ggCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGAc                                                                                                                                          <  2:8633956/150‑1 (MQ=255)

                ggAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCgg                                                                                                                                       >  1:1333090/1‑150 (MQ=255)

                   aGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCtt                                                                                                                                    >  1:5464597/1‑150 (MQ=255)

                    gTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTc                                                                                                                                   <  1:2103723/150‑1 (MQ=255)

                     tCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTcg                                                                                                                                  <  2:6541413/150‑1 (MQ=255)

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

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

                            ggCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAgg                                                                                                                           >  2:2143544/1‑150 (MQ=255)

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

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

                                   ctgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATa                                                                                                                    >  2:2963299/1‑150 (MQ=255)

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

                                   cgggCATTCGCGGTTCAGCGAGCACGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATa                                                                                                                    <  1:8738436/148‑1 (MQ=255)

                                    tgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAg                                                                                                                   <  1:414682/150‑1 (MQ=255)

                                    tgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAg                                                                                                                   >  1:9147825/1‑150 (MQ=255)

                                    tgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAg                                                                                                                   <  1:9034990/150‑1 (MQ=255)

                                    tgGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAg                                                                                                                   >  2:5911479/1‑150 (MQ=255)

                                      gCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGt                                                                                                                 <  1:6548512/150‑1 (MQ=255)

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

                                         ttCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGcc                                                                                                              <  1:51821/150‑1 (MQ=255)

                                           cgcgggtCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCGGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCt                                                                                                            <  1:7714886/144‑1 (MQ=255)

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

                                              tGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGa                                                                                                         <  2:7349150/150‑1 (MQ=255)

                                              tGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGa                                                                                                         <  1:1607652/150‑1 (MQ=255)

                                              tGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGa                                                                                                         <  2:7981863/150‑1 (MQ=255)

                                               gTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGAc                                                                                                        <  1:3156766/150‑1 (MQ=255)

                                                 tCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGa                                                                                                      <  1:8393492/150‑1 (MQ=255)

                                                  cAGCGAGCAAGCCCCTGCCTGGGTCGCGGCCTTTGGCTGGGCGCCGCCCATGCCGCCCAGGCCTGCGGTTTGTACCCACTTGCATTTCTGGGTGCCGCCGTTGTGCTGGCGACCGGCTTCGACTAAGTTTTCATTGGTGCCCTGGACTAt                                                                                                     <  2:4035223/150‑1 (MQ=255)

                                                      gAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGccc                                                                                                 <  2:3283382/150‑1 (MQ=255)

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

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

                                                         cAAGCACCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTgt                                                                                              <  2:4755473/150‑2 (MQ=255)

                                                         cAAGCACCAGCCATGACGGCTGCGAGAGGCTGTGCTCCGCCCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTgt                                                                                              <  2:5660658/150‑2 (MQ=255)

                                                              cccAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGcc                                                                                         >  1:4521663/1‑150 (MQ=255)

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

                                                                      gggTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGa                                                                                 <  1:3673032/150‑1 (MQ=255)

                                                                      gggTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGa                                                                                 >  2:7847212/1‑150 (MQ=255)

                                                                           gcgGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAg                                                                            >  2:6425485/1‑150 (MQ=255)

                                                                               ccAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGc                                                                        <  1:8644971/150‑1 (MQ=255)

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

                                                                                     gCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCcggcgg                                                                  <  1:2837000/150‑1 (MQ=255)

                                                                                     gCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCcggcgg                                                                  <  1:5626820/150‑1 (MQ=255)

                                                                                       tGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGt                                                                                    <  1:5904598/130‑1 (MQ=255)

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

                                                                                           cgccgcCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCt                                                            <  2:3455365/150‑1 (MQ=255)

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

                                                                                             ccgccCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCATTTCAGGCTGACGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTgg                                                          <  2:1888376/150‑1 (MQ=255)

                                                                                              cgccCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGc                                                         <  1:7082314/150‑1 (MQ=255)

                                                                                              cgccCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGc                                                         >  2:6969380/1‑150 (MQ=255)

                                                                                              cgccCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGc                                                         <  2:6132506/150‑1 (MQ=255)

                                                                                              cgccCATGCCGCCAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGc                                                         <  1:8714925/150‑1 (MQ=255)

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

                                                                                                 ccATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCAt                                                      <  1:3851752/150‑1 (MQ=255)

                                                                                                  cATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATa                                                     <  1:8083458/150‑1 (MQ=255)

                                                                                                     gccgccAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACAt                                                  <  2:4977158/150‑1 (MQ=255)

                                                                                                     gccgccAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACAt                                                  <  2:2701140/150‑1 (MQ=255)

                                                                                                       cgccAAGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACAtgg                                                >  2:4616760/1‑148 (MQ=255)

                                                                                                        gccCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTc                                                                >  2:1468995/1‑133 (MQ=255)

                                                                                                          caaGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGcga                                             <  1:6415481/148‑3 (MQ=255)

                                                                                                          caaGGCCACCGGTCAGAATCCATTTGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGcga                                             <  1:339731/148‑3 (MQ=255)

                                                                                                           cAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAg                                            <  1:6797181/150‑1 (MQ=255)

                                                                                                           cAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAg                                            <  2:7770814/150‑1 (MQ=255)

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

                                                                                                                cAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCt                                       <  1:1621628/150‑1 (MQ=255)

                                                                                                                  gCGGTGAGTACCCCCTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGTTCCAGCTGCCGGCGGTCATCTGGCCATCCCTCGCCAGGCCCTTg                                     <  2:2796005/150‑1 (MQ=255)

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

                                                                                                                         gTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCa                              <  1:6342318/150‑1 (MQ=255)

                                                                                                                            cccACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGtt                           <  1:7565275/150‑1 (MQ=255)

                                                                                                                            cccACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGtt                           <  2:828620/150‑1 (MQ=255)

                                                                                                                             ccACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTc                          >  1:5811779/1‑150 (MQ=255)

                                                                                                                             ccACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTc                          >  1:6172835/1‑150 (MQ=255)

                                                                                                                                tttGCCCTTGAGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGGCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGtt                       <  2:1758650/149‑1 (MQ=255)

                                                                                                                                cTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGtt                       >  2:6795487/1‑150 (MQ=255)

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

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

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

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

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

                                                                                                                                   gCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGaa                    >  1:4002670/1‑150 (MQ=255)

                                                                                                                                    cccttgaGGTTGCCGTCATAATGCTGGCGTCCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGaaa                   <  2:4262423/144‑2 (MQ=255)

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

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

                                                                                                                                      ctagaGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTGCTTATCCAGCTCGTTGAAatg                 >  2:5107956/5‑150 (MQ=255)

                                                                                                                                       ttgaGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAatga                <  2:3581004/147‑2 (MQ=255)

                                                                                                                                       ttgaGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAatga                >  1:3494321/4‑149 (MQ=255)

                                                                                                                                       ttgaGGTTGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAatga                >  2:8135004/4‑149 (MQ=255)

                                                                                                                                         cAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTc              <  1:8458294/150‑1 (MQ=255)

                                                                                                                                            gttGCCGTCATAATGCTGGCGACCGGCTTCCACGAAGGTTTCATAGGTGCCCTGAACAATGCCCTGTGCACCGATGTAAATCCACGAACCGGCGGTCATCTGGCCATACATGGCGAGACCCTTCTTATCCAGCTCGTTGAAATGATCCCa           >  2:4455515/3‑150 (MQ=255)

                                                                                                                                            gCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCa           >  1:7820081/1‑150 (MQ=255)

                                                                                                                                            gCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCa           <  2:1321047/150‑1 (MQ=255)

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

TAGCGGGTTTCCAGGCGGAAGTCGATACGGCTCTGCTGGCATTCGATGTTCAGCGAGCAAGCCCCAGCCAGGGTCGCGGCCAGTGGCTGGGCGCCGCCCATGCCGCCCAGGCCAGCGGTGAGTACCCACTTGCCTTTCAGGCTGCCGCCGTAGTGCTGGCGACCGGCTTCGACGAAGGTTTCATAGGTGCCCTGGACGATGCCCTGGCTGCCGATGTAGATCCAGCTGCCGGCGGTCATCTGGCCATACATCGCCAGGCCCTTGGCGTCCAGTTCGTTGAAGTGTTCCCAGTTGGCCCA  >  NC_002947/5738792‑5739090

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