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	99.3%	833.9 / ‑5.3	272	L288P (CTG→CCG)	pdxA	PdxA
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (139/130); total (140/132)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 283 total aligned reads displayed.

AGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCAGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCACGGCGTCG  >  NC_000913/52464‑52643

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

    ccGTAATAAACCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCa                                                                                     <  2:6964617/93‑1 (MQ=255)

    ccGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCa                                                                                     >  2:6795956/1‑93 (MQ=255)

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

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

      gTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACCGATGTGCGAATAAAGGGCGGGCCCAGc                                                                                   >  1:4394783/1‑93 (MQ=255)

       taataaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCg                                                                                  >  2:2377732/1‑93 (MQ=255)

       taataaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCg                                                                                  >  2:3856555/1‑93 (MQ=255)

         ataaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTa                                                                                >  2:4819600/1‑93 (MQ=255)

          taaAGCTGCCGACATCGGCTTTGCCACGCCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTaa                                                                               <  1:4208791/93‑1 (MQ=255)

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

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

                tGCNGACATCGGCTTTGCCACGTCNNGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCAGGCCCAGCGTAATGNTca                                                                         <  2:7547643/93‑1 (MQ=255)

                tGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTca                                                                         <  2:3606540/93‑1 (MQ=255)

                 gCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTcac                                                                        >  1:5658425/1‑93 (MQ=255)

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

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

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

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

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

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

                         cGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGcc                                                                <  1:1222445/93‑1 (MQ=255)

                          ggCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGccc                                                               >  1:2774649/1‑93 (MQ=255)

                          ggCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGccc                                                               >  1:1636404/1‑93 (MQ=255)

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

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

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

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

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

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

                             tttGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGaa                                                            <  2:1581387/93‑1 (MQ=255)

                              ttGCCACGTCCCGCCAGGTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAAAGGTCACACCGCGCCCGAAg                                                           >  2:7039993/1‑93 (MQ=255)

                                  cACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTTATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCt                                                       >  2:5633271/1‑93 (MQ=255)

                                  cACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCt                                                       >  1:1998081/1‑93 (MQ=255)

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

                                  cACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCt                                                       >  1:5436619/1‑93 (MQ=255)

                                   ccGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTg                                                      <  2:620785/92‑1 (MQ=255)

                                   aCGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTg                                                      >  2:2524614/1‑93 (MQ=255)

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

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

                                      tCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTa                                                   >  1:5960590/1‑93 (MQ=255)

                                      tCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTa                                                   <  2:1092145/93‑1 (MQ=255)

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

                                        ccGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTAtt                                                 >  2:4422065/1‑93 (MQ=255)

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

                                          gCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTAtttt                                               <  2:3893982/93‑1 (MQ=255)

                                          gCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTAtttt                                               <  1:2720592/93‑1 (MQ=255)

                                             aGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGc                                            <  1:7500077/93‑1 (MQ=255)

                                               ttCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCAc                                          <  1:5194422/93‑1 (MQ=255)

                                                 cAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACgg                                        >  1:493366/1‑93 (MQ=255)

                                                  aaGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACggg                                       <  1:3360223/93‑1 (MQ=255)

                                                  aaGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACggg                                       <  2:1734525/93‑1 (MQ=255)

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

                                                    gcgcGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGaa                                     <  2:3738927/93‑1 (MQ=255)

                                                    gcgcGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGaa                                     <  2:5479091/93‑1 (MQ=255)

                                                    gcgcGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGaa                                     <  1:1475069/93‑1 (MQ=255)

                                                    gcgcGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGaa                                     <  1:933778/93‑1 (MQ=255)

                                                    gcgcGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGaa                                     <  1:2364316/93‑1 (MQ=255)

                                                    gcgcGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGaa                                     <  1:2870224/93‑1 (MQ=255)

                                                    gcgcGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGCATTTTAGCACGGGaa                                     <  2:7859313/93‑1 (MQ=255)

                                                      gcgGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGa                                   <  2:7890947/93‑1 (MQ=255)

                                                      gcgGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGa                                   >  1:378240/1‑93 (MQ=255)

                                                      gcgGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGa                                   >  1:3040953/1‑93 (MQ=255)

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

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

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

                                                            ccGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGa                             <  2:1201978/93‑1 (MQ=255)

                                                             cGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGAt                            <  2:2028986/93‑1 (MQ=255)

                                                               tGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCg                          >  1:3886091/1‑93 (MQ=255)

                                                               tGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCg                          >  1:7626673/1‑93 (MQ=255)

                                                               tGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCg                          >  2:6421157/1‑93 (MQ=255)

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

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

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

                                                                  tCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTgg                       >  1:3680335/1‑93 (MQ=255)

                                                                   ccACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGt                      >  2:3209363/1‑93 (MQ=255)

                                                                     aCTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTAc                    <  1:7139559/93‑1 (MQ=255)

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

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

                                                                      cTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACa                   >  2:7637405/1‑93 (MQ=255)

                                                                       tGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACac                  >  1:6356338/1‑92 (MQ=255)

                                                                        gATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATc                 >  2:182681/1‑93 (MQ=255)

                                                                        gATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATc                 <  1:6973477/93‑1 (MQ=255)

                                                                        gATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATc                 >  1:1402028/1‑93 (MQ=255)

                                                                        gATGTGCGAATAAAGGGCGGCCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATc                 <  2:5851888/93‑1 (MQ=255)

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

                                                                            tgCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCa             >  1:4244956/1‑93 (MQ=255)

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

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

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

                                                                                aaTAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAc         >  1:6912194/1‑93 (MQ=255)

                                                                                 aTAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCCAAGCCCTGGGATTTTAGCACGGGAAGACCCTGATCGTGACAAATCGCCAGCAcg        >  2:7358162/1‑93 (MQ=255)

                                                                                  tAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcgg       >  1:6491837/1‑93 (MQ=255)

                                                                                  tAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcgg       <  2:6672907/93‑1 (MQ=255)

                                                                                  tAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcgg       >  1:7437791/1‑93 (MQ=255)

                                                                                   aaaGTGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggc      <  2:1047305/93‑1 (MQ=255)

                                                                                     aGGGCGGGCCCCGCGTAATGTTTACACCGCGCCCGAAGCCCCGGGATTTTAGCACGGGAAGACCCCGATCGTGGTACATCGCCCGCAcggcgt    >  2:2961803/1‑93 (MQ=255)

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

                                                                                      gggcgggcCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgtc   <  1:1263320/93‑1 (MQ=255)

                                                                                       ggCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgtcg  >  1:4115398/1‑93 (MQ=255)

AGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCAGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCACGGCGTCG  >  NC_000913/52464‑52643

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