BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,977,219	T→C	H59H (CAT→CAC)	lysR →	LysR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,977,219	0	T	C	99.3%	934.2 / ‑3.0	296	H59H (CAT→CAC)	lysR	LysR
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (135/159); total (137/159)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-02

Only 100 of 305 total aligned reads displayed.

CACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTGC  >  NC_000913/2977129‑2977307

cacCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACAccc                                                                                        >  2:2267857/1‑93 (MQ=255)

cacCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACAccc                                                                                        <  1:815419/93‑1 (MQ=255)

cacCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACAccc                                                                                        <  2:4249728/93‑1 (MQ=255)

cacCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACAccc                                                                                        <  1:6797906/93‑1 (MQ=255)

  ccTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTa                                                                                      >  2:2918529/1‑93 (MQ=255)

   cTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTAc                                                                                     <  2:2061095/93‑1 (MQ=255)

   cTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTAc                                                                                     <  2:5512368/93‑1 (MQ=255)

    tCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTAcc                                                                                    <  2:7360117/93‑1 (MQ=255)

    tCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTAcc                                                                                    >  1:2832744/1‑93 (MQ=255)

        aGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGc                                                                                >  2:3764350/1‑93 (MQ=255)

         gCCAACCGTCAGCCGCGAACTTGCGCACTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCa                                                                               >  1:3020488/1‑93 (MQ=255)

          ccAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCaa                                                                              <  2:530455/93‑1 (MQ=255)

           cAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAg                                                                             <  1:6556507/93‑1 (MQ=255)

             aCCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAggc                                                                           >  1:1682418/1‑92 (MQ=255)

             aCCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGa                                                                           <  1:5797581/93‑1 (MQ=255)

              ccGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGAc                                                                          <  1:7565168/93‑1 (MQ=255)

              ccGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGAc                                                                          <  2:3406201/93‑1 (MQ=255)

                 tCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGc                                                                       <  2:7299651/93‑1 (MQ=255)

                    gCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTc                                                                    <  1:2482978/93‑1 (MQ=255)

                     ccGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCt                                                                   <  1:2711895/93‑1 (MQ=255)

                       gcgAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGt                                                                 >  2:7195556/1‑93 (MQ=255)

                       ccgAACTTGCGCGCTTTGAGAAGGTGATCAGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGt                                                                 <  2:1886938/92‑1 (MQ=255)

                        cgAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGtt                                                                <  2:1511545/93‑1 (MQ=255)

                          aaCTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTg                                                              <  1:3901530/93‑1 (MQ=255)

                             ttGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTgaag                                                           <  2:4108589/93‑1 (MQ=255)

                              tGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTgaaga                                                          <  1:3750271/93‑1 (MQ=255)

                               gcgcgcTTTGAGAATGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTgaagaa                                                         <  2:5974037/93‑1 (MQ=255)

                               gcgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTgaagaa                                                         <  2:4568037/93‑1 (MQ=255)

                                cgcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGaagaag                                                        <  2:249240/93‑1 (MQ=255)

                                 gcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGt                                                       >  2:3020038/1‑93 (MQ=255)

                                  cgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTg                                                      <  1:6222128/93‑1 (MQ=255)

                                   gcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGc                                                     >  1:7813235/1‑93 (MQ=255)

                                   gcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGc                                                     >  2:1099273/1‑93 (MQ=255)

                                     tttGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCaa                                                   >  1:901295/1‑93 (MQ=255)

                                     tttGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCaa                                                   <  2:5425067/93‑1 (MQ=255)

                                     tttGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCaa                                                   <  2:3966237/93‑1 (MQ=255)

                                      ttGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAAc                                                  <  1:1830755/93‑1 (MQ=255)

                                      ttGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAAc                                                  <  1:3333269/93‑1 (MQ=255)

                                       tGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACg                                                 >  2:6951931/1‑93 (MQ=255)

                                        gagaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGa                                                >  2:4634902/1‑93 (MQ=255)

                                        gagaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGa                                                <  1:4142103/93‑1 (MQ=255)

                                         agaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGAt                                               <  1:5432020/93‑1 (MQ=255)

                                         agaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGAt                                               <  1:1458654/93‑1 (MQ=255)

                                           aaGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATcc                                             <  1:2291090/93‑1 (MQ=255)

                                             ggTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTg                                           <  1:4583427/93‑1 (MQ=255)

                                             ggTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTg                                           <  1:3290855/93‑1 (MQ=255)

                                             ggTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTg                                           <  2:2375725/93‑1 (MQ=255)

                                             ggTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTg                                           <  1:6324250/93‑1 (MQ=255)

                                              gTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTgg                                          <  1:4510088/93‑1 (MQ=255)

                                                gATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTa                                        <  1:3808271/93‑1 (MQ=255)

                                                gATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTa                                        >  1:5375864/1‑93 (MQ=255)

                                                gATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATACTGGTa                                        <  2:7679380/93‑1 (MQ=255)

                                                  tCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACg                                      <  2:3520079/93‑1 (MQ=255)

                                                  tCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACg                                      >  1:3208287/1‑93 (MQ=255)

                                                   cGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACgg                                     <  1:138474/93‑1 (MQ=255)

                                                    gggCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGa                                    >  1:7208909/1‑93 (MQ=255)

                                                    gggCTGAAATTGTTTGAGCGCGTACGTGAGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACggt                                    <  1:7045439/93‑2 (MQ=255)

                                                      gCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGGCTGCGTCTGTTTGAAGAAGTGCAAAGATACTGGTAGGGACt                                  >  1:2941231/1‑93 (MQ=255)

                                                       cTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTg                                 <  1:1542646/93‑1 (MQ=255)

                                                         gAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGa                               >  1:1969286/1‑93 (MQ=255)

                                                         gAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGa                               >  2:6510994/1‑93 (MQ=255)

                                                          aaaTTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGAt                              <  2:5601962/93‑1 (MQ=255)

                                                          aaaTTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGAt                              <  1:5833954/93‑1 (MQ=255)

                                                           aaTTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATc                             <  1:617941/93‑1 (MQ=255)

                                                           aaTTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATc                             >  1:7246429/1‑93 (MQ=255)

                                                            aTTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCg                            >  2:5749895/1‑93 (MQ=255)

                                                             ttGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGc                           >  2:1337359/1‑93 (MQ=255)

                                                             ttGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGc                           >  2:5817335/1‑93 (MQ=255)

                                                             ttGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGc                           >  1:4391285/1‑93 (MQ=255)

                                                             ttGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGc                           <  2:7479390/93‑1 (MQ=255)

                                                             atGTTTGGGCGCGTTCGGGGGCGATTACCCCCCACCGTGCAAGGACTGGGTCTGTTTGAAAAAGTGCAACGATCCTGGTACGGACTGGATCGc                           <  2:4986271/92‑1 (MQ=255)

                                                               gTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCAt                         >  2:1051107/1‑93 (MQ=255)

                                                               gTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCAt                         >  1:747586/1‑93 (MQ=255)

                                                               gTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCAt                         >  2:7749669/1‑93 (MQ=255)

                                                                tttGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCAtt                        >  2:7861837/1‑93 (MQ=255)

                                                                tttGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCAtt                        >  1:1919872/1‑93 (MQ=255)

                                                                  tGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGt                      <  1:1074310/93‑1 (MQ=255)

                                                                    aGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCa                    >  1:5202944/1‑93 (MQ=255)

                                                                     gcgcGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAg                   <  2:2680993/93‑1 (MQ=255)

                                                                     gcgcGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAg                   >  2:6741305/1‑93 (MQ=255)

                                                                         gTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGcgcc               >  2:7066794/1‑93 (MQ=255)

                                                                         gTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGcgcc               >  1:499309/1‑93 (MQ=255)

                                                                         gTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGcgcc               >  1:4833902/1‑93 (MQ=255)

                                                                          tACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGcgccg              >  2:3280231/1‑93 (MQ=255)

                                                                           tcgggggCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGcgccgc             <  1:7446127/89‑1 (MQ=255)

                                                                            cGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCa            <  2:4876111/93‑1 (MQ=255)

                                                                             gTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAg           >  1:958781/1‑93 (MQ=255)

                                                                              tGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGa          >  1:2038306/1‑93 (MQ=255)

                                                                              tGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGa          <  1:1041203/93‑1 (MQ=255)

                                                                               gggCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGaa         <  2:1970073/93‑1 (MQ=255)

                                                                                ggCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGaaa        >  1:40595/1‑93 (MQ=255)

                                                                                ggCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGaaa        <  1:4714254/93‑1 (MQ=255)

                                                                                ggCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGaaa        <  2:297370/93‑1 (MQ=255)

                                                                                 gCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAg       <  1:6661838/93‑1 (MQ=255)

                                                                                  cGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGt      <  2:1695206/93‑1 (MQ=255)

                                                                                   gATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTc     <  1:6863930/93‑1 (MQ=255)

                                                                                    aTTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCt    >  1:5985651/1‑93 (MQ=255)

                                                                                     ttACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGGGCCGCAGAAAGCCTg   >  1:5643320/1‑93 (MQ=255)

                                                                                      tACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTgc  >  2:6896048/1‑93 (MQ=255)

                                                                                      tACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTgc  <  1:7899362/93‑1 (MQ=255)

CACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTGC  >  NC_000913/2977129‑2977307

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`