BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,298,169	T→C	T38A (ACC→GCC)	napG ←	ferredoxin‑type protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,298,169	0	T	C	87.8%	586.6 / 62.7	239	T38A (ACC→GCC)	napG	ferredoxin‑type protein
Reads supporting (aligned to +/- strand): ref base T (14/15); new base C (103/106); total (118/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

Only 100 of 242 total aligned reads displayed.

AACACACTGACCGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAGCGGCGGCGACCATTTT  >  NC_000913/2298077‑2298258

    cacTGACCGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGc                                                                                       >  2:1289214/1‑93 (MQ=255)

     acTGACCGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCt                                                                                      <  1:2303556/93‑1 (MQ=255)

      cTGACCGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCtg                                                                                     <  2:1762906/93‑1 (MQ=255)

      cTGACCGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCtg                                                                                     <  1:3968724/93‑1 (MQ=255)

        gACCGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGTTTGCtgtt                                                                                   >  1:141143/1‑93 (MQ=255)

        gACCGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCtgtt                                                                                   <  2:1848103/93‑1 (MQ=255)

         aCCGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCtgttg                                                                                  >  1:5278505/1‑93 (MQ=255)

         aCCGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCtgttg                                                                                  >  1:5321972/1‑93 (MQ=255)

          ccGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCtgttgt                                                                                 <  1:1616673/93‑1 (MQ=255)

          ccGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCtgttgt                                                                                 <  2:5509709/93‑1 (MQ=255)

          ccGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCtgttgt                                                                                 <  1:5919455/93‑1 (MQ=255)

           cGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTa                                                                                >  2:6138216/1‑93 (MQ=255)

            gCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTaa                                                                               <  1:5240925/93‑1 (MQ=255)

            gCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTaa                                                                               <  2:7202802/93‑1 (MQ=255)

             cAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAAc                                                                              <  2:634472/93‑1 (MQ=255)

             cAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAAc                                                                              <  2:6513991/93‑1 (MQ=255)

              aaCGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGTTTGCTGTTGTAAcc                                                                             <  2:7002601/93‑1 (MQ=255)

              aaCGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAAcc                                                                             >  1:1116597/1‑93 (MQ=255)

              aaCGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAAcc                                                                             >  1:480945/1‑93 (MQ=255)

                 gCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCa                                                                          <  1:922603/93‑1 (MQ=255)

                       aGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCca                                                                    <  2:764761/93‑1 (MQ=255)

                         gCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCcaca                                                                  <  2:5299205/93‑1 (MQ=255)

                          cACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCacac                                                                 >  2:5859803/1‑93 (MQ=255)

                           aCTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACAcc                                                                <  2:4811481/93‑1 (MQ=255)

                            cTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACAccc                                                               <  1:4987659/93‑1 (MQ=255)

                             tGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCa                                                              >  1:7296983/1‑93 (MQ=255)

                             tGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAAACCCAGCGCCACACCCa                                                              >  1:4244639/1‑93 (MQ=255)

                               gCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACg                                                            >  2:1188417/1‑93 (MQ=255)

                               gCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACg                                                            >  2:3559505/1‑93 (MQ=255)

                               gCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACg                                                            >  2:6797805/1‑93 (MQ=255)

                               gCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACg                                                            >  1:447866/1‑93 (MQ=255)

                                 aaaGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGc                                                          <  2:3003095/93‑1 (MQ=255)

                                 aaaGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGc                                                          >  2:2653772/1‑93 (MQ=255)

                                 aaaGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGc                                                          >  1:4948314/1‑93 (MQ=255)

                                  aaGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGca                                                         >  2:1361167/1‑93 (MQ=255)

                                  aaGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGca                                                         <  1:682964/93‑1 (MQ=255)

                                    ggCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGcagc                                                       >  1:3264033/1‑93 (MQ=255)

                                    ggCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGcagc                                                       <  1:6228381/93‑1 (MQ=255)

                                     gCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGcagcc                                                      <  2:2444402/93‑1 (MQ=255)

                                       gTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGcagccag                                                    >  2:3110218/1‑93 (MQ=255)

                                        ttCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACGGcagccagc                                                   <  2:3316943/93‑1 (MQ=255)

                                        ttCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGcagccagc                                                   >  2:3928721/1‑93 (MQ=255)

                                        ttCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGcagccagc                                                   >  1:2714219/1‑93 (MQ=255)

                                         tctcGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACGGCagccagcc                                                  <  2:5019892/93‑1 (MQ=255)

                                         tctcGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACGGCagccagcc                                                  <  2:5991177/93‑1 (MQ=255)

                                         tctcGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCagccagcc                                                  >  2:6171635/1‑93 (MQ=255)

                                           tcGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCg                                                >  2:7414593/1‑93 (MQ=255)

                                            cGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGc                                               <  2:4499739/93‑1 (MQ=255)

                                            cGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGc                                               <  2:7028756/93‑1 (MQ=255)

                                             gTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGcc                                              >  1:3633554/1‑93 (MQ=255)

                                              tttATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGccc                                             >  1:2427275/1‑92 (MQ=255)

                                              tttATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCct                                             >  2:3489763/1‑93 (MQ=255)

                                              tttATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCct                                             >  2:3253898/1‑93 (MQ=255)

                                                   ggCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACCCCCACGGCAGCCAGCCCGCCTGCtgt                                        <  2:3215463/93‑1 (MQ=255)

                                                   ggCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCtgt                                        <  2:2329143/93‑1 (MQ=255)

                                                      tCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCg                                     >  1:6393122/1‑93 (MQ=255)

                                                       ccGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGa                                    <  1:4870729/93‑1 (MQ=255)

                                                        cgggcgggCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGaa                                   <  2:6699013/93‑1 (MQ=255)

                                                        cgggcgggCGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGaa                                   <  1:4393416/93‑1 (MQ=255)

                                                         gggcgggcGCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAAc                                  >  2:2155300/1‑93 (MQ=255)

                                                            cgggcgCAACCGCACGCCAGATGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGAc                               <  1:7097064/93‑1 (MQ=255)

                                                             gggcgCAACCGCACGCAAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACa                              <  1:50630/93‑1 (MQ=255)

                                                              ggcgCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACAt                             <  2:4731274/93‑1 (MQ=255)

                                                              ggcgCAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACAt                             >  1:2968961/1‑93 (MQ=255)

                                                               gcgcAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATc                            >  2:3759244/1‑93 (MQ=255)

                                                                 gcAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATcgc                          <  2:5489788/93‑1 (MQ=255)

                                                                 gcAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATcgc                          <  2:4977158/93‑1 (MQ=255)

                                                                  cAACCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATcgcg                         >  2:2516738/1‑93 (MQ=255)

                                                                   aaCCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCgcgc                        <  2:3346789/93‑1 (MQ=255)

                                                                    aCCGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCa                       <  1:4044037/93‑1 (MQ=255)

                                                                     ccGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAg                      >  1:830151/1‑93 (MQ=255)

                                                                     ccGCACGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAg                      <  2:2840548/93‑1 (MQ=255)

                                                                          cGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgc                 <  2:1001182/93‑1 (MQ=255)

                                                                          cGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgc                 <  1:6895757/93‑1 (MQ=255)

                                                                          cGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgc                 <  1:6100498/93‑1 (MQ=255)

                                                                          cGCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgc                 <  1:4719857/93‑1 (MQ=255)

                                                                           gCCAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcg                <  1:2121368/93‑1 (MQ=255)

                                                                            ccAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcgg               >  1:7502920/1‑93 (MQ=255)

                                                                             nngATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggc              <  1:3396897/91‑1 (MQ=255)

                                                                             cAGATGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggc              <  1:5367439/93‑1 (MQ=255)

                                                                             cAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggc              <  2:6946574/93‑1 (MQ=255)

                                                                             cAGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggc              >  2:5787575/1‑93 (MQ=255)

                                                                              aGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggcg             >  2:6440980/1‑93 (MQ=255)

                                                                              aGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggcg             >  1:1337523/1‑93 (MQ=255)

                                                                              aGATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggcg             <  1:3139087/93‑1 (MQ=255)

                                                                               gATGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggcgg            >  1:1661193/1‑93 (MQ=255)

                                                                               gATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggcgg            <  2:5366571/93‑1 (MQ=255)

                                                                               gATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggcgg            >  2:1671189/1‑93 (MQ=255)

                                                                               gATGCGCGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggcgg            >  1:794103/1‑93 (MQ=255)

                                                                                aTGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggcggc           <  1:4177482/93‑1 (MQ=255)

                                                                                 tGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAgcggcggcg          >  2:4233156/1‑93 (MQ=255)

                                                                                  gcgcGTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAGCGGCGGCGa         >  1:947946/1‑93 (MQ=255)

                                                                                   cgcgTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAGCGGCGGCGAc        >  2:1486304/1‑93 (MQ=255)

                                                                                     cgTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAGCGGCGGCGACCa      <  2:2570951/93‑1 (MQ=255)

                                                                                     cgTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAGCGGCGGCGACCa      >  1:5618422/1‑93 (MQ=255)

                                                                                     cgTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAGCGGCGGCGACCa      >  1:3474409/1‑93 (MQ=255)

                                                                                      gTGCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAGCGGCGGCGACCAt     <  2:7309448/93‑1 (MQ=255)

                                                                                       tGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAGCGGCGGCGACCAtt    >  2:4087593/1‑93 (MQ=255)

                                                                                        gCGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAGCGGCGGCGACCAttt   >  2:4311160/1‑93 (MQ=255)

                                                                                         cGGCTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAGCGGCGGCGACCAtttt  <  2:477276/93‑1 (MQ=255)

AACACACTGACCGCAACGCACACAGGCACTGGCAAAGGCGTTCTCGTTTATGGCTCCGGGCGGGCGCAACCGCACGCCAGATGCGCGTGCGGTTTGCTGTTGTAACCCCAGCGCCACACCCACGGCAGCCAGCCCGCCTGCTGTGCGAACGACATCGCGCAGAAAGCGGCGGCGACCATTTT  >  NC_000913/2298077‑2298258

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