BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,578,097	A→G	T148A (ACC→GCC)	tktB →	transketolase 2 isozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,578,097	0	A	G	85.0%	659.8 / 95.8	280	T148A (ACC→GCC)	tktB	transketolase 2 isozyme
Reads supporting (aligned to +/- strand): ref base A (21/21); new base G (117/121); total (138/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 280 total aligned reads displayed.

AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGA  >  NC_000913/2578006‑2578187

  gTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCt                                                                                         >  1:452344/1‑93 (MQ=255)

    ttGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTAt                                                                                       <  1:6105341/93‑1 (MQ=255)

     tGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTAtg                                                                                      <  2:6898683/93‑1 (MQ=255)

       gCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTAtgtg                                                                                    <  1:6438661/93‑1 (MQ=255)

        cGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATgtgt                                                                                   <  1:2270784/93‑1 (MQ=255)

         gAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGtt                                                                                  >  2:4313897/1‑93 (MQ=255)

          aaCGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGttt                                                                                 >  1:5158553/1‑93 (MQ=255)

           aCGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGTCTATGTGTTTa                                                                                >  1:4750541/1‑93 (MQ=255)

           aCGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTa                                                                                <  1:6270946/93‑1 (MQ=255)

           aCGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTa                                                                                <  2:5804180/93‑1 (MQ=255)

           aCGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTa                                                                                >  2:3757947/1‑93 (MQ=255)

           aCGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTa                                                                                >  1:3389412/1‑93 (MQ=255)

             gCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATg                                                                              >  1:3857207/1‑93 (MQ=255)

              ccGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATgg                                                                             <  2:3054223/93‑1 (MQ=255)

              ccGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATgg                                                                             <  1:385114/93‑1 (MQ=255)

              ccGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATgg                                                                             >  2:4905287/1‑93 (MQ=255)

                gTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGc                                                                           >  2:1941435/1‑93 (MQ=255)

                gTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGc                                                                           >  1:5424410/1‑93 (MQ=255)

                 tCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCg                                                                          >  1:6641961/1‑93 (MQ=255)

                   gggCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGAc                                                                        <  1:4371786/93‑1 (MQ=255)

                   gggCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGAc                                                                        >  1:4936449/1‑93 (MQ=255)

                   gggCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGAc                                                                        >  2:7096445/1‑93 (MQ=255)

                      cTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGc                                                                     <  2:2110664/93‑1 (MQ=255)

                      cTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGc                                                                     <  2:275785/93‑1 (MQ=255)

                       tGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCt                                                                    >  1:3134526/1‑93 (MQ=255)

                          cGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGcc                                                                 <  2:1272595/93‑1 (MQ=255)

                           gATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCt                                                                <  1:3535115/93‑1 (MQ=255)

                           gATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCt                                                                <  1:2337051/93‑1 (MQ=255)

                              aGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGAt                                                             <  1:1060424/93‑1 (MQ=255)

                               gCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATg                                                            <  2:3462883/93‑1 (MQ=255)

                               gCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATg                                                            >  1:1970628/1‑93 (MQ=255)

                               gCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATg                                                            <  2:1074869/93‑1 (MQ=255)

                                cAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATgg                                                           <  2:4149351/93‑1 (MQ=255)

                                 agagCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGa                                                          >  1:5643526/1‑93 (MQ=255)

                                 agagCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGa                                                          <  2:6102785/93‑1 (MQ=255)

                                  gagCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGaa                                                         <  2:1807978/93‑1 (MQ=255)

                                  gagCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGaa                                                         >  2:4026323/1‑93 (MQ=255)

                                    gCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAgg                                                       >  2:2650177/1‑93 (MQ=255)

                                    gCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAgg                                                       <  1:6310943/93‑1 (MQ=255)

                                    gCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAgg                                                       >  2:5118519/1‑93 (MQ=255)

                                      gTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTa                                                     >  1:4609475/1‑93 (MQ=255)

                                      gTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTa                                                     <  2:481461/93‑1 (MQ=255)

                                      gTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTa                                                     >  1:3934420/1‑93 (MQ=255)

                                      gTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTa                                                     <  2:6809016/93‑1 (MQ=255)

                                        acacTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTAtt                                                   <  1:6708669/93‑1 (MQ=255)

                                         cacTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTAttt                                                  >  2:5917741/1‑93 (MQ=255)

                                            tGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTccc                                               >  1:2151546/1‑93 (MQ=255)

                                            tGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTccc                                               >  1:2978186/1‑93 (MQ=255)

                                             ggcggcGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCa                                              <  1:5919201/93‑1 (MQ=255)

                                              gcggcgCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCAc                                             <  1:2646449/93‑1 (MQ=255)

                                                 gcgcAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGaa                                          <  1:3730668/93‑1 (MQ=255)

                                                  cgcAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAg                                         <  2:3158961/93‑1 (MQ=255)

                                                   gcAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGt                                        >  2:4876541/1‑93 (MQ=255)

                                                   gcAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGt                                        >  2:2442132/1‑93 (MQ=255)

                                                   gcAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGt                                        <  2:6057077/93‑1 (MQ=255)

                                                    cAGTTTAACCAGCCAGACCATGAGGTCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTc                                       >  2:2862236/1‑93 (MQ=255)

                                                    cAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTc                                       >  2:6506516/1‑93 (MQ=255)

                                                     aGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCt                                      >  1:4086655/1‑93 (MQ=255)

                                                     aGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCt                                      >  1:3185437/1‑93 (MQ=255)

                                                      gTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTg                                     <  2:5498623/93‑1 (MQ=255)

                                                      gTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTg                                     >  1:6958895/1‑93 (MQ=255)

                                                      gTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTg                                     >  2:4471428/1‑93 (MQ=255)

                                                       tttAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGt                                    <  1:7387912/93‑1 (MQ=255)

                                                       tttAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGt                                    >  2:7023760/1‑93 (MQ=255)

                                                       tttAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGt                                    <  2:4345993/93‑1 (MQ=255)

                                                        ttAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGtt                                   <  1:6071207/93‑1 (MQ=255)

                                                         tAACCAGCCAGACCATGAGATGGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCCCGAAGTCTGTTc                                  <  1:7443421/93‑1 (MQ=255)

                                                         tAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCTCACGAAGTCTGTTc                                  <  2:5196189/93‑1 (MQ=255)

                                                         tAACCAGCCAGACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTc                                  <  2:6445447/93‑1 (MQ=255)

                                                           aCCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGc                                >  2:5645573/1‑93 (MQ=255)

                                                             cagccagACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTg                              >  1:550561/1‑93 (MQ=255)

                                                             cagccagACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTg                              <  2:2964785/93‑1 (MQ=255)

                                                             cagccagACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCta                              <  1:58951/93‑2 (MQ=255)

                                                                ccagACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggca                           <  2:7199179/93‑1 (MQ=255)

                                                                ccagACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggca                           <  2:1010917/93‑1 (MQ=255)

                                                                 cagACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcag                          >  2:1499454/1‑93 (MQ=255)

                                                                   gACCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcaggc                        >  1:1527497/1‑93 (MQ=255)

                                                                    aCCATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcaggca                       >  1:7473687/1‑93 (MQ=255)

                                                                     ccATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCAc                      >  2:5610018/1‑93 (MQ=255)

                                                                     ccATGAGATCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCAc                      <  1:343446/93‑1 (MQ=255)

                                                                       aTGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGc                    >  1:1252491/1‑93 (MQ=255)

                                                                        tGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCt                   <  1:3268230/93‑1 (MQ=255)

                                                                            aTCGTCGATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGa               >  2:4526036/1‑93 (MQ=255)

                                                                             tcgtcgATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGAc              <  1:1021881/93‑1 (MQ=255)

                                                                             tcgtcgATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGAc              <  1:1459513/93‑1 (MQ=255)

                                                                             tcgtcgATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGAc              >  2:767937/1‑93 (MQ=255)

                                                                              cgtcgATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACt             <  2:3886544/93‑1 (MQ=255)

                                                                              cgtcgATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACt             >  1:1863826/1‑93 (MQ=255)

                                                                                tcgATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTgg           <  2:3309647/93‑1 (MQ=255)

                                                                                 cgATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTggg          >  2:906051/1‑93 (MQ=255)

                                                                                 cgATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTggg          >  2:1050696/1‑93 (MQ=255)

                                                                                 cgATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTggg          >  2:2407085/1‑93 (MQ=255)

                                                                                  gATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGc         <  2:5996748/93‑1 (MQ=255)

                                                                                  gATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGc         >  2:5442026/1‑93 (MQ=255)

                                                                                  gATCACTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGc         >  1:705649/1‑93 (MQ=255)

                                                                                      aCTTCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGc     >  2:6188982/1‑93 (MQ=255)

                                                                                        ttCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTg   >  2:6266392/1‑93 (MQ=255)

                                                                                        ttCGCCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTg   >  2:4064378/1‑93 (MQ=255)

                                                                                        ttCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTg   >  2:7470866/1‑93 (MQ=255)

                                                                                         tCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGa  >  1:3628173/1‑93 (MQ=255)

AGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGA  >  NC_000913/2578006‑2578187

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: `‑`