BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,101,591	C→T	A322A (GCG→GCA)	cpxA ←	CpxA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,101,591	0	C	T	83.7%	714.0 / 103.7	264	A322A (GCG→GCA)	cpxA	CpxA
Reads supporting (aligned to +/- strand): ref base C (16/27); new base T (121/100); total (137/127)
Fisher's exact test for biased strand distribution p-value = 4.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

Only 100 of 282 total aligned reads displayed.

CACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGTGACATCAC  >  NC_000913/4101500‑4101680

cACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAAcg                                                                                          >  2:1687056/1‑93 (MQ=255)

 aCTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAAtgc                                                                                         <  2:4483103/93‑1 (MQ=255)

  cTTTCCAGGTCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAAtgcc                                                                                        <  2:5848126/93‑1 (MQ=255)

      ccAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTATCTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGt                                                                                    <  2:4667389/93‑1 (MQ=255)

      ccAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGt                                                                                    <  1:3497928/93‑1 (MQ=255)

       cAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGtt                                                                                   >  2:5087831/1‑93 (MQ=255)

       cAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGtt                                                                                   >  1:5611438/1‑93 (MQ=255)

        aGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTa                                                                                  >  1:5204149/1‑93 (MQ=255)

        aGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTa                                                                                  >  2:2871424/1‑93 (MQ=255)

         gggCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTAt                                                                                 >  2:988164/1‑93 (MQ=255)

          ggCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATc                                                                                >  1:224509/1‑93 (MQ=255)

          ggCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATc                                                                                <  2:1870543/93‑1 (MQ=255)

          ggCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATc                                                                                >  1:2289637/1‑93 (MQ=255)

           gCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATcc                                                                               <  1:4618972/93‑1 (MQ=255)

           gCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATcc                                                                               <  2:401035/93‑1 (MQ=255)

             gTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGCTAACTGTCAACGACTTGGCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAg                                                                             >  2:6347688/1‑93 (MQ=255)

                cGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGGAGTTAACTGGCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCAc                                                                          >  2:6595207/1‑93 (MQ=255)

                  gATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACtt                                                                        >  1:2558207/1‑93 (MQ=255)

                  gATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACtt                                                                        <  2:2514251/93‑1 (MQ=255)

                  gATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACtt                                                                        >  1:2209744/1‑93 (MQ=255)

                  gATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACtt                                                                        >  2:3873209/1‑93 (MQ=255)

                   aTTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTTTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTc                                                                       >  2:417754/1‑93 (MQ=255)

                   aTTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTc                                                                       >  1:1609063/1‑93 (MQ=255)

                     tGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCAc                                                                     <  1:1509839/93‑1 (MQ=255)

                      gCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACt                                                                    <  1:4962663/93‑1 (MQ=255)

                       ccGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTc                                                                   >  2:6071413/1‑93 (MQ=255)

                        cGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTcc                                                                  <  2:5776004/93‑1 (MQ=255)

                         gTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCca                                                                 >  1:6434559/1‑93 (MQ=255)

                          tACAGCGGCCACGGCCCAGGCGGGATGGTAACTGTCAACGACTTGTCCATTTGCTCGGCAGCGAATGCCGCGTTTCTTAGCACTTCTCTCcac                                                                >  1:6640684/1‑93 (MQ=255)

                            cAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACaa                                                              <  1:6262655/93‑1 (MQ=255)

                            cAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACaa                                                              <  2:6215315/93‑1 (MQ=255)

                             aGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAAc                                                             >  2:6355902/1‑93 (MQ=255)

                              gCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACt                                                            >  2:3048693/1‑93 (MQ=255)

                              gCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACt                                                            >  1:3158253/1‑93 (MQ=255)

                                ggCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTgg                                                          <  1:5300256/93‑1 (MQ=255)

                                ggCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTgg                                                          <  2:2297210/93‑1 (MQ=255)

                                   cACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTg                                                       >  2:4295689/1‑93 (MQ=255)

                                   cACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTg                                                       <  1:6004795/93‑1 (MQ=255)

                                     cGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGc                                                     <  1:6147057/93‑1 (MQ=255)

                                      ggCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCt                                                    <  1:1015621/93‑1 (MQ=255)

                                      ggCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCt                                                    <  2:6149074/93‑1 (MQ=255)

                                       gCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCtt                                                   <  1:1256679/93‑1 (MQ=255)

                                        cccAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCttt                                                  <  1:3017719/93‑1 (MQ=255)

                                         ccAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGTCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTg                                                 >  1:1070462/1‑93 (MQ=255)

                                          cAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGa                                                <  1:4644520/93‑1 (MQ=255)

                                          cAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGa                                                >  1:3659146/1‑93 (MQ=255)

                                           aGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGAt                                               <  2:5585608/93‑1 (MQ=255)

                                           aGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGAt                                               <  1:960816/93‑1 (MQ=255)

                                           aGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGAt                                               <  1:2840158/93‑1 (MQ=255)

                                           aGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGAt                                               <  1:2687534/93‑1 (MQ=255)

                                            ggCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATg                                              >  2:277332/1‑93 (MQ=255)

                                              cGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGt                                            >  1:5764418/1‑93 (MQ=255)

                                               gggAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGtt                                           >  2:2128013/1‑93 (MQ=255)

                                               gggAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGtt                                           >  2:3382432/1‑93 (MQ=255)

                                               gggAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGtt                                           >  1:384594/1‑93 (MQ=255)

                                               gggAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGtt                                           <  1:2736112/93‑1 (MQ=255)

                                                ggAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGttt                                          <  2:83344/93‑1 (MQ=255)

                                                ggAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGttt                                          >  2:5344304/1‑93 (MQ=255)

                                                  aaGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGtttac                                        >  2:4652195/1‑91 (MQ=255)

                                                  aaGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCg                                        >  2:4906357/1‑93 (MQ=255)

                                                  aaGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCg                                        >  1:1914589/1‑93 (MQ=255)

                                                   aGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGc                                       <  2:5190716/93‑1 (MQ=255)

                                                    gTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCt                                      >  2:3677326/1‑93 (MQ=255)

                                                    gTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCt                                      >  2:6309265/1‑93 (MQ=255)

                                                      tAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTaa                                    <  1:3697206/93‑1 (MQ=255)

                                                        aCTGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAAcc                                  >  2:2691420/1‑93 (MQ=255)

                                                        aCTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAAcc                                  >  1:3522808/1‑93 (MQ=255)

                                                        aCTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAAcc                                  >  2:5391033/1‑93 (MQ=255)

                                                          tGTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAg                                <  1:5047540/93‑1 (MQ=255)

                                                           gTCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAgc                               <  1:4066418/93‑1 (MQ=255)

                                                            tCAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAgcg                              <  1:3298318/93‑1 (MQ=255)

                                                            tCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAgcg                              <  1:4899782/93‑1 (MQ=255)

                                                             cAACGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCTCTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAgcgc                             <  2:1441564/93‑1 (MQ=255)

                                                              aaCGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGcgcg                            <  2:4091338/93‑1 (MQ=255)

                                                               aCGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGt                           >  2:5487294/1‑93 (MQ=255)

                                                               aCGACTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGt                           <  1:3725144/93‑1 (MQ=255)

                                                               aCGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGt                           <  2:3216233/93‑1 (MQ=255)

                                                                  aCTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGtttt                        >  2:2503654/1‑93 (MQ=255)

                                                                   cTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGttttt                       >  1:4741607/1‑93 (MQ=255)

                                                                   cTTGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGttttt                       <  1:14406/93‑1 (MQ=255)

                                                                    ttGCCCATTTGCTCTGCTTCGAATGCCTCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGtttttt                      <  1:5639742/93‑1 (MQ=255)

                                                                    ttGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGtttttt                      <  2:5887409/93‑1 (MQ=255)

                                                                     tGCCCATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGttttttt                     <  1:1389673/93‑1 (MQ=255)

                                                                        ccATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCt                  <  2:1196120/93‑1 (MQ=255)

                                                                        ccATTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCt                  >  1:1288024/1‑93 (MQ=255)

                                                                          aTTTGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGa                <  2:293331/93‑1 (MQ=255)

                                                                           tttGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGAt               <  1:5569019/93‑1 (MQ=255)

                                                                           tttGCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGAt               >  1:6505831/1‑93 (MQ=255)

                                                                              gCTCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTAc            >  2:6512378/1‑93 (MQ=255)

                                                                                tCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGt          >  1:2461221/1‑93 (MQ=255)

                                                                                tCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGt          >  1:5933314/1‑93 (MQ=255)

                                                                                tCGGCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGt          <  1:2314475/93‑1 (MQ=255)

                                                                                   gCTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGTGAc       <  1:5721586/93‑1 (MQ=255)

                                                                                   gCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGTGAc       <  2:1793773/93‑1 (MQ=255)

                                                                                    cTTCGAATGCCGCGTTATCCAGCACTTCCCTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGTcaca      <  1:4437831/93‑1 (MQ=255)

                                                                                    cTTCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGTGACa      >  1:4217418/1‑93 (MQ=255)

                                                                                      tCGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGTGACATc    >  1:3326350/1‑93 (MQ=255)

                                                                                       cGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGTGACATCa   >  1:2813416/1‑93 (MQ=255)

                                                                                       cGAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGTGACATCa   <  1:1192535/93‑1 (MQ=255)

                                                                                        gAATGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGTGACATCAc  <  2:1272778/93‑1 (MQ=255)

CACTTTCCAGGGCGTTCGGATTGCCGTACAGCGGCCACGGCCCAGGCGGGAAGTTAACTGTCAACGACTTGCCCATTTGCTCGGCTTCGAACGCCGCGTTATCCAGCACTTCACTCCACAACTGGTTGGCTTTGATGGTTTCGCTAACCAGCGCGTTTTTTTGCTGATTACGTGACATCAC  >  NC_000913/4101500‑4101680

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