BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,064,969	C→G	6.1%	E446Q (GAG→CAG) ‡	PP_3581 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,064,969	0	C	G	6.1%	603.6 / 21.4	181	E446Q (GAG→CAG) ‡	PP_3581	membrane protein
Reads supporting (aligned to +/- strand): ref base C (82/88); new base G (5/6); total (87/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 187 total aligned reads displayed.

CCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTGGTCCT  >  NC_002947/4064821‑4065118

ccGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTcc                                                                                                                                                      <  1:4916237/150‑1 (MQ=255)

   gTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAgc                                                                                                                                                   <  2:1081342/150‑1 (MQ=255)

       tttCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCa                                                                                                                                               <  2:548951/150‑1 (MQ=255)

         tCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCAgt                                                                                                                                             <  1:1641922/150‑2 (MQ=255)

         tCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCAgt                                                                                                                                             <  1:1641870/150‑2 (MQ=255)

         tCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAgg                                                                                                                                             <  2:345312/150‑1 (MQ=255)

          ccGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCAgtt                                                                                                                                            <  2:3529837/150‑3 (MQ=255)

          ccGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGgc                                                                                                                                            >  2:615710/1‑150 (MQ=255)

              cAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCaa                                                                                                                                        >  2:4952538/1‑150 (MQ=255)

              cAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCaa                                                                                                                                        <  1:2373580/150‑1 (MQ=255)

              cAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCaa                                                                                                                                        <  1:2221301/150‑1 (MQ=255)

                   cGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCt                                                                                                                                   >  1:3294571/1‑150 (MQ=255)

                    gATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTg                                                                                                                                  <  2:3213622/150‑1 (MQ=255)

                     aTGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgc                                                                                                                                 <  2:2241083/150‑1 (MQ=255)

                          gCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACCTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGcagcag                                                                                                                            <  2:4667151/150‑1 (MQ=255)

                           cGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCagcagc                                                                                                                           >  2:773910/1‑150 (MQ=255)

                             gTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCa                                                                                                                                               <  1:2913903/128‑1 (MQ=255)

                             gTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCa                                                                                                                                               >  2:2913903/1‑128 (MQ=255)

                                  gcagcagcaGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAgg                                                                                                                                             >  1:2648152/1‑125 (MQ=255)

                                  gcagcagcaGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAgg                                                                                                                                             <  2:2648152/125‑1 (MQ=255)

                                  gcagcagcaGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAAct                                                                                                                    >  1:3426787/1‑150 (MQ=255)

                                  gcagcagcaGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAAct                                                                                                                    >  2:1163067/1‑150 (MQ=255)

                                    agcagcagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAActgc                                                                                                                  <  2:319374/150‑1 (MQ=255)

                                    agcagcagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAActgc                                                                                                                  <  1:5189824/150‑1 (MQ=255)

                                         cagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGcc                                                                                                                                                <  2:1778744/115‑1 (MQ=255)

                                         cagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGcc                                                                                                                                                >  1:1778744/1‑115 (MQ=255)

                                          agGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGgcgc                                                                                                            <  2:2086365/150‑1 (MQ=255)

                                             cccTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTTGCGCGCa                                                                                                         <  1:1008008/150‑1 (MQ=255)

                                             cccTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCa                                                                                                         <  1:1008009/150‑1 (MQ=255)

                                                tgctgcGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTg                                                                                                      >  1:4084197/1‑150 (MQ=255)

                                                tgctgcGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTg                                                                                                      >  2:5161564/1‑150 (MQ=255)

                                                  ctgcGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCt                                                                                                    <  2:798892/150‑1 (MQ=255)

                                                        gTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCAt                                                                                              >  1:2147770/1‑150 (MQ=255)

                                                         tGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATc                                                                                             <  2:952851/150‑1 (MQ=255)

                                                          gCTCCCAGCGCAGCAACGCTTCGGCACCGTTCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCg                                                                                            >  1:449935/1‑150 (MQ=255)

                                                          gCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCg                                                                                            >  2:1571666/1‑150 (MQ=255)

                                                           cTCCCAGCGCAGCAACGCTTCGGCAACGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGa                                                                                           >  1:4750028/1‑150 (MQ=255)

                                                             cccAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGaca                                                                                         >  2:2727436/1‑150 (MQ=255)

                                                              ccAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGacac                                                                                        >  1:1602790/1‑150 (MQ=255)

                                                              ccAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGacac                                                                                        >  2:2860244/1‑150 (MQ=255)

                                                                     agcaACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCgaag                                                                                 >  1:3738560/1‑150 (MQ=255)

                                                                       caACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCgaagaa                                                                               >  1:1039454/1‑150 (MQ=255)

                                                                        aaCGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACATATCGaagaag                                                                              >  2:3699320/1‑150 (MQ=255)

                                                                        aaCGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGaagaag                                                                              >  1:769824/1‑150 (MQ=255)

                                                                        aaCGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGaagaag                                                                              >  1:4657679/1‑150 (MQ=255)

                                                                          cGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCt                                                                            <  1:840662/150‑1 (MQ=255)

                                                                          cGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCt                                                                            >  2:1362202/1‑150 (MQ=255)

                                                                           gCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTg                                                                           >  1:897210/1‑150 (MQ=255)

                                                                           gCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTg                                                                           >  1:3439545/1‑150 (MQ=255)

                                                                            cTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGaagaag                                                                              >  1:145649/1‑146 (MQ=255)

                                                                            cTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGaagaag                                                                              <  2:145649/146‑1 (MQ=255)

                                                                             ttCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTa                                                                         <  2:728272/150‑1 (MQ=255)

                                                                             ttCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTa                                                                         >  2:4843020/1‑150 (MQ=255)

                                                                               cGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAAc                                                                       <  2:585359/150‑1 (MQ=255)

                                                                                 gCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCg                                                                     >  1:689798/1‑150 (MQ=255)

                                                                                   aCCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGtt                                                                   <  1:4260802/150‑1 (MQ=255)

                                                                                    ccGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTc                                                                  >  1:309866/1‑150 (MQ=255)

                                                                                    ccGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTc                                                                  <  2:1632286/150‑1 (MQ=255)

                                                                                    ccGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTc                                                                  <  1:4952538/150‑1 (MQ=255)

                                                                                        tCGGTTGGCAGGCGTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAAGTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGc                                                              <  2:5076846/150‑1 (MQ=255)

                                                                                           gTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGt                                                                          <  2:4050899/135‑1 (MQ=255)

                                                                                           gTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCgg                                                           >  1:3896207/1‑150 (MQ=255)

                                                                                             tGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGgcg                                                         <  2:2162670/150‑1 (MQ=255)

                                                                                             tGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGgcg                                                         <  1:3662398/150‑1 (MQ=255)

                                                                                               gCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTc                                                                                                                                     <  1:4117803/72‑1 (MQ=255)

                                                                                               gCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCagcagc                                                                                                                           >  1:5137229/1‑82 (MQ=255)

                                                                                               gCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGAGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTc                                                                                                                                     >  2:4117803/1‑72 (MQ=255)

                                                                                                 aGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCttt                                                     <  2:4365636/150‑1 (MQ=255)

                                                                                                 aGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCttt                                                     <  1:4020306/150‑1 (MQ=255)

                                                                                                 aGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCttt                                                     >  2:3087243/1‑150 (MQ=255)

                                                                                                  ggCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCtttt                                                    <  1:4366550/150‑1 (MQ=255)

                                                                                                   gtcTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTg                                                   <  1:4511102/148‑1 (MQ=255)

                                                                                                   gCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGa                                                                                           <  1:4787185/110‑1 (MQ=255)

                                                                                                   gCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGa                                                                                           >  2:4787185/1‑110 (MQ=255)

                                                                                                       gggCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCAt                                               <  1:5029557/150‑1 (MQ=255)

                                                                                                        ggCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATg                                              <  1:1364630/150‑1 (MQ=255)

                                                                                                         gCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACATATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATgg                                             <  1:3699320/150‑1 (MQ=255)

                                                                                                            tCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTAc                                          >  2:5357605/1‑150 (MQ=255)

                                                                                                              gAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACAt                                        >  1:619164/1‑150 (MQ=255)

                                                                                                              gAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACAt                                        >  2:4815990/1‑150 (MQ=255)

                                                                                                                    gggCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGc                                  <  1:626405/150‑1 (MQ=255)

                                                                                                                      gCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGt                                >  2:1269379/1‑150 (MQ=255)

                                                                                                                       cTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTc                               <  2:1512540/150‑1 (MQ=255)

                                                                                                                           tAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCg                           >  2:5371187/1‑150 (MQ=255)

                                                                                                                               gCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTAc                       >  2:4069743/1‑150 (MQ=255)

                                                                                                                                cAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGc                            <  1:5102722/144‑1 (MQ=255)

                                                                                                                                cAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACg                      >  2:3895128/1‑150 (MQ=255)

                                                                                                                                   ttgcgccTGGCGCTGGAGCAGCGCCAGTCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACgcag                   >  1:2761363/7‑150 (MQ=255)

                                                                                                                                       gccTGGCGCTGGAGCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACAt                                        <  1:4033949/123‑1 (MQ=255)

                                                                                                                                         cTGGCGCTGGAGCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTc             >  1:3921616/1‑150 (MQ=255)

                                                                                                                                         cTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGTCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACACAGCAGCTc             <  2:3791828/150‑1 (MQ=255)

                                                                                                                                         cTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTc             <  2:888604/150‑1 (MQ=255)

                                                                                                                                         cTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTc             >  1:888388/1‑150 (MQ=255)

                                                                                                                                         cTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTc             <  2:1903301/150‑1 (MQ=255)

                                                                                                                                           ggCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGt           <  1:1348841/150‑1 (MQ=255)

                                                                                                                                            gcgcTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCgg                             <  2:5109024/131‑1 (MQ=255)

                                                                                                                                            gcgcTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTg          >  2:3698891/1‑150 (MQ=255)

                                                                                                                                             cgcTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGc         <  2:4023656/150‑1 (MQ=255)

                                                                                                                                             cgcTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGc         <  2:2078047/150‑1 (MQ=255)

                                                                                                                                                tgctCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTgg      >  2:2094720/1‑150 (MQ=255)

CCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTGGTCCT  >  NC_002947/4064821‑4065118

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