BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,654,325	C→A	5.3%	E576* (GAG→TAG) ‡	PP_1446 ←	TonB‑dependent receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,654,325	0	C	A	5.3%	1273.0 / 42.0	374	E576* (GAG→TAG) ‡	PP_1446	TonB‑dependent receptor
Reads supporting (aligned to +/- strand): ref base C (187/167); new base A (7/13); total (194/180)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 381 total aligned reads displayed.

CAAGGTGGTGTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGC  >  NC_002947/1654176‑1654474

cAAGGTGGTGTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTc                                                                                                                                                       <  2:2126175/150‑1 (MQ=255)

  aGGTGGTGTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGc                                                                                                                                                     <  2:3205814/150‑1 (MQ=255)

  aGGTGGTGTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGc                                                                                                                                                     >  2:7022481/1‑150 (MQ=255)

       gtgtAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTg                                                                                                                                                <  1:7090859/150‑1 (MQ=255)

        tgtAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGt                                                                                                                                               <  1:5615551/150‑1 (MQ=255)

        tgtAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGACAGTGGCGAGCCACTgc                                                                                                                                               <  2:3051832/150‑2 (MQ=255)

         gtAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTc                                                                                                                                              >  2:4160283/1‑150 (MQ=255)

         gtAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTc                                                                                                                                              >  1:6520591/1‑150 (MQ=255)

         gtAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCGAGCCACTgcc                                                                                                                                              <  2:8176735/150‑3 (MQ=255)

         gtAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGACAGTGGCGAGCCACTgcc                                                                                                                                              <  2:837992/150‑3 (MQ=255)

          tAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCGAGCCACTgccc                                                                                                                                             <  2:8139569/150‑4 (MQ=255)

            gCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAg                                                                                                                                           <  2:1213138/150‑1 (MQ=255)

            gCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAg                                                                                                                                           <  1:7226157/150‑1 (MQ=255)

            gCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAg                                                                                                                                           >  1:1213138/1‑150 (MQ=255)

             ccATCGGTGCTGGTTTCGTTGGCGGCTTTGCGGGGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCCGGCCGCTGTTCAGGCGCAGCGGGGCGATGAGGGGAAGAGGCTCGCCCCAGTCGAgg                                                                                                                                          >  1:6726160/1‑150 (MQ=255)

                  ggTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCtt                                                                                                                                     <  2:2432927/150‑1 (MQ=255)

                    tGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTgg                                                                                                                                   >  2:1677130/1‑150 (MQ=255)

                    tGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTgg                                                                                                                                   >  1:2764764/1‑150 (MQ=255)

                      cTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCt                                                                                                                                 >  1:2273770/1‑150 (MQ=255)

                         gTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCgg                                                                                                                              >  1:6869106/1‑150 (MQ=255)

                         gTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCgg                                                                                                                              <  1:7505309/150‑1 (MQ=255)

                         gTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCgg                                                                                                                              <  1:1492780/150‑1 (MQ=255)

                         gTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCgg                                                                                                                              <  2:4883186/150‑1 (MQ=255)

                         gTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCgg                                                                                                                              >  2:8028284/1‑150 (MQ=255)

                          tttCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCggg                                                                                                                             <  1:8171782/150‑1 (MQ=255)

                          tttCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCggg                                                                                                                             >  1:5950124/1‑150 (MQ=255)

                             cGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGtat                                                                                                                          <  2:4547955/150‑1 (MQ=255)

                              gTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGtata                                                                                                                         >  1:5659490/1‑150 (MQ=255)

                              gTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGtata                                                                                                                         <  2:4945627/150‑1 (MQ=255)

                              gTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGtata                                                                                                                         <  1:7791109/150‑1 (MQ=255)

                              gTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGtata                                                                                                                         <  2:3807176/150‑1 (MQ=255)

                              gTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGtata                                                                                                                         <  2:391924/150‑1 (MQ=255)

                              gTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGtata                                                                                                                         <  2:7244513/150‑1 (MQ=255)

                              gTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGtata                                                                                                                         >  1:587402/1‑150 (MQ=255)

                               ttGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAg                                                                                                                        >  2:3696314/1‑150 (MQ=255)

                                  gCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCa                                                                                                                     <  1:2578151/150‑1 (MQ=255)

                                   cGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCAc                                                                                                                    <  1:8025944/150‑1 (MQ=255)

                                     ggTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCg                                                                                                                  >  1:4526041/1‑150 (MQ=255)

                                     ggTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCg                                                                                                                  <  1:2197845/150‑1 (MQ=255)

                                       tttGCGGTGCTGCGACGCAGCGTGCTGCACGTCTCCCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGa                                                                                                                <  2:3170809/150‑1 (MQ=255)

                                         tGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACa                                                                                                              >  2:2839405/1‑150 (MQ=255)

                                            ggTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACt                                                                                                           <  2:3780347/150‑1 (MQ=255)

                                                     gtgCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGa                                                                                                  >  2:3379735/3‑150 (MQ=255)

                                                         aGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCt                                                                                              <  2:6741906/150‑1 (MQ=255)

                                                            gTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGcc                                                                                           <  2:6277776/150‑1 (MQ=255)

                                                             tgctgcACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCg                                                                                          <  2:5462435/150‑1 (MQ=255)

                                                             tgctgcACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCg                                                                                          <  2:6456938/150‑1 (MQ=255)

                                                                tgcACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGGTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTaa                                                                                       >  1:2109633/1‑150 (MQ=255)

                                                                        tACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTtctc                                                                               >  1:945648/1‑150 (MQ=255)

                                                                        tACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTtctc                                                                               <  2:5641285/150‑1 (MQ=255)

                                                                         aCCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCg                                                                              <  2:5308300/150‑1 (MQ=255)

                                                                            cgcgCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTACTTGGCTCGGGTAAAGTAACCGGACAACTCCAGCGCGAAGCTGCCTTAACGGTTCTCGGcc                                                                           >  1:5242442/1‑150 (MQ=255)

                                                                             gcgcCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACg                                                                                     >  2:7872017/1‑139 (MQ=255)

                                                                             gcgcCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCa                                                                          >  2:2009042/1‑150 (MQ=255)

                                                                              cgcCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCaa                                                                         <  1:8787767/150‑1 (MQ=255)

                                                                                ccTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACt                                                                       >  1:8477538/1‑150 (MQ=255)

                                                                                    ccAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCa                                                                   >  2:6540150/1‑150 (MQ=255)

                                                                                    ccAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCa                                                                   <  1:1677130/150‑1 (MQ=255)

                                                                                      aGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGc                                                                 >  2:2270319/1‑150 (MQ=255)

                                                                                          gTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTc                                                             >  1:4656776/1‑150 (MQ=255)

                                                                                             cAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTg                                                          >  2:5229288/1‑150 (MQ=255)

                                                                                             cAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTg                                                          >  2:7459043/1‑150 (MQ=255)

                                                                                             cAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTg                                                          >  1:2957882/1‑150 (MQ=255)

                                                                                              aGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACGGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCACCTGCCCGCGGTCCTgg                                                         >  2:6573687/1‑150 (MQ=255)

                                                                                               gCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTggg                                                        <  1:7747571/150‑1 (MQ=255)

                                                                                               gCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTggg                                                        >  2:7771126/1‑150 (MQ=255)

                                                                                                 tCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGcc                                                      >  1:4276836/1‑150 (MQ=255)

                                                                                                 tCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGcc                                                      >  2:807416/1‑150 (MQ=255)

                                                                                                  cccAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCt                                                     >  2:573679/1‑150 (MQ=255)

                                                                                                  cccAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCt                                                     <  2:1141018/150‑1 (MQ=255)

                                                                                                  cccAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCt                                                     >  1:3335408/1‑150 (MQ=255)

                                                                                                   ccAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTc                                                    >  2:347415/1‑150 (MQ=255)

                                                                                                   ccAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTc                                                    >  1:6321939/1‑150 (MQ=255)

                                                                                                   ccAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAAAGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTc                                                    <  2:8563339/150‑1 (MQ=255)

                                                                                                    cAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCg                                                   >  1:6635997/1‑150 (MQ=255)

                                                                                                        ccAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGc                                               <  1:8363875/150‑1 (MQ=255)

                                                                                                         cAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGcc                                              >  1:3802732/1‑150 (MQ=255)

                                                                                                          aGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCg                                             <  1:1836427/150‑1 (MQ=255)

                                                                                                            gCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGAc                                                                                                               >  2:4823323/1‑82 (MQ=255)

                                                                                                            gCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGAc                                                                                                               <  1:4823323/82‑1 (MQ=255)

                                                                                                            gCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTa                                           >  1:6957542/1‑150 (MQ=255)

                                                                                                            gCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTa                                           <  1:3098230/150‑1 (MQ=255)

                                                                                                            gCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTa                                           >  2:2604754/1‑150 (MQ=255)

                                                                                                             ccGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAg                                          <  2:5153282/150‑1 (MQ=255)

                                                                                                                 tGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAgc                                      <  2:7619043/150‑1 (MQ=255)

                                                                                                                   ttCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAgcgc                                    >  1:863016/1‑150 (MQ=255)

                                                                                                                   ttCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAgcgc                                    >  2:5883702/1‑150 (MQ=255)

                                                                                                                   ttCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAgcgc                                    >  2:1874288/1‑150 (MQ=255)

                                                                                                                        gcgcAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACCACTCCAGCGCGAAGCTGCCGTAACGGGTCTCGGCCCACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCg                               >  2:6198278/1‑150 (MQ=255)

                                                                                                                            aGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTTTAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACg                           <  2:4287679/150‑1 (MQ=255)

                                                                                                                             gCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGc                          >  2:4640063/1‑150 (MQ=255)

                                                                                                                               ggggCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCt                        <  1:7575408/150‑1 (MQ=255)

                                                                                                                               ggggCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCt                        <  1:2410137/150‑1 (MQ=255)

                                                                                                                                  gCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTgg                     <  1:7359119/150‑1 (MQ=255)

                                                                                                                                      tGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACt                 >  2:504579/1‑150 (MQ=255)

                                                                                                                                         tcgacAGTGGCGAGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGt              >  1:3978602/5‑150 (MQ=255)

                                                                                                                                          gggCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTa             >  1:4416968/1‑150 (MQ=255)

                                                                                                                                          cgacAGTGGCGAGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTa             >  2:7336555/4‑150 (MQ=255)

                                                                                                                                              agaaccTCGACAGTGGCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTcc         >  1:6235486/6‑150 (MQ=255)

                                                                                                                                                tGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCgg       <  1:5009328/150‑1 (MQ=255)

CAAGGTGGTGTAGCCATCGGTGCTGGTTTCGTTGGCGGGTTTGCGGTGCTGCGACGCAGCGTGCTGCACGTCTACCCGCGCCTGCCAGCGGTCCAGCTCCCAGACCAGGCCGCTGTTCAGGCGCAGCGGGGCGATGCGGGGCAGTGGCTCGCCACTGTCGAGGTTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGC  >  NC_002947/1654176‑1654474

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