BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	52,554	A→G	L288P (CTG→CCG)	pdxA ←	PdxA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	52,554	0	A	G	98.2%	345.7 / ‑2.3	114	L288P (CTG→CCG)	pdxA	PdxA
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (54/58); total (55/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-01

Only 100 of 126 total aligned reads displayed.

AGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCAGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCACGGCGTCG  >  NC_000913/52464‑52643

aGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCggg                                                                                         <  2:2317475/93‑1 (MQ=255)

 gcgcCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAgggcgggc                                                                                        >  2:2507672/1‑93 (MQ=255)

 gcgcCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAgggcgggc                                                                                        >  2:557746/1‑93 (MQ=255)

   gcCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGccc                                                                                      >  2:4251471/1‑93 (MQ=255)

   gcCGCAATAAAACTACCGACATCGCCATTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGCCCACTGACGTGCGACTAAAGTGCGGGccc                                                                                      >  2:376142/1‑93 (MQ=255)

    ccGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCAGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGTGCGTGCCCa                                                                                     >  2:4201492/1‑93 (MQ=255)

     cGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAg                                                                                    >  2:2259777/1‑93 (MQ=255)

     cGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAg                                                                                    <  2:3266674/93‑1 (MQ=255)

      gTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGc                                                                                   >  1:2663094/1‑93 (MQ=255)

      gTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGc                                                                                   >  2:3366644/1‑93 (MQ=255)

      gTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGc                                                                                   >  2:764729/1‑93 (MQ=255)

      gTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGc                                                                                   <  1:540722/93‑1 (MQ=255)

        aataaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGt                                                                                 >  1:1930017/1‑93 (MQ=255)

        aataaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGCCGGGCCCAGCGt                                                                                 >  2:1577548/1‑93 (MQ=255)

          taaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTaa                                                                               <  1:1608518/93‑1 (MQ=255)

          taaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTaa                                                                               >  2:1453033/1‑93 (MQ=255)

           aaagCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGGCCACTGAGGTGCGAATAAAGGGCGGGCCCAGCGTAAt                                                                              <  2:4061158/93‑1 (MQ=255)

           aaaaCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAAt                                                                              >  1:2395588/1‑93 (MQ=255)

            aaaCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATg                                                                             >  2:2898415/1‑93 (MQ=255)

             aaCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAACGTAATGt                                                                            >  1:3796466/1‑93 (MQ=255)

               cTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTc                                                                          >  2:3112490/1‑93 (MQ=255)

                tGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCAGGCCCAGCGTAATGTTca                                                                         <  1:4510479/93‑1 (MQ=255)

                  ccGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTcaca                                                                       <  1:1822792/93‑1 (MQ=255)

                   cGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCacac                                                                      <  2:1938622/93‑1 (MQ=255)

                   cGACATCGGCTTTGCCACGTCCCGCCAGTGCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCacac                                                                      <  2:1938623/93‑1 (MQ=255)

                    gACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACAcc                                                                     <  1:2041455/93‑1 (MQ=255)

                    gACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACAcc                                                                     <  2:436604/93‑1 (MQ=255)

                      cATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACcgc                                                                   <  1:3252811/93‑1 (MQ=255)

                      cATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACcgc                                                                   >  1:631114/1‑93 (MQ=255)

                       aTCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACcgcg                                                                  >  1:3327309/1‑93 (MQ=255)

                       aTCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACcgcg                                                                  >  1:3811363/1‑93 (MQ=255)

                        tCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCgcgc                                                                 >  1:3959303/1‑93 (MQ=255)

                        tCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGACGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCgcgc                                                                 >  1:3178799/1‑93 (MQ=255)

                         cGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGCGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGcc                                                                <  2:2998947/93‑1 (MQ=255)

                          ggCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGccc                                                               <  1:4282796/93‑1 (MQ=255)

                           gCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCg                                                              >  2:2070694/1‑93 (MQ=255)

                           gCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCg                                                              <  2:3164392/93‑1 (MQ=255)

                           gCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCg                                                              <  2:3164393/93‑1 (MQ=255)

                            cTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGa                                                             <  2:1344073/93‑1 (MQ=255)

                              ttGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAg                                                           <  1:3607049/93‑1 (MQ=255)

                              ttGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAg                                                           <  1:1100156/93‑1 (MQ=255)

                                gCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGcc                                                         >  1:1033735/1‑93 (MQ=255)

                                gCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGcc                                                         >  2:3964981/1‑93 (MQ=255)

                                 ccACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCATTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGccc                                                        <  1:3284656/93‑1 (MQ=255)

                                 ccACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGccc                                                        <  2:1966706/93‑1 (MQ=255)

                                  cACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCt                                                       <  1:1393820/93‑1 (MQ=255)

                                    cGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTgg                                                     <  2:2194877/93‑1 (MQ=255)

                                     gTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGt                                                    >  2:519630/1‑93 (MQ=255)

                                     gTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGt                                                    <  2:1918477/93‑1 (MQ=255)

                                      tCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCTAAGCCCTGGtt                                                   >  1:2875269/1‑92 (MQ=255)

                                       cccGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTAt                                                  >  1:1418187/1‑93 (MQ=255)

                                        ccGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCAGACCGCGCCCGAAGCCCTGGTAtt                                                 >  2:2591218/1‑93 (MQ=255)

                                         cGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTAttt                                                >  1:3257275/1‑93 (MQ=255)

                                         cGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTAttt                                                <  2:2454203/93‑1 (MQ=255)

                                           ccAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTa                                              >  2:934088/1‑93 (MQ=255)

                                           ccAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTa                                              >  2:1216319/1‑93 (MQ=255)

                                             aGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGc                                            >  1:1136604/1‑93 (MQ=255)

                                             aGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGc                                            >  1:401699/1‑93 (MQ=255)

                                              gTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCa                                           >  1:862083/1‑93 (MQ=255)

                                              gTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCa                                           >  2:2511107/1‑93 (MQ=255)

                                               ttCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCAc                                          >  2:382646/1‑93 (MQ=255)

                                                   aGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGa                                      <  2:1649590/93‑1 (MQ=255)

                                                    gcgcGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGaa                                     >  1:564217/1‑93 (MQ=255)

                                                      gcgGTGCCGTGGTCCACTGATGTGCGAATCAAGCGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGa                                   <  1:2575097/93‑1 (MQ=255)

                                                       cgGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGAc                                  <  2:262201/93‑1 (MQ=255)

                                                         gTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGAccc                                >  2:2358095/1‑93 (MQ=255)

                                                          tGCCGTGGTCCACTGCTGTGCGCCTAGAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGACGCCCTGGTATTTTAGCACGGGAAGACCCt                               <  1:4225534/93‑1 (MQ=255)

                                                           gCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTg                              <  2:1840403/93‑1 (MQ=255)

                                                            ccGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGCATTTTAGCACGGGAAGACCCTGa                             >  2:4331032/1‑93 (MQ=255)

                                                             cGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGAt                            >  2:453490/1‑93 (MQ=255)

                                                                ggTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGt                         >  2:3219663/1‑93 (MQ=255)

                                                                ggTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGt                         <  1:2358522/93‑1 (MQ=255)

                                                                 gTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTg                        <  1:2977808/93‑1 (MQ=255)

                                                                 gTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTg                        <  2:1057275/93‑1 (MQ=255)

                                                                 gTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTg                        <  2:4181780/93‑1 (MQ=255)

                                                                     aCTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTAc                    >  1:1127130/1‑93 (MQ=255)

                                                                       tGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACCCCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACAt                  <  2:3452993/93‑1 (MQ=255)

                                                                       tGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACAt                  >  2:2309285/1‑93 (MQ=255)

                                                                          tgtgCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGc               >  2:2724638/1‑93 (MQ=255)

                                                                            tgCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCa             <  1:3638874/93‑1 (MQ=255)

                                                                            tgCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCa             <  1:2943598/93‑1 (MQ=255)

                                                                            tgCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCa             <  1:1454780/93‑1 (MQ=255)

                                                                             gCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAg            <  1:4644999/93‑1 (MQ=255)

                                                                              cGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGc           <  2:881747/93‑1 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          <  2:4408551/93‑1 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          >  1:3441179/1‑93 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          <  2:1353923/93‑1 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          <  2:4616341/93‑1 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          <  2:3850811/93‑1 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          <  2:3850810/93‑1 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          >  1:2035963/1‑93 (MQ=255)

                                                                                aaTAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAc         <  2:3269524/93‑1 (MQ=255)

                                                                                 aTAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcg        >  2:4457831/1‑93 (MQ=255)

                                                                                  tAAAGGGCAGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGACGACCCTGATCGTGGTACATCGCCAGCAcgg       >  1:1567331/1‑93 (MQ=255)

                                                                                  aaaaGTGCGGGCCCAGCGTAATGTTCACACCGCTCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcgg       <  1:1457065/92‑1 (MQ=255)

                                                                                   aaaGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggc      >  2:356274/1‑93 (MQ=255)

                                                                                     ggggcgggcCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTGGCGCGGGAAGACCCTGATCGTGGTACATCACCAGCAcggcgt    <  2:4223820/92‑1 (MQ=255)

                                                                                     aGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgt    >  2:2077934/1‑93 (MQ=255)

                                                                                     aGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGNGGTACATCGCCAGCAcggcgt    >  1:651173/1‑93 (MQ=255)

                                                                                      gggcgggcCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgtc   >  2:259174/1‑93 (MQ=255)

AGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCAGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCACGGCGTCG  >  NC_000913/52464‑52643

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