BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	474,643	G→T	78.7%	V173V (GTC→GTA)	tesB ←	acyl‑CoA thioesterase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	474,643	0	G	T	78.7%	263.1 / 78.9	127	V173V (GTC→GTA)	tesB	acyl‑CoA thioesterase 2
Reads supporting (aligned to +/- strand): ref base G (17/10); new base T (56/44); total (73/54)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

Only 100 of 127 total aligned reads displayed.

AAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGC  >  NC_000913/474547‑474741

aaGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACg                                                                                                 <  1:192978/100‑1 (MQ=255)

  gATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTg                                                                                               >  1:1412360/1‑100 (MQ=255)

  gATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTg                                                                                               >  1:868545/1‑100 (MQ=255)

    tCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGAc                                                                                             <  2:2383471/100‑1 (MQ=255)

    tCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCCGACGATGTGGTTCTGCTACGTGAc                                                                                             <  2:4124667/100‑1 (MQ=255)

      aGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCt                                                                                           >  1:1382183/1‑100 (MQ=255)

       gAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCtt                                                                                          >  1:3433621/1‑100 (MQ=255)

         aGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTc                                                                                        >  1:2466381/1‑100 (MQ=255)

         aGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTc                                                                                        <  1:1930506/100‑1 (MQ=255)

            gTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGt                                                                                     >  2:3436776/1‑100 (MQ=255)

             tgaCCGCGCCGATACTGAGGAACGCGCAGGCCATACGGCACGCTACCATTTGCGCGCATCCACCGCTGACGCTGTAGTTCTGCTACGGGACATTTCAg                                                                                      >  2:2585276/3‑98 (MQ=255)

             tAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTg                                                                                    >  1:3990444/1‑100 (MQ=255)

              acccGAGCAGATACTGATGGACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCTCACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTg                                                                                    <  2:3849772/97‑1 (MQ=255)

               aCCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTggg                                                                                  >  1:1855260/1‑100 (MQ=255)

               aCCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTggg                                                                                  <  1:3095452/100‑1 (MQ=255)

                ccGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGt                                                                                 <  2:2397166/100‑1 (MQ=255)

                 cGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGtt                                                                                >  2:577748/1‑100 (MQ=255)

                 cGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGtt                                                                                >  2:404250/1‑100 (MQ=255)

                  gAGCAGATACTGATGAACGCGCAGGTAATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTggg                                                                                  <  2:2874635/97‑1 (MQ=255)

                   aGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTAt                                                                              <  1:3721409/100‑1 (MQ=255)

                    gCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTAt                                                                              >  1:3119140/1‑99 (MQ=255)

                    gCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATg                                                                             <  1:3420570/100‑1 (MQ=255)

                     cAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGa                                                                            <  2:4094642/100‑1 (MQ=255)

                       gATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGaaa                                                                          >  2:443924/1‑100 (MQ=255)

                       gATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGaa                                                                           <  2:1277897/99‑1 (MQ=255)

                        aTACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAAc                                                                         <  1:2246586/100‑1 (MQ=255)

                        aTACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAAc                                                                         <  1:4065467/100‑1 (MQ=255)

                        aTACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAAc                                                                         <  1:3079847/100‑1 (MQ=255)

                         tACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACt                                                                        >  1:1943018/1‑100 (MQ=255)

                           cTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTcc                                                                      <  2:489238/100‑1 (MQ=255)

                             gatgaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCAc                                                                    >  1:2715141/1‑100 (MQ=255)

                             gatgaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCAc                                                                    >  1:164675/1‑100 (MQ=255)

                             gatgaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCAc                                                                    >  1:3610217/1‑100 (MQ=255)

                               tgaACGCGCAGGTCATGCGGCACGCTACCATTTGCTCGGATCTACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATCAACGTCCACcg                                                                  <  2:2942071/100‑1 (MQ=255)

                               tgaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACcg                                                                  >  1:1684424/1‑100 (MQ=255)

                               tgaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACcg                                                                  >  1:2167551/1‑100 (MQ=255)

                                gaACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACcgg                                                                 >  2:2777689/1‑100 (MQ=255)

                                 aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACcgg                                                                 >  1:154489/1‑99 (MQ=255)

                                 aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACcgga                                                                >  2:1326876/1‑100 (MQ=255)

                                 aaCGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACcgga                                                                <  1:2305277/100‑1 (MQ=255)

                                   cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACcggacg                                                              <  2:1925746/100‑1 (MQ=255)

                                   cgcgCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACcggacg                                                              >  1:402265/1‑100 (MQ=255)

                                        aGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACtt                                                         >  1:3631298/1‑100 (MQ=255)

                                         ggTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTc                                                        >  2:3839046/1‑100 (MQ=255)

                                         ggTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGCTATGAAACTCCACCGGACGGACTTc                                                        >  2:1003496/1‑100 (MQ=255)

                                           tCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCa                                                      <  1:629330/100‑1 (MQ=255)

                                            cATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAg                                                     <  1:2616286/100‑1 (MQ=255)

                                             aTCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGc                                                    >  1:3648196/1‑100 (MQ=255)

                                              tCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCg                                                   >  2:4100792/1‑100 (MQ=255)

                                              tCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCg                                                   >  1:2319918/1‑100 (MQ=255)

                                              tCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCg                                                   <  1:7995/100‑1 (MQ=255)

                                              tCCGACACCCAACCATTTGTGCGCATCCACACCTGACTATGTGGTACTTCTACGTCACCGTTCAGTGGGTTATTAAACTCCACCGGACGGACTTCCAGCg                                                   <  2:1057544/100‑1 (MQ=255)

                                               ccGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCgg                                                  >  2:3264353/1‑100 (MQ=255)

                                                cGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGa                                                 <  1:2394547/100‑1 (MQ=255)

                                                 ggCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGAc                                                >  2:3888985/1‑100 (MQ=255)

                                                 ggCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGAc                                                >  2:1557436/1‑100 (MQ=255)

                                                  gCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACg                                               >  2:2608801/1‑100 (MQ=255)

                                                    aCGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGAt                                             >  1:3274015/1‑100 (MQ=255)

                                                     cGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATc                                            >  1:2143893/1‑100 (MQ=255)

                                                      gCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCg                                           <  1:3345016/100‑1 (MQ=255)

                                                      gCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCg                                           >  1:1540198/1‑100 (MQ=255)

                                                       cTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGc                                          >  2:1312571/1‑100 (MQ=255)

                                                       cTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGc                                          <  1:3866291/100‑1 (MQ=255)

                                                           cATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAg                                        >  2:2193652/1‑98 (MQ=255)

                                                            aTTTGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATg                                     >  1:1917846/1‑100 (MQ=255)

                                                             tttGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGa                                    >  2:3105905/1‑100 (MQ=255)

                                                               tGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAAt                                  >  1:3796577/1‑100 (MQ=255)

                                                               tGCGCGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAAt                                  <  1:2442447/100‑1 (MQ=255)

                                                                gcgcGGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAAtt                                 <  1:2555757/100‑1 (MQ=255)

                                                                gcgcGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAAtt                                 >  2:1632824/1‑100 (MQ=255)

                                                                 cgcgGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAAttt                                >  1:1424793/1‑100 (MQ=255)

                                                                  gcgGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGCTTCTGAAACTCCACCGGCCGGACTTCCAGCGGACGACCGCAGATGAATTTa                               >  2:1549262/1‑100 (MQ=255)

                                                                   cgGATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTAt                              <  2:1623679/100‑1 (MQ=255)

                                                                   cgGATCCACACCTGACGATGTGGTTCCTCTACGTGCCCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGAACGCAGATGAATTTAt                              <  2:3481140/100‑1 (MQ=255)

                                                                     gATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTAt                              >  1:3476139/1‑98 (MQ=255)

                                                                     gATCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCt                            >  1:511050/1‑100 (MQ=255)

                                                                      aTCCACACCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATAAATTTAt                              >  2:3452655/1‑97 (MQ=255)

                                                                        ccACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTc                         <  1:469220/100‑1 (MQ=255)

                                                                         cacaCCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTCGGTTATGAAACTCCCCCGGACGGACTTCCACCGGACGATCGCAGATGCATTTATCTTTCa                        >  2:1851642/1‑100 (MQ=255)

                                                                          acacCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAg                       <  2:2783581/100‑1 (MQ=255)

                                                                          acacCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAg                       <  1:3400123/100‑1 (MQ=255)

                                                                          acacCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAg                       <  1:553731/100‑1 (MQ=255)

                                                                           cacCTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGc                      <  1:2347940/100‑1 (MQ=255)

                                                                             ccTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCAc                    >  2:3495323/1‑100 (MQ=255)

                                                                             ccTGACGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCAc                    <  1:3435347/100‑1 (MQ=255)

                                                                                  cGATGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTggcg               >  2:3096874/1‑100 (MQ=255)

                                                                                  cGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTggcg               >  2:887306/1‑100 (MQ=255)

                                                                                    aTGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTggcggc             <  2:2725145/100‑1 (MQ=255)

                                                                                    aTGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTggcggc             >  1:2725161/1‑100 (MQ=255)

                                                                                    aTGTGGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTggcggc             <  1:1672536/100‑1 (MQ=255)

                                                                                     tgtgGTTCTGCTACGTGGTCTTTCAGTGGGTTATGAAACTCCGCCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGgca            <  2:2935545/100‑1 (MQ=255)

                                                                                     tgtgGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGgca            >  1:2176135/1‑100 (MQ=255)

                                                                                       tgGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGgcagc          >  2:1941109/1‑100 (MQ=255)

                                                                                       tgGTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGgcagc          >  1:1090014/1‑100 (MQ=255)

                                                                                        ggTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGgcagca         <  1:3839075/100‑1 (MQ=255)

                                                                                        ggTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGgcagca         >  1:1313238/1‑100 (MQ=255)

                                                                                         gTTCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGcagcag        >  1:116065/1‑100 (MQ=255)

                                                                                          ttCTGCTACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAgg       <  2:2988880/100‑1 (MQ=255)

                                                                                              gCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTgcg   >  1:3294718/1‑100 (MQ=255)

                                                                                               cGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTgcgc  >  1:3070810/1‑100 (MQ=255)

AAGATCAGAAGCGTAACCGAGCAGATACTGATGAACGCGCAGGTCATCCGGCACGCTACCATTTGCGCGGATCCACACCTGACGATGTGGTTCTGCGACGTGACCTTTCAGTGGGTTATGAAACTCCACCGGACGGACTTCCAGCGGACGATCGCAGATGAATTTATCTTTCAGCACTGGCGGCAGCAGGTGCGC  >  NC_000913/474547‑474741

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