BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,612,407	A→T	5.5%	K29* (AAG→TAG) ‡	PP_4079 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,612,407	0	A	T	5.5%	966.9 / 29.7	291	K29* (AAG→TAG) ‡	PP_4079	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (143/130); new base T (6/10); total (149/142)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 306 total aligned reads displayed.

GGACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGC  >  NC_002947/4612259‑4612556

 gACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAg                                                                                                                                                     >  1:4829299/1‑150 (MQ=255)

     aGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTgg                                                                                                                                                 <  2:4719875/150‑1 (MQ=255)

     aGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTgg                                                                                                                                                 <  2:2870860/150‑1 (MQ=255)

       cTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGgcg                                                                                                                                               <  2:1849079/150‑1 (MQ=255)

           gtcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTgc                                                                                                                                           <  2:2525657/150‑2 (MQ=255)

           gtcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTgc                                                                                                                                           <  2:2222553/150‑2 (MQ=255)

           gtcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTgc                                                                                                                                           <  1:3959404/150‑2 (MQ=255)

           gtcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTgg                                                                                                                                           >  1:3806962/1‑150 (MQ=255)

            tcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGa                                                                                                                                          <  2:1690268/150‑1 (MQ=255)

            tcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGa                                                                                                                                          <  1:2497605/150‑1 (MQ=255)

             cgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGAc                                                                                                                                         <  2:1126518/150‑1 (MQ=255)

             cgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCACTGGAc                                                                                                                                         <  2:832947/150‑1 (MQ=255)

              gATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTgctca                                                                                                                                        <  2:5585439/150‑5 (MQ=255)

                 gACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACt                                                                                                                                     <  1:297845/150‑1 (MQ=255)

                 gACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACt                                                                                                                                     <  1:1326661/150‑1 (MQ=255)

                    cACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGt                                                                                                                                  >  1:1122400/1‑150 (MQ=255)

                    cACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGt                                                                                                                                  >  1:5504350/1‑150 (MQ=255)

                     aCGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGtt                                                                                                                                 <  1:6266855/150‑1 (MQ=255)

                       gATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTAc                                                                                                                               <  2:5363638/150‑1 (MQ=255)

                         tGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCa                                                                                                                             >  1:5425620/1‑150 (MQ=255)

                            ccAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATc                                                                                                                          >  1:72537/1‑150 (MQ=255)

                             cAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCg                                                                                                                         >  1:3599530/1‑150 (MQ=255)

                                    gTCCCGTTGGCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGgccgcc                                                                                                                  <  2:5757938/150‑1 (MQ=255)

                                     tCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGccc                                                                                                                 >  2:1254909/1‑150 (MQ=255)

                                      cccGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCa                                                                                                                <  2:6209003/150‑1 (MQ=255)

                                        cGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAg                                                                                                              >  1:5064605/1‑150 (MQ=255)

                                             tCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaa                                                                                                         <  2:3807279/150‑1 (MQ=255)

                                                  ttCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                    >  2:5061985/1‑150 (MQ=255)

                                                  ttCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGAcaacccaa                                                                                                    >  2:607951/1‑150 (MQ=255)

                                                       ggggCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAAcccaacccaa                                                                                               <  2:3764085/150‑1 (MQ=255)

                                                         ggCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGt                                                                                             <  1:3459539/150‑1 (MQ=255)

                                                              tCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGt                                                                                        >  2:4308365/1‑150 (MQ=255)

                                                                aTGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTc                                                                                      <  2:4902460/150‑1 (MQ=255)

                                                                 tGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCa                                                                                     <  1:3019699/150‑1 (MQ=255)

                                                                   aGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGc                                                                                   <  2:366371/150‑1 (MQ=255)

                                                                      cgccTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATc                                                                                <  2:6334921/150‑1 (MQ=255)

                                                                          tGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACg                                                                            >  2:2459727/1‑150 (MQ=255)

                                                                            gCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTc                                                                          >  2:1002350/1‑150 (MQ=255)

                                                                               gTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGc                                                                       >  1:6105462/1‑150 (MQ=255)

                                                                               gTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGc                                                                       >  2:1980561/1‑150 (MQ=255)

                                                                                 cTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGccc                                                                     >  1:3905723/1‑150 (MQ=255)

                                                                                   ggATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGGCCTTGAGCAGAGCCACCCAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCACGTCGCGTTCAGCATCAACGTCAGCCCaa                                                                   <  1:372303/150‑1 (MQ=255)

                                                                                      ttGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTTGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGc                                                                >  1:4688636/1‑150 (MQ=255)

                                                                                      ttGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGc                                                                >  2:3539682/1‑150 (MQ=255)

                                                                                        gcgcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTc                                                              <  2:3176199/150‑1 (MQ=255)

                                                                                         cgcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCa                                                             <  1:1011624/150‑1 (MQ=255)

                                                                                         cgcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGT‑CTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCaa                                                            <  2:1436085/150‑1 (MQ=255)

                                                                                          gcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCaa                                                            >  2:1589263/1‑150 (MQ=255)

                                                                                          gcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCaa                                                            <  1:6386004/150‑1 (MQ=255)

                                                                                          gcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCaa                                                            >  2:6016317/1‑150 (MQ=255)

                                                                                          gcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCaa                                                            <  1:2738039/150‑1 (MQ=255)

                                                                                           ccTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAAc                                                           <  1:501303/150‑1 (MQ=255)

                                                                                            ctctGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACa                                                          >  2:4585972/1‑150 (MQ=255)

                                                                                            ctctGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACa                                                          <  1:4632742/150‑1 (MQ=255)

                                                                                             tctGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAg                                                         >  2:5101219/1‑150 (MQ=255)

                                                                                             tctGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAg                                                         >  2:4342724/1‑150 (MQ=255)

                                                                                                gctggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCaa                                                      <  2:6681080/150‑1 (MQ=255)

                                                                                                 ctggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCacc                                                     >  2:4847230/1‑148 (MQ=255)

                                                                                                 ctggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAAc                                                     >  1:1945698/1‑150 (MQ=255)

                                                                                                 ctggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAAc                                                     >  1:6719577/1‑150 (MQ=255)

                                                                                                  tggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAAcc                                                    <  2:6501095/150‑1 (MQ=255)

                                                                                                    gtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCg                                                  >  2:5917040/1‑150 (MQ=255)

                                                                                                        tgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACa                                              >  2:3424601/1‑150 (MQ=255)

                                                                                                         gctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAg                                             >  2:2289784/1‑150 (MQ=255)

                                                                                                          ctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAgt                                            <  1:6202939/150‑1 (MQ=255)

                                                                                                             gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTc                                         >  2:6867725/1‑150 (MQ=255)

                                                                                                             gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTc                                         <  2:4472537/150‑1 (MQ=255)

                                                                                                                ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGAc                                      >  1:704181/1‑150 (MQ=255)

                                                                                                                ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGAc                                      >  1:1483016/1‑150 (MQ=255)

                                                                                                                      ccAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggcg                                >  2:5209385/1‑150 (MQ=255)

                                                                                                                       cAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggcgg                               >  2:5252139/1‑150 (MQ=255)

                                                                                                                       cAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggcgg                               >  2:4820558/1‑150 (MQ=255)

                                                                                                                       cAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTggcgg                               >  1:2235163/1‑150 (MQ=255)

                                                                                                                          gTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCa                            >  1:4054591/1‑150 (MQ=255)

                                                                                                                            tGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATc                          <  1:3150987/150‑1 (MQ=255)

                                                                                                                                cGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTg                      >  2:1164787/1‑150 (MQ=255)

                                                                                                                                 gTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTga                     >  1:5653268/1‑150 (MQ=255)

                                                                                                                                 gTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTga                     <  2:5832310/150‑1 (MQ=255)

                                                                                                                                 gTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTga                     <  1:2784255/150‑1 (MQ=255)

                                                                                                                                 gTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTga                     <  2:3599530/150‑1 (MQ=255)

                                                                                                                                 gTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTga                     >  1:6325399/1‑150 (MQ=255)

                                                                                                                                   ccTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCACCCCGAACAGTGTCGACGTGGCGGCCATCCCTgaga                   <  2:4965035/150‑1 (MQ=255)

                                                                                                                                   ccTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTgaga                   <  2:4208128/150‑1 (MQ=255)

                                                                                                                                   ccTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTgaga                   <  1:6056792/150‑1 (MQ=255)

                                                                                                                                   ccTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTgaga                   >  2:2649651/1‑150 (MQ=255)

                                                                                                                                    cTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGagag                  <  1:5950483/150‑1 (MQ=255)

                                                                                                                                    cTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGagag                  >  1:772841/1‑150 (MQ=255)

                                                                                                                                     ttGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGc                 >  1:2347284/1‑150 (MQ=255)

                                                                                                                                     ttGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGc                 >  1:5186086/1‑150 (MQ=255)

                                                                                                                                        aGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATc              <  1:846668/150‑1 (MQ=255)

                                                                                                                                         gCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATcc             <  1:4933525/150‑1 (MQ=255)

                                                                                                                                         gCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATcc             >  1:1580919/1‑150 (MQ=255)

                                                                                                                                           agagCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTg           <  2:5348975/150‑1 (MQ=255)

                                                                                                                                           agagCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTg           <  2:1903959/150‑1 (MQ=255)

                                                                                                                                            gcgCCACCTTGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTgg          >  2:3146834/3‑150 (MQ=255)

                                                                                                                                             agCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGa         >  2:1763803/1‑150 (MQ=255)

                                                                                                                                              gCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGaa        <  1:4214481/150‑1 (MQ=255)

                                                                                                                                               ccaccaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAc       >  2:4149166/1‑150 (MQ=255)

                                                                                                                                                 accaAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAAccc     >  2:3014886/1‑150 (MQ=255)

                                                                                                                                                   caAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAg   >  2:4804720/1‑150 (MQ=255)

GGACCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTGGAACCCAGC  >  NC_002947/4612259‑4612556

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