BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,728,403	T→C	R431R (CGT→CGC)	lhr →	Lhr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,728,403	0	T	C	99.1%	731.0 / ‑4.9	221	R431R (CGT→CGC)	lhr	Lhr
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (112/107); total (112/109)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 229 total aligned reads displayed.

AACGCCGCTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCT  >  NC_000913/1728313‑1728492

aaCGCCGCTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCAc                                                                                         <  2:4925263/93‑1 (MQ=255)

 aCGCCGCTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCAcc                                                                                        >  2:3105731/1‑93 (MQ=255)

  cgccgcTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCAccc                                                                                       >  2:4981565/1‑93 (MQ=255)

  cgccgcTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCAccc                                                                                       <  1:2492145/93‑1 (MQ=255)

    ccgcTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCgg                                                                                     >  2:2493366/1‑93 (MQ=255)

    ccgcTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCgg                                                                                     <  2:678932/93‑1 (MQ=255)

     cgcTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGc                                                                                    <  2:2744848/93‑1 (MQ=255)

     cgcTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGc                                                                                    >  1:7796800/1‑93 (MQ=255)

      gcTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTTCCCCGCACCCGGCg                                                                                   >  1:6465342/1‑93 (MQ=255)

      gcTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCg                                                                                   <  1:8644934/93‑1 (MQ=255)

      gcTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCg                                                                                   >  1:5104664/1‑93 (MQ=255)

      gcTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCg                                                                                   >  1:1827046/1‑93 (MQ=255)

       cTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGt                                                                                  <  1:6850553/93‑1 (MQ=255)

          tCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGAt                                                                               >  2:5118157/1‑93 (MQ=255)

          tCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGAt                                                                               <  2:632976/93‑1 (MQ=255)

             gTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTa                                                                            <  2:2804809/93‑1 (MQ=255)

               tGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGt                                                                          <  1:5512256/93‑1 (MQ=255)

                gCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTc                                                                         >  1:6951731/1‑93 (MQ=255)

                gCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTc                                                                         >  2:4034328/1‑93 (MQ=255)

                gCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTc                                                                         <  2:7635331/93‑1 (MQ=255)

                   aGTGTGTTACAACGCATTGGTCTCGCCGGAAATCAGGTTGGCGGTGTATCTAAAGGTCTGTTTTTCACCCGTACCCGGCGTGATTTAGTCGAt                                                                      <  1:6008343/93‑1 (MQ=255)

                    gTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGAtt                                                                     <  1:2896999/93‑1 (MQ=255)

                    gTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGAtt                                                                     >  2:3741185/1‑93 (MQ=255)

                     tGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTc                                                                    <  1:7010042/93‑1 (MQ=255)

                     tGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTc                                                                    <  2:4419721/93‑1 (MQ=255)

                      gggTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTcc                                                                   <  2:2859484/93‑1 (MQ=255)

                       ggTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCg                                                                  >  2:7183584/1‑93 (MQ=255)

                       ggTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCg                                                                  >  1:7576411/1‑93 (MQ=255)

                       ggTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCg                                                                  >  2:1980303/1‑93 (MQ=255)

                        gTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGc                                                                 <  2:525030/93‑1 (MQ=255)

                         ttACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCa                                                                <  2:8669994/93‑1 (MQ=255)

                         ttACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCa                                                                >  1:2745141/1‑93 (MQ=255)

                         taaCAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCa                                                                >  2:6102306/3‑93 (MQ=255)

                          tACAGCGCATTGGGCCCGCCGGACATCAGGTTGGCGGTGTATCTAAGGGGCTGTTTTGCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAg                                                               <  2:6744864/93‑1 (MQ=255)

                            cAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTc                                                             >  2:7950101/1‑93 (MQ=255)

                            cAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTc                                                             >  2:6598211/1‑93 (MQ=255)

                            cAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTc                                                             >  1:766280/1‑93 (MQ=255)

                              aCGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCAt                                                           <  2:4830311/93‑1 (MQ=255)

                              aCGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCAt                                                           >  2:1425882/1‑93 (MQ=255)

                              aCGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTCCCCCCCGCACGCGGCGTGATTTAGTCGATTCCGCAGTCAt                                                           <  1:1346138/93‑1 (MQ=255)

                               cGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCAtt                                                          <  1:6655784/93‑1 (MQ=255)

                               cGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCAtt                                                          <  2:1094177/93‑1 (MQ=255)

                                 cATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGt                                                        >  1:7428268/1‑93 (MQ=255)

                                   ttGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTag                                                      <  2:8378400/93‑1 (MQ=255)

                                   ttGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTag                                                      <  2:5400331/93‑1 (MQ=255)

                                   ttGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTag                                                      >  2:574570/1‑93 (MQ=255)

                                   ttGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTag                                                      >  1:3670090/1‑93 (MQ=255)

                                     ggTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTagag                                                    >  1:801211/1‑93 (MQ=255)

                                     ggTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTagag                                                    <  2:815941/93‑1 (MQ=255)

                                     ggTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTagag                                                    <  1:4412095/93‑1 (MQ=255)

                                        cgcgCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAgtgt                                                 <  1:2122185/93‑1 (MQ=255)

                                        cgcgCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAgtgt                                                 >  1:7118014/1‑93 (MQ=255)

                                        cgcgCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTACCCCGTACTCGCCGTGCTTTAGTCGATTCCGCAGTCCTTGCagagaga                                                 >  2:5584683/1‑90 (MQ=255)

                                         gcgcCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTa                                                <  1:2791624/93‑1 (MQ=255)

                                            gcGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGt                                             <  1:6745949/92‑1 (MQ=255)

                                            ccGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGt                                             <  2:7533277/93‑1 (MQ=255)

                                              ggACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTc                                           <  1:1648727/93‑1 (MQ=255)

                                              ggACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTc                                           >  1:4540863/1‑93 (MQ=255)

                                              ggACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTc                                           >  1:2600049/1‑93 (MQ=255)

                                              ggACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGGGTATGttg                                           >  1:1205186/1‑92 (MQ=255)

                                               gACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCg                                          >  2:5851430/1‑93 (MQ=255)

                                               gACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCg                                          <  1:1727459/93‑1 (MQ=255)

                                                  aTCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcag                                       <  2:4890019/93‑1 (MQ=255)

                                                   tCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTTCGCAGGCATTGTAGAGTGTATGTTCgcagg                                      >  2:8080697/1‑93 (MQ=255)

                                                    cAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcaggc                                     <  1:1985031/93‑1 (MQ=255)

                                                    cAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcaggc                                     <  1:1925292/93‑1 (MQ=255)

                                                     aGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcaggca                                    <  1:2730312/93‑1 (MQ=255)

                                                      ggTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTTCGCAGTCATTGTAGAGTGTATGTTCgcaggcag                                   >  1:5155796/1‑93 (MQ=255)

                                                      ggTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcaggcag                                   >  1:7244133/1‑93 (MQ=255)

                                                      ggTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcaggcag                                   >  1:4797502/1‑93 (MQ=255)

                                                      ggTTGGCGGTGTATCTAAAGGGCTGTTTTCCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcaggcag                                   <  1:974018/93‑1 (MQ=255)

                                                        ttGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCaggcaggc                                 <  2:7325076/93‑1 (MQ=255)

                                                        ttGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCaggcaggc                                 <  1:5746500/93‑1 (MQ=255)

                                                         tGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCt                                >  2:7267558/1‑93 (MQ=255)

                                                          ggCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTg                               <  1:3368175/93‑1 (MQ=255)

                                                           gCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTgg                              <  1:5445512/93‑1 (MQ=255)

                                                            cGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGa                             <  1:8663970/93‑1 (MQ=255)

                                                            cGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGa                             >  1:1076139/1‑93 (MQ=255)

                                                            cGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGa                             <  2:1825005/93‑1 (MQ=255)

                                                             ggTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGaa                            >  1:8380246/1‑93 (MQ=255)

                                                              gtgtATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGaaa                           <  2:7359381/93‑1 (MQ=255)

                                                                    cTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGa                     >  2:3547381/1‑93 (MQ=255)

                                                                    cTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGa                     <  2:4413363/93‑1 (MQ=255)

                                                                     tAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGac                    >  1:8510422/1‑93 (MQ=255)

                                                                     tAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGac                    >  2:6717609/1‑93 (MQ=255)

                                                                       aaGGGCTGTTTTTCCCCCGCTCCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGacac                  <  2:6410702/93‑1 (MQ=255)

                                                                        aGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGAcacc                 <  1:1550509/93‑1 (MQ=255)

                                                                         gggCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGAcacca                >  2:295791/1‑93 (MQ=255)

                                                                           gCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACaccacc              <  1:8751846/93‑1 (MQ=255)

                                                                             tGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTTACACCACCGc            >  2:4289214/1‑93 (MQ=255)

                                                                              gTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCa           >  1:388582/1‑93 (MQ=255)

                                                                              gTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGAAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCa           >  1:4122701/1‑93 (MQ=255)

                                                                                  ttCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAAt       >  1:1465694/1‑93 (MQ=255)

                                                                                   tCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATc      >  2:7091235/1‑93 (MQ=255)

                                                                                   tCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATc      <  1:4230696/93‑1 (MQ=255)

                                                                                   tCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATc      <  2:330766/93‑1 (MQ=255)

                                                                                      cccGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCctc   >  2:1921609/1‑93 (MQ=255)

                                                                                      cccGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCctc   <  2:8800621/93‑1 (MQ=255)

                                                                                       ccGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCctct  >  1:2334485/1‑93 (MQ=255)

                                                                                       ccGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCctct  <  2:8633067/93‑1 (MQ=255)

AACGCCGCTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCT  >  NC_000913/1728313‑1728492

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