BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	732,469	A→G	10.1%	K669R (AAG→AGG) ‡	clpB →	chaperone protein ClpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	732,469	0	A	G	10.1%	509.4 / 33.9	178	K669R (AAG→AGG) ‡	clpB	chaperone protein ClpB
Reads supporting (aligned to +/- strand): ref base A (87/73); new base G (4/14); total (91/87)
Fisher's exact test for biased strand distribution p-value = 1.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 179 total aligned reads displayed.

AGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCC  >  NC_002947/732322‑732617

  ttCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCt                                                                                                                                                   >  2:923043/1‑150 (MQ=255)

  ttCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCt                                                                                                                                                   <  2:2253484/150‑1 (MQ=255)

  ttCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCt                                                                                                                                                   <  1:5071062/150‑1 (MQ=255)

  ttCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCt                                                                                                                                                   <  1:4203569/150‑1 (MQ=255)

     cTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTAc                                                                                                                                                >  2:1554488/1‑150 (MQ=255)

      tGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACt                                                                                                                                               >  1:1343606/1‑150 (MQ=255)

      tGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACt                                                                                                                                               >  2:4399040/1‑150 (MQ=255)

      tGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACt                                                                                                                                               >  2:4160058/1‑150 (MQ=255)

      tGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACt                                                                                                                                               >  1:2265685/1‑150 (MQ=255)

        ttCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCg                                                                                                                                             <  2:4557516/150‑1 (MQ=255)

        ttCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCg                                                                                                                                             >  1:3419804/1‑150 (MQ=255)

        ttCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAAATCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAATAG‑GGCGGTTACCTGACCGAGTCCGTGCGGCGCAAGCCTTACTCg                                                                                                                                             >  1:3419683/1‑150 (MQ=255)

         tCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCgg                                                                                                                                            <  1:2925512/150‑1 (MQ=255)

           gACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTAGGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCggtg                                                                                                                                          <  2:2238062/150‑1 (MQ=255)

              acCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGgtggtg                                                                                                                                       >  2:139721/1‑150 (MQ=255)

                 gAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGtgctg                                                                                                                                    <  1:626792/150‑1 (MQ=255)

                   agagGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTgctgct                                                                                                                                  <  1:2731434/150‑1 (MQ=255)

                   agagGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTgctgct                                                                                                                                  <  1:4859703/150‑1 (MQ=255)

                    gagGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGctgctg                                                                                                                                 >  2:3729258/1‑150 (MQ=255)

                    gagGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGctgctg                                                                                                                                 <  2:3992806/150‑1 (MQ=255)

                       gCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGAc                                                                                                                              <  2:3794680/150‑1 (MQ=255)

                         gATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGa                                                                                                                            >  2:2162069/1‑150 (MQ=255)

                         gATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGa                                                                                                                            <  1:3298093/150‑1 (MQ=255)

                           tGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAgg                                                                                                                          <  2:1723059/150‑1 (MQ=255)

                           tGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAgg                                                                                                                          <  2:573633/150‑1 (MQ=255)

                           tGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCGTACTCGGTGGTGCTGCTGGACGAgg                                                                                                                          <  1:3655168/150‑1 (MQ=255)

                               gcgcATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGa                                                                                                                      >  2:571949/1‑150 (MQ=255)

                                   aTCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGaaaaa                                                                                                                  <  1:3973111/150‑1 (MQ=255)

                                   aTCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGaaaaa                                                                                                                  <  2:309701/150‑1 (MQ=255)

                                           gTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGG‑AAAAGCCCACCCg                                                                                                         >  2:1681878/1‑150 (MQ=255)

                                            tCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCg                                                                                                         >  1:4327562/1‑150 (MQ=255)

                                              cGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGa                                                                                                       >  2:3915543/1‑150 (MQ=255)

                                               gAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGAt                                                                                                      >  1:2520963/1‑150 (MQ=255)

                                                     aTGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAGGCCCACCCGGATGTGTTc                                                                                                >  2:116913/1‑150 (MQ=255)

                                                        gagaAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAAc                                                                                             >  2:262133/1‑150 (MQ=255)

                                                            aGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCAGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGc                                                                                         <  2:1938107/150‑1 (MQ=255)

                                                            aGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCAGCGCAAGCCGTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGc                                                                                         <  2:1938059/150‑1 (MQ=255)

                                                             gCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCt                                                                                        <  2:4692978/150‑1 (MQ=255)

                                                                 tCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTg                                                                                    >  1:5057831/1‑150 (MQ=255)

                                                                 tCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTg                                                                                    >  2:2187866/1‑150 (MQ=255)

                                                                      cgcgcgCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGt                                                                               <  2:2694433/150‑1 (MQ=255)

                                                                        gccgtcTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAG‑GGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGttc                                                                             <  1:2539504/145‑3 (MQ=255)

                                                                         ccgacTGATTGGTGCACCTCCCGTATATGTCGGCTATGAACAGCGTCGCGTA‑CTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCt                                                                            <  2:3408224/146‑1 (MQ=255)

                                                                          cgtcTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAG‑GGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCtt                                                                           <  2:1006418/147‑2 (MQ=255)

                                                                           gcCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTgg                                                                          >  2:4397567/1‑150 (MQ=255)

                                                                            tcTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAG‑GGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGa                                                                         <  2:2857147/149‑1 (MQ=255)

                                                                             cTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAG‑GGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGat                                                                        <  1:2354814/150‑2 (MQ=255)

                                                                             cTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGaa                                                                        >  2:504895/1‑150 (MQ=255)

                                                                             cTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGaa                                                                        >  2:3059935/1‑150 (MQ=255)

                                                                               gATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAG‑GGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTtgatga                                                                      >  1:599450/1‑150 (MQ=255)

                                                                                 tCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACg                                                                    <  2:4856873/150‑1 (MQ=255)

                                                                                   ggTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgc                                                                  >  2:148460/1‑150 (MQ=255)

                                                                                    gTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgcc                                                                 >  2:2042496/1‑150 (MQ=255)

                                                                                    gTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgcc                                                                 <  1:5041745/150‑1 (MQ=255)

                                                                                      gCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgccgc                                                               <  2:4337168/150‑1 (MQ=255)

                                                                                      gCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgccgc                                                               <  2:3209446/150‑1 (MQ=255)

                                                                                       caccacCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgccgcc                                                              >  2:708257/1‑150 (MQ=255)

                                                                                          caccaGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGa                                                           >  1:3386690/1‑150 (MQ=255)

                                                                                          caccaGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGa                                                           >  2:4945908/1‑150 (MQ=255)

                                                                                           accaGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGAc                                                          <  1:2463938/150‑1 (MQ=255)

                                                                                             ccggttATGTCGGCTATGAAGAG‑GGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCg                                                        >  1:3994997/6‑150 (MQ=255)

                                                                                                  tATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGc                                                   >  2:4581149/1‑150 (MQ=255)

                                                                                                    tGTCGGCTATGAAGAG‑GGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGAcggtca                                                 <  1:3067712/150‑6 (MQ=255)

                                                                                                       aGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACgg                                              <  2:2527804/150‑1 (MQ=255)

                                                                                                        gggTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACgg‑g                                            <  2:21374/150‑1 (MQ=255)

                                                                                                        gggTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACgg‑g                                            >  2:3700267/1‑150 (MQ=255)

                                                                                                         ggTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑gc                                           >  1:4562179/1‑150 (MQ=255)

                                                                                                           tATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑gcgc                                         >  1:3155617/1‑150 (MQ=255)

                                                                                                            aTGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCa                                        >  1:912260/1‑150 (MQ=255)

                                                                                                              gAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCAcc                                      >  2:3625296/1‑150 (MQ=255)

                                                                                                                 gag‑GGCTGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTg                                   <  1:1720666/150‑1 (MQ=255)

                                                                                                                 gag‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTg                                   <  1:1720542/150‑1 (MQ=255)

                                                                                                                 gag‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTg                                   <  1:212111/150‑1 (MQ=255)

                                                                                                                     ggCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGTGCGC‑CCGTggct                                <  2:3419804/150‑3 (MQ=255)

                                                                                                                         gcgttcTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGG‑CCGTACGGTCGATTTcc                            <  2:1928968/145‑1 (MQ=255)

                                                                                                                         gTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTcc                            <  2:1860933/150‑1 (MQ=255)

                                                                                                                          cgttcTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGG‑CCGTACGGTCGATTTCCg                           <  1:4710943/146‑1 (MQ=255)

                                                                                                                            aCCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCa                         >  2:2079319/1‑150 (MQ=255)

                                                                                                                             ccTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCaa                        >  1:4714913/1‑150 (MQ=255)

                                                                                                                              cTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAac                       >  2:1796831/1‑150 (MQ=255)

                                                                                                                              cTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAac                       >  1:259370/1‑150 (MQ=255)

                                                                                                                              cTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAac                       >  1:3920785/1‑150 (MQ=255)

                                                                                                                                gACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAacac                     >  1:3154591/1‑150 (MQ=255)

                                                                                                                                gACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAacac                     >  1:896612/1‑150 (MQ=255)

                                                                                                                                  ccGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcg                   >  1:823569/1‑150 (MQ=255)

                                                                                                                                   cGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcgt                  >  2:2905976/1‑150 (MQ=255)

                                                                                                                                   cGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcgt                  >  2:2345434/1‑150 (MQ=255)

                                                                                                                                   cGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcgt                  <  2:1343606/150‑1 (MQ=255)

                                                                                                                                     aGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcgtga                >  1:1129830/1‑150 (MQ=255)

                                                                                                                                       gCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcgtgatc              <  1:1471410/150‑1 (MQ=255)

                                                                                                                                            gcggcgCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcgtgatcgtgat         <  1:4976058/150‑1 (MQ=255)

                                                                                                                                               gcgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGa                                                                                                       >  1:1553147/1‑54 (MQ=255)

                                                                                                                                               gcgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGa                                                                                                       <  2:1553147/54‑1 (MQ=255)

                                                                                                                                               gcgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGAc      <  1:285762/150‑1 (MQ=255)

                                                                                                                                               gcgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGAc      >  2:823378/1‑150 (MQ=255)

                                                                                                                                               gcgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGAc      >  2:33188/1‑150 (MQ=255)

                                                                                                                                                 gcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCt    <  1:4012391/150‑1 (MQ=255)

                                                                                                                                                 gcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCt    >  1:2698472/1‑150 (MQ=255)

                                                                                                                                                 gcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCt    >  1:2698313/1‑150 (MQ=255)

                                                                                                                                                   aaGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTcc  >  2:2274709/1‑150 (MQ=255)

AGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCC  >  NC_002947/732322‑732617

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