BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,158,281	C→T	L294L (CTG→CTA)	ygiRQ ←	Fe‑S oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,158,281	0	C	T	85.3%	812.9 / 96.3	312	L294L (CTG→CTA)	ygiRQ	Fe‑S oxidoreductase
Reads supporting (aligned to +/- strand): ref base C (24/22); new base T (127/139); total (151/161)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-02

Only 100 of 325 total aligned reads displayed.

GCGCGGGCACAG‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTT  >  NC_000913/3158190‑3158370

gcgcGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgc                                                                                           <  2:7212298/93‑1 (MQ=255)

  gcgGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgcag                                                                                         >  1:3180237/1‑93 (MQ=255)

    gggCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgtagca                                                                                       <  1:6377585/93‑1 (MQ=255)

      gCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAac                                                                                     <  1:624806/93‑1 (MQ=255)

         cagccTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACg                                                                                  <  1:4200687/89‑1 (MQ=255)

           gccTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTa                                                                                <  1:6687598/91‑1 (MQ=255)

           gccTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTa                                                                                >  2:1462632/3‑93 (MQ=255)

            ccTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAg                                                                               >  1:2357758/2‑93 (MQ=255)

            ccTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAg                                                                               <  2:5579216/92‑1 (MQ=255)

             cTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAgg                                                                              <  2:3461206/93‑1 (MQ=255)

               gggTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGtt                                                                            >  1:2262067/1‑93 (MQ=255)

                ggttggttTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGttt                                                                           <  1:866687/93‑1 (MQ=255)

                ggttggttTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGttt                                                                           <  1:1013626/93‑1 (MQ=255)

                ggttggttTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGttt                                                                           >  2:5744242/1‑93 (MQ=255)

                ggttggttTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGttt                                                                           <  1:1959722/93‑1 (MQ=255)

                  ttggttTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGttttt                                                                         <  2:526789/93‑1 (MQ=255)

                  ttggctTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGttttt                                                                         >  2:6997054/1‑93 (MQ=255)

                   tggttTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGtttttt                                                                        >  2:4982025/1‑93 (MQ=255)

                     gttTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTcc                                                                      >  1:5441891/1‑93 (MQ=255)

                          gtggtgCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACgg                                                                 >  2:3294614/1‑93 (MQ=255)

                            ggtgCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGtt                                                               >  2:46585/1‑93 (MQ=255)

                               gCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCg                                                            >  2:468762/1‑93 (MQ=255)

                                cAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCgg                                                           <  1:4877607/93‑1 (MQ=255)

                                cAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCgg                                                           <  2:3883406/93‑1 (MQ=255)

                                cAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCgg                                                           <  1:4400042/93‑1 (MQ=255)

                                cAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCgg                                                           <  2:7461326/93‑1 (MQ=255)

                                 aGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCggg                                                          <  1:7333106/93‑1 (MQ=255)

                                  gAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGgc                                                         <  1:2608967/93‑1 (MQ=255)

                                  gAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGgc                                                         >  2:265926/1‑93 (MQ=255)

                                  gAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGgc                                                         >  2:2186988/1‑93 (MQ=255)

                                       cgcgAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGt                                                    <  2:4312685/93‑1 (MQ=255)

                                       cgcgAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGt                                                    <  1:5892113/93‑1 (MQ=255)

                                       cgcgAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGt                                                    >  1:6923631/1‑93 (MQ=255)

                                       cgcgAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGt                                                    >  2:2407435/1‑93 (MQ=255)

                                       cgcgAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGt                                                    <  1:7234936/93‑1 (MQ=255)

                                        gcgAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTg                                                   >  1:3680640/1‑93 (MQ=255)

                                        gcgAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTg                                                   <  1:7096354/93‑1 (MQ=255)

                                          gAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGc                                                 >  2:6611228/1‑93 (MQ=255)

                                          gAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGc                                                 <  1:6523381/93‑1 (MQ=255)

                                          gAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGc                                                 <  2:5984168/93‑1 (MQ=255)

                                               aTGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCAc                                            <  1:1504432/93‑1 (MQ=255)

                                                tGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACg                                           >  1:5223024/1‑93 (MQ=255)

                                                tGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACg                                           <  2:7335184/93‑1 (MQ=255)

                                                tGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACg                                           <  1:2055353/93‑1 (MQ=255)

                                                tGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACg                                           >  2:5911825/1‑93 (MQ=255)

                                                 gggCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACgg                                          <  1:5992683/93‑1 (MQ=255)

                                                 gggCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACgg                                          >  2:631801/1‑93 (MQ=255)

                                                 gggCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACgg                                          <  2:6112525/93‑1 (MQ=255)

                                                  ggCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGt                                         >  1:4154553/1‑93 (MQ=255)

                                                  ggCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGt                                         <  2:41786/93‑1 (MQ=255)

                                                  ggCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGt                                         <  1:3636403/93‑1 (MQ=255)

                                                   gCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGtt                                        >  1:2856268/1‑93 (MQ=255)

                                                    cGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTa                                       <  2:1059993/93‑1 (MQ=255)

                                                     gTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTAc                                      <  2:3975179/93‑1 (MQ=255)

                                                     gTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTAc                                      <  1:2875207/93‑1 (MQ=255)

                                                      tACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGATTCGGGCGCGGTGGCTGCACGGTTACg                                     <  1:5661056/93‑1 (MQ=255)

                                                       aCAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACgg                                    <  1:1234458/93‑1 (MQ=255)

                                                        cAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGc                                   <  2:3011615/93‑1 (MQ=255)

                                                        cAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGc                                   >  2:4078189/1‑93 (MQ=255)

                                                           cACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCttt                                <  2:2070756/93‑1 (MQ=255)

                                                           cACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCttt                                >  2:4761395/1‑93 (MQ=255)

                                                            aCTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTg                               <  1:34940/93‑1 (MQ=255)

                                                            aCTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTg                               <  1:5318941/93‑1 (MQ=255)

                                                            aCTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTg                               >  2:1258459/1‑93 (MQ=255)

                                                             cTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTgg                              >  1:2554018/1‑93 (MQ=255)

                                                              tttATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGc                             <  2:3611218/93‑1 (MQ=255)

                                                                tATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCtt                           <  2:5337121/93‑1 (MQ=255)

                                                                tATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCtt                           >  1:3167165/1‑93 (MQ=255)

                                                                 aTCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTc                          >  1:3811517/1‑93 (MQ=255)

                                                                 aTCGCCCTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTc                          >  2:561662/1‑93 (MQ=255)

                                                                  tCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTcc                         <  2:2941742/93‑1 (MQ=255)

                                                                     cccTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCAGc                      <  2:2564223/93‑1 (MQ=255)

                                                                      ccTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCt                     >  2:3805003/1‑93 (MQ=255)

                                                                       cTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCtt                    >  2:2349960/1‑93 (MQ=255)

                                                                       cTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCtt                    >  1:7181828/1‑93 (MQ=255)

                                                                       cTTCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCtt                    >  1:1857361/1‑93 (MQ=255)

                                                                       cTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCtt                    >  2:2701730/1‑93 (MQ=255)

                                                                        ttCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCttt                   >  1:7518612/1‑93 (MQ=255)

                                                                        ttCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCttt                   >  2:6214357/1‑93 (MQ=255)

                                                                        ttCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCttt                   >  1:5823506/1‑93 (MQ=255)

                                                                         tCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCtttt                  <  1:5697483/93‑1 (MQ=255)

                                                                         tCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCtttt                  >  1:4696803/1‑93 (MQ=255)

                                                                         tCACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCtttt                  >  1:4465750/1‑93 (MQ=255)

                                                                          cACTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCttttt                 <  1:2098091/93‑1 (MQ=255)

                                                                           aCTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTc                >  1:5952982/1‑93 (MQ=255)

                                                                           aCTTTCTCGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTc                <  2:1667084/93‑1 (MQ=255)

                                                                                ctcGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc           >  1:1619020/1‑93 (MQ=255)

                                                                                ctcGAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc           >  2:1732774/1‑93 (MQ=255)

                                                                                   gAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCAc        >  2:932720/1‑93 (MQ=255)

                                                                                   gAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCAc        >  2:893425/1‑93 (MQ=255)

                                                                                   gAAAGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCAc        >  2:5238357/1‑93 (MQ=255)

                                                                                    aaaGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCAcc       >  1:4565497/1‑93 (MQ=255)

                                                                                    aaaGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCAcc       >  2:511621/1‑93 (MQ=255)

                                                                                    aaaGAAGGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCAcc       >  2:3199382/1‑93 (MQ=255)

                                                                                     aagaagGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCg      >  1:4074328/1‑93 (MQ=255)

                                                                                     aagaagGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCg      >  2:4457683/1‑93 (MQ=255)

                                                                                      agaagGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCgg     <  1:3381597/93‑1 (MQ=255)

                                                                                       gaagGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGt    >  2:2244841/1‑93 (MQ=255)

                                                                                        aagGTAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGtt   >  1:5621327/1‑93 (MQ=255)

                                                                                        aagGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGtt   <  2:323501/93‑1 (MQ=255)

GCGCGGGCACAG‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTT  >  NC_000913/3158190‑3158370

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