BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	361,946	A→G	L153P (CTG→CCG)	lacY ←	galactoside permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	361,946	0	A	G	97.2%	531.6 / 1.0	176	L153P (CTG→CCG)	lacY	galactoside permease
Reads supporting (aligned to +/- strand): ref base A (2/3); new base G (78/93); total (80/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 185 total aligned reads displayed.

AAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCT  >  NC_000913/361856‑362035

aaaaCGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGgc                                                                                         <  1:48480/93‑1 (MQ=255)

  aaCGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGgcgc                                                                                       >  2:53309/1‑93 (MQ=255)

   aCGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGcc                                                                                      <  2:3888954/93‑1 (MQ=255)

   aCGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGACGACAATAGAGGCACAGGGCGcc                                                                                      >  1:725859/1‑93 (MQ=255)

      gCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAg                                                                                   >  2:95432/1‑93 (MQ=255)

       cGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGc                                                                                  <  1:4376017/93‑1 (MQ=255)

        gAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGcc                                                                                 <  1:4720518/93‑1 (MQ=255)

           gATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAac                                                                              >  2:2849140/1‑93 (MQ=255)

           gATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAac                                                                              >  1:263483/1‑93 (MQ=255)

            aTGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAaca                                                                             <  2:3232988/93‑1 (MQ=255)

             tGATTGCACTGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAacac                                                                            <  2:2602385/93‑1 (MQ=255)

              gAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAAcaca                                                                           >  1:2278595/1‑93 (MQ=255)

              gAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAAcaca                                                                           >  1:4116395/1‑93 (MQ=255)

                    acaGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCaaa                                                                     <  1:664194/93‑1 (MQ=255)

                      aGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACCCAGCCAAACa                                                                   <  2:1537239/93‑1 (MQ=255)

                      aGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACa                                                                   >  2:1957613/1‑93 (MQ=255)

                        ccAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATc                                                                 >  1:2886036/1‑93 (MQ=255)

                        ccAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATc                                                                 <  1:4650395/93‑1 (MQ=255)

                        ccAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATc                                                                 <  2:2229982/93‑1 (MQ=255)

                         cAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATcc                                                                >  2:1054084/1‑93 (MQ=255)

                         cAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATcc                                                                <  2:4470374/93‑1 (MQ=255)

                         cAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATcc                                                                <  1:3416024/93‑1 (MQ=255)

                          agagCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCcg                                                               <  2:4439826/93‑1 (MQ=255)

                            agCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCcgcg                                                             >  1:1894463/1‑93 (MQ=255)

                            agCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCcgcg                                                             >  2:4564761/1‑93 (MQ=255)

                             gCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCcgcgc                                                            <  2:2317552/93‑1 (MQ=255)

                             gCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCcgcgc                                                            <  1:2121144/93‑1 (MQ=255)

                             gCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCcgcgc                                                            >  1:3581135/1‑93 (MQ=255)

                             gCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCcgcgc                                                            >  2:3837316/1‑93 (MQ=255)

                              cccAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCcgcgcg                                                           <  1:440021/93‑1 (MQ=255)

                               ccagccagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCcgcgcgc                                                          >  2:4094212/1‑93 (MQ=255)

                                 agccagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGa                                                        <  1:2545909/93‑1 (MQ=255)

                                 agccagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGa                                                        <  1:236923/93‑1 (MQ=255)

                                  gccagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGAc                                                       >  2:2352125/1‑93 (MQ=255)

                                    cagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCa                                                     <  1:1438169/93‑1 (MQ=255)

                                    cagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCa                                                     >  2:936355/1‑93 (MQ=255)

                                    cagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCa                                                     >  2:1842712/1‑93 (MQ=255)

                                      gAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCaaa                                                   <  1:275774/93‑1 (MQ=255)

                                      gAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCaaa                                                   <  2:2671933/93‑1 (MQ=255)

                                       aaaaCAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAAt                                                  <  2:2750691/93‑1 (MQ=255)

                                       aaaaCAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAAt                                                  <  2:2801394/93‑1 (MQ=255)

                                        aaacaaacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAAtt                                                 <  2:3270696/93‑1 (MQ=255)

                                        aaacaaacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAAtc                                                 >  1:2854482/1‑92 (MQ=255)

                                         aacaaacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTc                                                >  1:4171477/1‑93 (MQ=255)

                                         aacaaacTGATTATTGATGGGGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTc                                                <  1:1342211/93‑1 (MQ=255)

                                          acaaacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCg                                               >  2:370767/1‑93 (MQ=255)

                                             aacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGaaa                                            <  2:872765/93‑1 (MQ=255)

                                             aacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGaaa                                            >  2:1807097/1‑93 (MQ=255)

                                             aacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGaaa                                            <  1:1750267/93‑1 (MQ=255)

                                             aacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGaaa                                            <  2:2769175/93‑1 (MQ=255)

                                             aacTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGaaa                                            >  1:3797571/1‑93 (MQ=255)

                                               cTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCGGCCAACACAGCCAAACATCCGCGCGCGCCCAAATTCGAAAtt                                          <  2:2644681/93‑1 (MQ=255)

                                                tGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTa                                         >  2:2880272/1‑93 (MQ=255)

                                                 gATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTAc                                        >  1:3195613/1‑93 (MQ=255)

                                                 gATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTAc                                        <  1:3747541/93‑1 (MQ=255)

                                                 gATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTAc                                        >  2:3448326/1‑93 (MQ=255)

                                                 gATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTAc                                        <  2:1637624/93‑1 (MQ=255)

                                                 gATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTAc                                        >  2:1850667/1‑93 (MQ=255)

                                                  attattGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACt                                       <  2:3132427/93‑1 (MQ=255)

                                                  attatgGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACt                                       <  1:2909746/93‑1 (MQ=255)

                                                    tattGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGc                                     <  2:2180736/93‑1 (MQ=255)

                                                    tatcGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGc                                     <  2:2835538/93‑1 (MQ=255)

                                                     attGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCg                                    <  2:2085665/93‑1 (MQ=255)

                                                      ttGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGa                                   >  1:1919422/1‑93 (MQ=255)

                                                       tGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGAc                                  <  1:2968819/93‑1 (MQ=255)

                                                       tGATGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGAc                                  <  1:1955742/93‑1 (MQ=255)

                                                        tgtGGTGAACAGGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACg                                 <  2:1330635/91‑1 (MQ=255)

                                                         aTGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACgg                                >  1:3298468/1‑93 (MQ=255)

                                                         aTGGTGAACATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACgc                                >  1:4550514/1‑92 (MQ=255)

                                                               aaCATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACt                          >  1:707872/1‑93 (MQ=255)

                                                               aaCATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACt                          >  2:2126281/1‑93 (MQ=255)

                                                                 cATGATGCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACttt                        >  2:3777274/1‑93 (MQ=255)

                                                                  atgatgCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTtc                       >  2:1159306/1‑93 (MQ=255)

                                                                  atgatgCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTtc                       <  2:3148716/93‑1 (MQ=255)

                                                                   tgatgCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTtct                      >  1:1263619/1‑93 (MQ=255)

                                                                      tgCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAGTTCGAAATTACTGCGACGGCTGACTTTCTCaa                   >  2:1121771/1‑93 (MQ=255)

                                                                      tgCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCaa                   <  1:547565/93‑1 (MQ=255)

                                                                      tgCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCaa                   >  2:666361/1‑93 (MQ=255)

                                                                       gCCGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACCTTCTCAAt                  <  2:637003/93‑1 (MQ=255)

                                                                        ccGACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATa                 >  1:3074749/1‑93 (MQ=255)

                                                                          gACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATaaa               >  2:286545/1‑93 (MQ=255)

                                                                          gACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATaaa               <  2:2842842/93‑1 (MQ=255)

                                                                          gACAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATaaa               <  2:3522787/93‑1 (MQ=255)

                                                                           aCAATCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAAt              <  1:3468044/93‑1 (MQ=255)

                                                                             aaTCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGc            <  1:1656118/93‑1 (MQ=255)

                                                                             aaTCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTAAATAAATGc            >  2:3378215/1‑93 (MQ=255)

                                                                               tCGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCct          <  2:2038262/93‑1 (MQ=255)

                                                                                cGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCctc         <  2:4294369/93‑1 (MQ=255)

                                                                                cGAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCctc         <  1:4511673/93‑1 (MQ=255)

                                                                                cGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCctc         <  2:4018006/93‑1 (MQ=255)

                                                                                 gAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCctct        <  2:297469/93‑1 (MQ=255)

                                                                                 gAGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCctct        >  1:2536005/1‑93 (MQ=255)

                                                                                  aGGCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTa       <  2:3841383/93‑1 (MQ=255)

                                                                                   ggCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTAc      <  1:2987847/93‑1 (MQ=255)

                                                                                   ggCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTAc      <  1:4429103/93‑1 (MQ=255)

                                                                                   ggCACACGGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTAc      <  2:1077240/93‑1 (MQ=255)

                                                                                     cacacgGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTActg    >  1:4207354/1‑93 (MQ=255)

                                                                                      acacgGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTActgc   >  2:2469767/1‑93 (MQ=255)

                                                                                      acacaGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTActgc   <  1:1590619/93‑1 (MQ=255)

                                                                                       cacgGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTActgct  >  1:2184265/1‑93 (MQ=255)

AAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCT  >  NC_000913/361856‑362035

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