BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,002,601	T→G	5.4%	Q166P (CAG→CCG)	PP_4405 ←	sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,002,601	0	T	G	5.4%	1217.5 / 38.0	346	Q166P (CAG→CCG)	PP_4405	sensory box protein
Reads supporting (aligned to +/- strand): ref base T (143/183); new base G (10/9); total (153/193)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 351 total aligned reads displayed.

CCATGTCCATGGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCACGGAACCGTCC  >  NC_002947/5002452‑5002750

 catgtccatgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTg                                                                                                                                                      <  1:449641/150‑1 (MQ=255)

  atgtccatgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGc                                                                                                                                                     >  1:3203734/1‑150 (MQ=255)

    gtccatgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCaa                                                                                                                                                   <  1:6725341/150‑1 (MQ=255)

     tccatgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCAGCGGCAAc                                                                                                                                                  <  2:2231635/150‑1 (MQ=255)

       catgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTACCGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCGGTTGCCGCTGCACctg                                                                                                                                                >  2:4252768/1‑150 (MQ=255)

       catgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAActg                                                                                                                                                >  2:6509605/1‑150 (MQ=255)

         tgGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAActgct                                                                                                                                              >  2:6257263/1‑150 (MQ=255)

          ggCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAActgctg                                                                                                                                             <  1:7142021/150‑1 (MQ=255)

          ggCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAActgctg                                                                                                                                             <  2:5771596/150‑1 (MQ=255)

          ggCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAActgctg                                                                                                                                             <  2:2337536/150‑1 (MQ=255)

             aaTCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGtt                                                                                                                                          >  2:6671605/1‑150 (MQ=255)

              aTCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCAGCGGCAActggacgaga                                                                                                                                         <  1:3996703/150‑8 (MQ=255)

               tCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCa                                                                                                                                        <  2:3895782/150‑1 (MQ=255)

                 aGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCAcc                                                                                                                                      >  2:3135709/1‑150 (MQ=255)

                   cGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCa                                                                                                                                    <  2:4778200/150‑1 (MQ=255)

                   cGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCa                                                                                                                                    >  1:1361295/1‑150 (MQ=255)

                     aCAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAgg                                                                                                                                  <  2:2626476/150‑1 (MQ=255)

                      cAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGa                                                                                                                                 <  1:5564405/150‑1 (MQ=255)

                            cTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTc                                                                                                                           <  1:249888/150‑1 (MQ=255)

                             tGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCa                                                                                                                          >  1:5906056/1‑150 (MQ=255)

                              ggCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCAc                                                                                                                         <  1:3966528/150‑1 (MQ=255)

                                cGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACgg                                                                                                                       <  1:2711041/150‑1 (MQ=255)

                                cGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACgg                                                                                                                       <  2:278316/150‑1 (MQ=255)

                                     gcgcTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCg                                                                                                                  >  1:6755233/1‑150 (MQ=255)

                                      cgcTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGt                                                                                                                 >  1:5705066/1‑150 (MQ=255)

                                        cTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAActgctg                                                                                                                                             >  2:1436837/1‑120 (MQ=255)

                                        cTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTc                                                                                                               <  1:6449187/150‑1 (MQ=255)

                                        cTGGGCGCGCGGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAActgctg                                                                                                                                             <  1:1436837/120‑1 (MQ=255)

                                         tGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCg                                                                                                              >  2:862234/1‑150 (MQ=255)

                                          gggCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCgg                                                                                                             <  1:5104736/150‑1 (MQ=255)

                                           ggCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTGCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGc                                                                                                            >  2:1371815/1‑150 (MQ=255)

                                            gcgcgcTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGcc                                                                                                           <  1:1807878/150‑1 (MQ=255)

                                            gcgcgcTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGcc                                                                                                           >  1:4834747/1‑150 (MQ=255)

                                             cgcgcTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCa                                                                                                          >  2:5872744/1‑150 (MQ=255)

                                             cgcgcTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGGTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCa                                                                                                          <  1:7204826/150‑1 (MQ=255)

                                              gcgcTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCAc                                                                                                         <  1:5175018/150‑1 (MQ=255)

                                                gcTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCt                                                                                                       >  2:4166151/1‑150 (MQ=255)

                                                gcTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCt                                                                                                       <  2:3977247/150‑1 (MQ=255)

                                                    gCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGaaa                                                                                                   >  1:1004272/1‑150 (MQ=255)

                                                       cTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTc                                                                                                >  1:5398/1‑150 (MQ=255)

                                                              tGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGAc                                                                                         >  1:2594204/1‑150 (MQ=255)

                                                                cAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTg                                                                                       <  1:7378825/150‑1 (MQ=255)

                                                                    gcgcTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGttc                                                                                   >  2:678305/1‑150 (MQ=255)

                                                                    gcgcTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGttc                                                                                   >  2:7055243/1‑150 (MQ=255)

                                                                        tctcCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTtcttct                                                                               <  1:1586715/150‑1 (MQ=255)

                                                                        tctcCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTtcttct                                                                               >  1:3504348/1‑150 (MQ=255)

                                                                        tctcCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTtcttct                                                                               >  1:6320301/1‑150 (MQ=255)

                                                                        tctcCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTtcttct                                                                               <  2:505581/150‑1 (MQ=255)

                                                                        tctcCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTtcttct                                                                               >  1:5511103/1‑150 (MQ=255)

                                                                          tcCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGt                                                                             >  1:2048613/1‑150 (MQ=255)

                                                                                 ggTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCTGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGtt                                                                      <  1:3135479/150‑1 (MQ=255)

                                                                                   tAACGTTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                                    <  1:2903615/150‑1 (MQ=255)

                                                                                   tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCCCGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                                    <  1:1872188/150‑1 (MQ=255)

                                                                                   tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                                    <  2:5651598/150‑1 (MQ=255)

                                                                                   tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                                    >  1:989008/1‑150 (MQ=255)

                                                                                   tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                                    >  2:3287198/1‑150 (MQ=255)

                                                                                   tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                                    <  1:6553482/150‑1 (MQ=255)

                                                                                   tAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                                    <  2:296685/150‑1 (MQ=255)

                                                                                   caaCGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTga                                                                    <  2:4921028/149‑1 (MQ=255)

                                                                                    aaCGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTgag                                                                   <  1:1174766/150‑1 (MQ=255)

                                                                                       ggTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGacc                                                                >  2:1944283/1‑149 (MQ=255)

                                                                                       ggTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGaca                                                                >  1:7164397/1‑150 (MQ=255)

                                                                                       ggTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGaca                                                                >  2:442236/1‑150 (MQ=255)

                                                                                        gTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGacac                                                               <  1:6082638/150‑1 (MQ=255)

                                                                                        gTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGacac                                                               <  2:4602392/150‑1 (MQ=255)

                                                                                               acacGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGc                                                        >  2:3426029/1‑150 (MQ=255)

                                                                                                 acGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTTGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCtt                                                      <  2:1103773/150‑1 (MQ=255)

                                                                                                 acGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCtt                                                      <  1:7440184/150‑1 (MQ=255)

                                                                                                 acGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCtt                                                      <  2:6875602/150‑1 (MQ=255)

                                                                                                  cGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTg                                                     <  1:4318738/150‑1 (MQ=255)

                                                                                                        aaCGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCg                                               >  1:4550130/1‑150 (MQ=255)

                                                                                                        aaCGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCg                                               <  1:4068958/150‑1 (MQ=255)

                                                                                                         aCGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGt                                              <  1:4112618/150‑1 (MQ=255)

                                                                                                         aCGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGt                                              >  2:3812689/1‑150 (MQ=255)

                                                                                                                  gTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGc                                     <  2:7189210/150‑1 (MQ=255)

                                                                                                                    tCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGa                                   <  2:2444630/150‑1 (MQ=255)

                                                                                                                      ggCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATg                                 <  2:6580086/150‑1 (MQ=255)

                                                                                                                      ggCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATg                                 >  1:122080/1‑150 (MQ=255)

                                                                                                                      ggCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATg                                 >  2:4519713/1‑150 (MQ=255)

                                                                                                                      ggCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATg                                 >  1:2154049/1‑150 (MQ=255)

                                                                                                                      ggCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATg                                 >  2:2699695/1‑150 (MQ=255)

                                                                                                                         cAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATgcgc                              <  2:6099889/150‑1 (MQ=255)

                                                                                                                          cccgcccGGGTCTCGTCCAGTGGCCGCAGCAACGGCTGGTTCCCCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCa                             <  2:6002442/144‑1 (MQ=255)

                                                                                                                          aaTTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCa                             >  1:3784364/1‑150 (MQ=255)

                                                                                                                               aCGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAg                        >  1:5361005/1‑150 (MQ=255)

                                                                                                                               aCGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAg                        >  1:4696081/1‑150 (MQ=255)

                                                                                                                               aCGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAg                        >  1:4575808/1‑150 (MQ=255)

                                                                                                                               aCGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAg                        >  1:2339469/1‑150 (MQ=255)

                                                                                                                                    tctcGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGc                   <  1:3843244/150‑1 (MQ=255)

                                                                                                                                      tcGTCCAGTTGCCGCTGCAACTGCTGGTTCTCCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGt                 >  1:3576045/1‑150 (MQ=255)

                                                                                                                                      tcGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGt                 >  2:7347144/1‑150 (MQ=255)

                                                                                                                                      tcGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGt                 <  1:6509743/150‑1 (MQ=255)

                                                                                                                                        gTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGc               >  1:5997346/1‑150 (MQ=255)

                                                                                                                                        gTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGc               >  1:2246509/1‑150 (MQ=255)

                                                                                                                                         agcAGTTGCAGCGGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGcc              >  1:420254/3‑150 (MQ=255)

                                                                                                                                          gcAGTTGCAGCGGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGt                 >  1:4780982/2‑146 (MQ=255)

                                                                                                                                             gTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCACgg          >  1:6927031/1‑150 (MQ=255)

                                                                                                                                                  cGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCACGGAACCg     >  1:6692808/1‑150 (MQ=255)

                                                                                                                                                   gCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCACGGAACCGt    >  2:2016497/1‑150 (MQ=255)

                                                                                                                                                   gCTGCAACTGCAGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCAAGCGTCGATATTGCGATCCGCACCCAGCATGCGTGCCACGGAACCGt    >  2:4681186/1‑150 (MQ=255)

CCATGTCCATGGCAATCAGCGACAACGGCTGGCGGAAGCGCTGGGCGCGCTGGCACTCCAGCTGCAACGCGCTCTCCAGCCGGTAACGGTTGGCGACACGGGTCAACGCATCGCGTTCGGCCAATTCACGGTTCTCGTCCAGTTGCCGCTGCAACTGCTGGTTCACCCAGGACAACTCACGGGTGCGTTCGGCCACCTGAAACTCCAGCGACTGGTTCTTCTGTTCCAGTTGAGACACCAAACGCTTGCCAGCGTCGATATTGCGATGCGCACCCAGCATGCGTGCCACGGAACCGTCC  >  NC_002947/5002452‑5002750

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