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	85.0%	516.6 / 75.7	227	R431R (CGT→CGC)	lhr	Lhr
Reads supporting (aligned to +/- strand): ref base T (17/17); new base C (104/89); total (121/106)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 227 total aligned reads displayed.

AACGCCGCTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCT  >  NC_000913/1728313‑1728492

aaCGCCGCTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTAc                                                                                         <  2:1324888/93‑1 (MQ=255)

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

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

 aCGCCGCTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCAcc                                                                                        <  2:4452941/93‑1 (MQ=255)

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

   gccgcTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCg                                                                                      >  2:5106710/1‑93 (MQ=255)

   gccgcTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCg                                                                                      >  1:1863899/1‑93 (MQ=255)

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

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

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

       cTTTCTGCTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGCCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGt                                                                                  >  1:946908/1‑93 (MQ=255)

         ttCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGtta                                                                                >  2:5406248/1‑91 (MQ=255)

         ttCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGa                                                                                <  2:5994037/93‑1 (MQ=255)

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

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

           cTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGAtt                                                                              <  2:3875556/93‑1 (MQ=255)

           cTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGAtt                                                                              >  1:1082759/1‑93 (MQ=255)

            tGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGAttt                                                                             >  1:2235423/1‑93 (MQ=255)

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

                 ccAGTGGGTTACAACGCATTGGTCGCGCCGGTCATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCg                                                                        >  2:5438814/1‑93 (MQ=255)

                  cAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGa                                                                       >  1:3491595/1‑93 (MQ=255)

                  cAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGa                                                                       <  1:4235189/93‑1 (MQ=255)

                  cAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGa                                                                       <  1:3973819/93‑1 (MQ=255)

                   aGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGAt                                                                      <  2:6210995/93‑1 (MQ=255)

                     tGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTc                                                                    <  1:3003878/93‑1 (MQ=255)

                     tGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTc                                                                    >  2:1760452/1‑93 (MQ=255)

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

                        gTTACAACACATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGc                                                                 >  2:4858073/1‑93 (MQ=255)

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

                           aCAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGt                                                              >  1:6984403/1‑93 (MQ=255)

                             aaCGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCa                                                            >  1:222896/1‑93 (MQ=255)

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

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

                                 cATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGt                                                        <  1:5055332/93‑1 (MQ=255)

                                  aTTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTa                                                       >  2:2693092/1‑93 (MQ=255)

                                  aTTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTa                                                       <  2:6796458/93‑1 (MQ=255)

                                  aTTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTa                                                       >  1:5743140/1‑93 (MQ=255)

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

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

                                    tGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTaga                                                     <  1:4504866/93‑1 (MQ=255)

                                      gTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAgt                                                   <  1:5859068/93‑1 (MQ=255)

                                      gTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAgt                                                   <  2:4270270/93‑1 (MQ=255)

                                       tCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAgtg                                                  <  2:7460576/93‑1 (MQ=255)

                                        cgcgCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAgtgt                                                 <  1:3224957/93‑1 (MQ=255)

                                         gcgcCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCTGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTa                                                <  2:3642590/93‑1 (MQ=255)

                                         gcgcCGGACATCAGGTTGGCGGCGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTa                                                >  2:422517/1‑93 (MQ=255)

                                           gcCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATg                                              <  2:3949893/93‑1 (MQ=255)

                                           gcCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATg                                              <  1:3423326/93‑1 (MQ=255)

                                           gcCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTT‑CCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGt                                             <  1:3742167/93‑1 (MQ=255)

                                             cGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGtt                                            >  2:6870470/1‑93 (MQ=255)

                                             cGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGtt                                            >  2:5084864/1‑93 (MQ=255)

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

                                                aCATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgc                                         >  1:2479600/1‑93 (MQ=255)

                                                 cATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCGCCCGGCATGATTTACTCGATTCCGCAGGCATTGTAGAGTGTATGTTCgca                                        <  1:3626993/93‑1 (MQ=255)

                                                 cATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgca                                        <  2:1880009/93‑1 (MQ=255)

                                                 cATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgca                                        >  1:5229958/1‑93 (MQ=255)

                                                  aTCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcag                                       >  1:1138273/1‑93 (MQ=255)

                                                  aTCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcag                                       >  2:6659662/1‑93 (MQ=255)

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

                                                     aGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcaggca                                    <  1:1002312/93‑1 (MQ=255)

                                                     aGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcaggca                                    <  2:1072512/93‑1 (MQ=255)

                                                     aGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcaggca                                    >  2:2707888/1‑93 (MQ=255)

                                                      ggTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCgcaggcag                                   <  1:2085847/93‑1 (MQ=255)

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

                                                       gTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGcaggcagg                                  >  2:4395436/1‑93 (MQ=255)

                                                        ttGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCaggcaggc                                 >  2:4035823/1‑93 (MQ=255)

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

                                                          nncGGTGTATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTg                               <  2:2542261/91‑1 (MQ=255)

                                                              gtgtATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGaaa                           >  2:6618712/1‑93 (MQ=255)

                                                               tgtATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGaaaa                          >  1:1642088/1‑93 (MQ=255)

                                                               tgtATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGaaaa                          <  1:6451650/93‑1 (MQ=255)

                                                                 tATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAAcc                        >  2:2178856/1‑93 (MQ=255)

                                                                 tATCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAAcc                        <  1:916626/93‑1 (MQ=255)

                                                                  aTCTAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCt                       <  2:5802405/93‑1 (MQ=255)

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

                                                                     tAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGac                    <  2:27272/93‑1 (MQ=255)

                                                                     tAAAGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGac                    <  2:5032577/93‑1 (MQ=255)

                                                                       aaGGGCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGacac                  >  2:2963121/1‑93 (MQ=255)

                                                                        aGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGAcacc                 >  2:757642/1‑93 (MQ=255)

                                                                         gggCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGAcacca                >  1:4960625/1‑93 (MQ=255)

                                                                         gggCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGa                                                                       <  2:1480027/38‑1 (MQ=255)

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

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

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

                                                                          ggCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTGGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGAcaccac               <  1:2888155/93‑1 (MQ=255)

                                                                          ggCTGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGAcaccac               >  2:1391958/1‑93 (MQ=255)

                                                                          agCTGTTTTTCCCCCGCGCCCGGCGTGATTTAGTAGATTCCGCAGTCATTGTAGAGTGTAAGTTCGCAGGCAGGCTGGCAAACCCGGcaccac               >  1:7173046/2‑93 (MQ=255)

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

                                                                             tGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGc            >  2:1250355/1‑93 (MQ=255)

                                                                             tGTTTTTCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGc            <  1:2963065/93‑1 (MQ=255)

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

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

                                                                               tttttCCCCCGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCAt          <  1:1777069/93‑1 (MQ=255)

                                                                                 tttCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATaa        <  1:5419600/93‑1 (MQ=255)

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

                                                                                     ccccGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCct    >  1:1687485/1‑93 (MQ=255)

                                                                                     ccccGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCct    >  1:67601/1‑93 (MQ=255)

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

                                                                                      cccGCACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCCCATAATCctc   >  1:2535879/1‑93 (MQ=255)

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

AACGCCGCTTTCTGTTGCCAGTGGGTTACAACGCATTGGTCGCGCCGGACATCAGGTTGGCGGTGTATCTAAAGGGCTGTTTTTCCCCCGTACCCGGCGTGATTTAGTCGATTCCGCAGTCATTGTAGAGTGTATGTTCGCAGGCAGGCTGGAAAACCTGACACCACCGCATAATCCTCT  >  NC_000913/1728313‑1728492

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