BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,737,951	G→A	5.2%	G491G (GGC→GGT)	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,737,951	0	G	A	5.2%	991.8 / 25.7	308	G491G (GGC→GGT)	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand): ref base G (116/176); new base A (7/9); total (123/185)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 311 total aligned reads displayed.

GATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACG  >  NC_002947/5737804‑5738100

 aTGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGTTGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCTTGCGCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCa                                                                                                                                                    <  1:4948766/150‑1 (MQ=255)

   gTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCTAGAAGCCCATGCCCACGCCACCGCCATg                                                                                                                                                  <  2:8683406/150‑1 (MQ=255)

    tCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATgg                                                                                                                                                 <  1:7906470/150‑1 (MQ=255)

    tCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATgg                                                                                                                                                 <  2:4876178/150‑1 (MQ=255)

    tCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTAGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATgg                                                                                                                                                 >  2:2039034/1‑150 (MQ=255)

     cATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGt                                                                                                                                                <  1:741824/150‑1 (MQ=255)

     cATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCAGCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGt                                                                                                                                                >  1:8282390/1‑150 (MQ=255)

        aaCCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCa                                                                                                                                             >  2:5885824/1‑150 (MQ=255)

        aaCCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCa                                                                                                                                             <  1:6632408/150‑1 (MQ=255)

          ccGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGc                                                                                                                                           <  2:1901998/150‑1 (MQ=255)

          ccGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGc                                                                                                                                           <  1:2894225/150‑1 (MQ=255)

           cGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCg                                                                                                                                          >  2:7447832/1‑150 (MQ=255)

            ggCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGa                                                                                                                                         <  1:578694/150‑1 (MQ=255)

                 cGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATAccc                                                                                                                                    <  1:2390169/150‑1 (MQ=255)

                          gCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGgcgt                                                                                                                           <  2:4470983/150‑2 (MQ=255)

                          gCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGgcgc                                                                                                                           >  1:6243929/1‑150 (MQ=255)

                          gCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGgcgc                                                                                                                           >  1:1580495/1‑150 (MQ=255)

                          gCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGgcgc                                                                                                                           >  1:6720509/1‑150 (MQ=255)

                           cATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGcgcc                                                                                                                          >  2:370538/1‑150 (MQ=255)

                           cATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATG‑CCACGCCACCGCCATGGTGCAg                                                                                                                                            <  1:8939663/131‑1 (MQ=255)

                           cATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATG‑CCACGCCACCGCCATGGTGCAg                                                                                                                                            >  2:8939663/1‑131 (MQ=255)

                            aTGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGcgccg                                                                                                                         >  1:7600190/1‑150 (MQ=255)

                              gACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCgccgcc                                                                                                                       >  2:517174/1‑150 (MQ=255)

                                  ccAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCg                                                                                                                   >  2:2866790/1‑150 (MQ=255)

                                       cccTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGtt                                                                                                              <  1:300807/150‑1 (MQ=255)

                                         cTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcc                                                                                                            >  2:3422476/1‑149 (MQ=255)

                                          tGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcag                                                                                                           >  1:5092987/1‑150 (MQ=255)

                                           gggTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagc                                                                                                          <  2:7067729/150‑1 (MQ=255)

                                           gggTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagc                                                                                                          <  1:4837222/150‑1 (MQ=255)

                                           gggTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagc                                                                                                          >  2:5756310/1‑150 (MQ=255)

                                            ggTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCATGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagca                                                                                                         >  1:8090504/1‑150 (MQ=255)

                                            ggTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagca                                                                                                         <  2:7714497/150‑1 (MQ=255)

                                               cGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAggg                                                                                                      <  1:2257951/150‑1 (MQ=255)

                                                gTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGc                                                                                                     <  2:6766879/150‑1 (MQ=255)

                                                gTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGc                                                                                                     <  2:3894310/150‑1 (MQ=255)

                                                   ggTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGtt                                                                                                  <  1:5985719/150‑1 (MQ=255)

                                                   ggTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGtt                                                                                                  <  2:3276678/150‑1 (MQ=255)

                                                    gTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTg                                                                                                 <  1:2564506/150‑1 (MQ=255)

                                                     tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCGTGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                                                <  2:6134304/150‑1 (MQ=255)

                                                     tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                                                <  2:1323531/150‑1 (MQ=255)

                                                     tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                                                <  2:1689175/150‑1 (MQ=255)

                                                        gTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagca                                                                                             <  1:7783212/150‑1 (MQ=255)

                                                          aCCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagcagc                                                                                           <  1:6121518/150‑1 (MQ=255)

                                                          aCCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagcagc                                                                                           <  1:7040913/150‑1 (MQ=255)

                                                           cccGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCCGACGATGACCATGCCCGAGCGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCGTGCGGTGTTCAGCAGGGCGGTGAGCcgcg                                                                                          >  1:5246002/1‑150 (MQ=255)

                                                           cccGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCg                                                                                          <  1:6399343/150‑1 (MQ=255)

                                                           cccGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCg                                                                                          <  2:6957662/150‑1 (MQ=255)

                                                            ccGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCgg                                                                                         <  2:267059/150‑1 (MQ=255)

                                                            ccGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCgg                                                                                         <  2:3972070/150‑1 (MQ=255)

                                                             cGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGc                                                                                        >  2:5986859/1‑150 (MQ=255)

                                                               ggCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCa                                                                                      <  2:4783124/150‑1 (MQ=255)

                                                                 cGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGt                                                                                    <  1:3436192/150‑1 (MQ=255)

                                                                  gATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTc                                                                                   <  2:4719174/150‑1 (MQ=255)

                                                                  gATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTc                                                                                   <  1:3879961/150‑1 (MQ=255)

                                                                     gcgcTGGTCGGCATCTTCCGAACCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGa                                                                                <  1:4150559/150‑1 (MQ=255)

                                                                     gcgcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGa                                                                                <  2:1434003/150‑1 (MQ=255)

                                                                      cgcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaa                                                                               <  1:3603490/150‑1 (MQ=255)

                                                                       gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                                              >  1:3796765/1‑150 (MQ=255)

                                                                       gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                                              <  2:2424623/150‑1 (MQ=255)

                                                                       gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                                              <  2:3805655/150‑1 (MQ=255)

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

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

                                                                              ggcCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTc                                                                       <  2:5685931/150‑1 (MQ=255)

                                                                              ggcCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGCTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCggc                                                                       >  2:3824794/1‑148 (MQ=255)

                                                                                 cTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGa                                                                    >  1:2774908/1‑150 (MQ=255)

                                                                                    ttccgaaCCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGcc                                                                 >  1:1125367/7‑150 (MQ=255)

                                                                                    aTCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGGGAAGCCCATGCCCACGCCACCGCCATGGTGCAg                                                                                                                                            <  2:3110702/75‑1 (MQ=255)

                                                                                       ggTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATc                                                              <  1:6592507/150‑1 (MQ=255)

                                                                                        gaaCCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTct                                                             >  2:9014121/3‑149 (MQ=255)

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

                                                                                             gTCGCAGTCGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCTCCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCTGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCAt                                                        <  1:3170976/150‑1 (MQ=255)

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

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

                                                                                                 cAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCGCCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGcc                                                    >  2:3531098/1‑150 (MQ=255)

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

                                                                                                   gACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTc                                                  <  2:1113309/150‑1 (MQ=255)

                                                                                                       aTGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGtt                                              <  1:5414446/150‑1 (MQ=255)

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

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

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

                                                                                                              tGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGt                                       >  1:5677200/1‑150 (MQ=255)

                                                                                                                 ccGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCg                                    >  2:547738/1‑150 (MQ=255)

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

                                                                                                                 ccGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTCTTTCATGGATTCGGTTTCACGGtttg                                    <  1:4253139/150‑3 (MQ=255)

                                                                                                                   gAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCggg                                  <  1:2107994/150‑1 (MQ=255)

                                                                                                                     gTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCt                                <  2:8282390/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                                            gagaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATAcc                         >  1:1562617/1‑150 (MQ=255)

                                                                                                                              gaAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGa                       >  2:1599766/1‑150 (MQ=255)

                                                                                                                                  cccATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCg                   >  1:4034197/1‑150 (MQ=255)

                                                                                                                                    cATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGa                 >  1:6835043/1‑150 (MQ=255)

                                                                                                                                          cACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAg           >  2:8378699/1‑150 (MQ=255)

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

                                                                                                                                          cACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAg           <  1:7986159/150‑1 (MQ=255)

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

GATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCTGGATACCGAGCCGGAGTCGAGGTGGTCACG  >  NC_002947/5737804‑5738100

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