BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	732,476	C→T	11.8%	Y671Y (TAC→TAT)	clpB →	chaperone protein ClpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	732,476	0	C	T	11.8%	669.6 / 62.5	246	Y671Y (TAC→TAT)	clpB	chaperone protein ClpB
Reads supporting (aligned to +/- strand): ref base C (115/102); new base T (11/18); total (126/120)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 251 total aligned reads displayed.

CTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGG  >  NC_002947/732327‑732625

cTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTAc                                                                                                                                                        >  1:8577412/1‑150 (MQ=255)

cTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTAc                                                                                                                                                        <  2:3019948/150‑1 (MQ=255)

cTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTAc                                                                                                                                                        <  1:3184286/150‑1 (MQ=255)

 tGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGTTCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGTGGCCGTGCGGCGCAAGCCTTACt                                                                                                                                                       >  2:7466987/1‑150 (MQ=255)

     cGACACCGAAGAGGCGATGGTGCGCATCGTCATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCggt                                                                                                                                                   >  1:4850753/1‑150 (MQ=255)

      gACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCggtg                                                                                                                                                  <  1:4659799/150‑1 (MQ=255)

      gACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCggtg                                                                                                                                                  <  2:2170430/150‑1 (MQ=255)

      gACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCggtg                                                                                                                                                  <  2:2691074/150‑1 (MQ=255)

         acCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGgtggtg                                                                                                                                               >  2:1668960/1‑150 (MQ=255)

           cGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGtgct                                                                                                                                             >  1:6304960/1‑150 (MQ=255)

           cGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGtgct                                                                                                                                             >  1:212283/1‑150 (MQ=255)

            gAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGtgctg                                                                                                                                            >  2:5323297/1‑150 (MQ=255)

              agagGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTgctgct                                                                                                                                          <  2:2264429/150‑1 (MQ=255)

                agGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTgg                                                                                                                                        <  2:8577412/150‑1 (MQ=255)

                  gCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGAc                                                                                                                                      >  2:4673659/1‑150 (MQ=255)

                     aTGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAg                                                                                                                                   <  2:8128938/150‑1 (MQ=255)

                     aTGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAg                                                                                                                                   <  1:2353201/150‑1 (MQ=255)

                      tGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAgg                                                                                                                                  <  2:8894103/150‑1 (MQ=255)

                      tGGTGCGCAACGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAgg                                                                                                                                  <  2:8347083/150‑1 (MQ=255)

                         tGCGCGTCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTgg                                                                                                                               <  2:5033179/150‑1 (MQ=255)

                         tGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTgg                                                                                                                               >  2:7629523/1‑150 (MQ=255)

                            acATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGaaa                                                                                                                            <  2:8637448/149‑1 (MQ=255)

                             cATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGaaaa                                                                                                                           >  1:8873738/1‑150 (MQ=255)

                               tCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAg                                                                                                                         <  1:7629523/150‑1 (MQ=255)

                                    aTGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCAc                                                                                                                    >  2:2279442/1‑150 (MQ=255)

                                       tCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCg                                                                                                                 >  1:4589575/1‑150 (MQ=255)

                                       tCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCg                                                                                                                 >  1:357323/1‑150 (MQ=255)

                                         cGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGa                                                                                                               >  1:419801/1‑150 (MQ=255)

                                         cGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGa                                                                                                               >  1:3499686/1‑150 (MQ=255)

                                          gAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGAt                                                                                                              <  1:1076046/150‑1 (MQ=255)

                                              tCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATgtgt                                                                                                          >  2:2233501/1‑150 (MQ=255)

                                                aTGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTc                                                                                                                                                      <  1:361252/104‑1 (MQ=255)

                                                aTGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTc                                                                                                        <  1:874444/150‑1 (MQ=255)

                                                aTGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTc                                                                                                        >  1:1742086/1‑150 (MQ=255)

                                                 tGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCa                                                                                                       >  2:6434251/1‑150 (MQ=255)

                                                   gagaAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAAc                                                                                                     >  2:495342/1‑150 (MQ=255)

                                                   gagaAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAAc                                                                                                     <  1:6470458/150‑1 (MQ=255)

                                                      aaGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTg                                                                                                  >  1:2614323/1‑150 (MQ=255)

                                                        gCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCt                                                                                                >  1:7049506/1‑150 (MQ=255)

                                                              cGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCa                                                                                          <  2:7191334/150‑1 (MQ=255)

                                                              cGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCa                                                                                          <  1:5157677/150‑1 (MQ=255)

                                                                 gagccgtcTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGt                                                                                       <  1:9162723/143‑1 (MQ=255)

                                                                  gcgcgcCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTg                                                                                      <  2:3750873/150‑1 (MQ=255)

                                                                    gcgcCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCt                                                                                    <  1:6095662/150‑1 (MQ=255)

                                                                     cgtcTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCtt                                                                                   <  1:5411381/147‑2 (MQ=255)

                                                                       tcTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGa                                                                                 >  2:3416529/2‑150 (MQ=255)

                                                                        cTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGaa                                                                                <  1:2369348/150‑1 (MQ=255)

                                                                           aTCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGAc                                                                             <  1:5560604/150‑1 (MQ=255)

                                                                              ggagCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTAAACGTGCTGTTGCAGGTGCTGGAAGACGgc                                                                          >  1:2783514/4‑150 (MQ=255)

                                                                               gTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACggtc                                                                         <  2:5980784/150‑3 (MQ=255)

                                                                               gTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgcc                                                                         >  2:4048676/1‑150 (MQ=255)

                                                                                 gCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgccgc                                                                       <  2:1751409/150‑1 (MQ=255)

                                                                                      tcccggttATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                                                  >  1:5040665/8‑150 (MQ=255)

                                                                                      accaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGAc                                                                  >  2:2872114/1‑150 (MQ=255)

                                                                                       ccaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGAcc                                                                 <  2:2963930/150‑1 (MQ=255)

                                                                                        caGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCg                                                                >  1:4069668/1‑150 (MQ=255)

                                                                                         aGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGa                                                               <  1:8373731/150‑1 (MQ=255)

                                                                                          ggCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGAc                                                              >  1:1260647/1‑150 (MQ=255)

                                                                                           gttATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGAcg                                                             <  1:8600650/148‑2 (MQ=255)

                                                                                           gCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACa                                                             >  1:2669523/1‑150 (MQ=255)

                                                                                                 tAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACg                                                       >  2:8823311/1‑150 (MQ=255)

                                                                                                  aGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCTGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACgg                                                      <  2:342868/150‑1 (MQ=255)

                                                                                                  aGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACgg                                                      >  2:5740745/1‑150 (MQ=255)

                                                                                                  aGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACgg                                                      <  1:7884358/150‑1 (MQ=255)

                                                                                                   gggTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACgg‑g                                                    >  2:4413209/1‑150 (MQ=255)

                                                                                                   ggctATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                                                  >  2:7825273/4‑137 (MQ=255)

                                                                                                     gTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACCGCCACGG‑gcg                                                  >  2:3900902/1‑150 (MQ=255)

                                                                                                     gTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑gcg                                                  >  2:4725646/1‑150 (MQ=255)

                                                                                                       aTGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCa                                                <  1:4354373/150‑1 (MQ=255)

                                                                                                        tGAAGAGGTCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCAc                                               >  1:8197970/1‑150 (MQ=255)

                                                                                                        tGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCAc                                               >  1:1728307/1‑150 (MQ=255)

                                                                                                         gAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCAcc                                              >  1:3048433/1‑150 (MQ=255)

                                                                                                          aaGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCg                                             >  2:8734498/1‑150 (MQ=255)

                                                                                                            gagGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTg                                           >  1:3684047/1‑150 (MQ=255)

                                                                                                             agGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGG‑AAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGTGCGC‑CCGTggc                                         <  1:1520500/150‑2 (MQ=255)

                                                                                                              gggCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGa                                         <  2:360506/150‑1 (MQ=255)

                                                                                                               ggCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGACAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGTGCGC‑CCGTggct                                        <  2:5646240/150‑3 (MQ=255)

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

                                                                                                                   gTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTcc                                    >  2:2480633/1‑150 (MQ=255)

                                                                                                                   gTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTcc                                    >  2:6873872/1‑150 (MQ=255)

                                                                                                                    ttACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCg                                   >  2:4451493/1‑150 (MQ=255)

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

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

                                                                                                                       tcTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGG‑CCGTACGGTCGATTTCCGCaa                                >  1:5331113/2‑150 (MQ=255)

                                                                                                                       ccTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCaa                                <  2:357332/150‑1 (MQ=255)

                                                                                                                         tGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAaca                              <  1:6653072/150‑1 (MQ=255)

                                                                                                                         tGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAaca                              <  2:6724865/150‑1 (MQ=255)

                                                                                                                           aCCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACAcc                            >  1:7286880/1‑150 (MQ=255)

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

                                                                                                                              gAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcgtg                         >  1:2508804/1‑150 (MQ=255)

                                                                                                                                     tGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACcgtgatcgtga                  <  2:5481138/150‑1 (MQ=255)

                                                                                                                                        ggcgCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGtgatga               <  2:8671526/150‑1 (MQ=255)

                                                                                                                                         gcgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGCCTTCCGCAACACCGTGATCGTGATGAc              <  2:1465243/150‑1 (MQ=255)

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

                                                                                                                                          cgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGAcc             >  1:8068007/1‑150 (MQ=255)

                                                                                                                                            cAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTc           >  2:5554078/1‑150 (MQ=255)

                                                                                                                                               gCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCaa        >  1:7448634/1‑150 (MQ=255)

                                                                                                                                                 cTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAAcc      >  2:110466/1‑150 (MQ=255)

                                                                                                                                                     cTCGGTGGTGTTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTggg  <  2:4850753/150‑1 (MQ=255)

                                                                                                                                                     cTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTggg  >  1:1339939/1‑150 (MQ=255)

CTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGG  >  NC_002947/732327‑732625

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