BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	732,491	G→C	12.5%	L676L (CTG→CTC) ‡	clpB →	chaperone protein ClpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	732,491	0	G	C	12.5%	474.3 / 47.6	159	L676L (CTG→CTC) ‡	clpB	chaperone protein ClpB
Reads supporting (aligned to +/- strand): ref base G (75/64); new base C (5/15); total (80/79)
Fisher's exact test for biased strand distribution p-value = 1.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

GAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCA  >  NC_002947/732342‑732640

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

 agGCGATGGTGCGCATCGACATGTCCGAAGTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTgg                                                                                                                                                       <  2:225991/150‑1 (MQ=255)

   gCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGAc                                                                                                                                                     >  1:3393983/1‑150 (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)

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

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

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

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

         gTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTg                                                                                                                                               <  2:3470866/150‑1 (MQ=255)

          tGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTgg                                                                                                                                              <  2:3183670/150‑1 (MQ=255)

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

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

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

                   aCATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGccc                                                                                                                                     <  1:564603/150‑1 (MQ=255)

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

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

                        tCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCGATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCg                                                                                                                                >  1:1511557/1‑150 (MQ=255)

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

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

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

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

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

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

                                  tGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCa                                                                                                                      >  1:4872457/1‑150 (MQ=255)

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

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

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

                                            cTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGtt                                                                                                            <  2:2638010/150‑1 (MQ=255)

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

                                              ccGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGc                                                                                                          >  2:3694092/1‑150 (MQ=255)

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

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

                                                       gcCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTgg                                                                                                 <  2:595353/150‑1 (MQ=255)

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

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

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

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

                                                           gATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGa                                                                                             >  1:3025363/1‑150 (MQ=255)

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

                                                              cGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACgg                                                                                          <  2:3261169/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                            gCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACa                                                                            <  1:4840670/150‑1 (MQ=255)

                                                                             cTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAg                                                                           <  2:3839826/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                              ag‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTgg                                                         >  2:1879422/1‑150 (MQ=255)

                                                                                                 gtcgtttACCTGACCGAGGCAGTGCGTCGCAAGCCTTTCTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACTGCCGCCTGACCGACAGCCTCGGTGCGC‑CCGTggct                                                       <  2:3419683/145‑3 (MQ=37)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                  ccTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAAc                      >  1:254600/1‑150 (MQ=255)

                                                                                                                                     tACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTg                   <  1:2162069/150‑1 (MQ=255)

                                                                                                                                      aCTCGGTGGTGCAGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTgg                  <  2:1370649/150‑1 (MQ=255)

                                                                                                                                         aggtgatCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGG‑CCGTACGGTCGATTTCCGCAACACGGTCATCATCATGACCTCGAACCTCGGcg               >  2:866040/2‑149 (MQ=255)

                                                                                                                                          ggtggtGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTc              <  2:4936704/150‑1 (MQ=255)

                                                                                                                                                gctgctGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCa        <  2:2971601/150‑1 (MQ=255)

                                                                                                                                                gctgctGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCa        <  2:5056906/150‑1 (MQ=255)

                                                                                                                                                     tgGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTc   <  1:4399040/150‑1 (MQ=255)

                                                                                                                                                      ggACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCa  <  2:5045291/150‑1 (MQ=255)

GAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAG‑GGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGG‑GCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCA  >  NC_002947/732342‑732640

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