BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,597,539	T→C	intergenic (‑250/+21)	b3460 ← / ← rpoH	LivJ/sigma(32) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,597,539	0	T	C	99.8%	1735.4 / NA	537	intergenic (‑250/+21)	b3460/rpoH	LivJ/sigma(32) factor
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (236/300); minor base G (1/0); total (237/300)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 553 total aligned reads displayed.

ACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGA  >  NC_000913/3597449‑3597628

 cAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGctct                                                                                        <  2:4719649/93‑1 (MQ=255)

  aaTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGTGGCCCAAAAAACAAAAACCCCGTACTTTCATCCAGGGCtctc                                                                                       >  2:8966961/1‑93 (MQ=255)

  aaTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtctc                                                                                       <  2:5890640/93‑1 (MQ=255)

  aaTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtctc                                                                                       <  1:5777259/93‑1 (MQ=255)

   aTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGctctct                                                                                      <  2:5287356/93‑1 (MQ=255)

     cTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGc                                                                                    <  1:9988146/93‑1 (MQ=255)

      tGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCt                                                                                   <  2:761290/93‑1 (MQ=255)

      tGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCt                                                                                   >  2:3235831/1‑93 (MQ=255)

       gCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCtt                                                                                  <  1:8252719/93‑1 (MQ=255)

         gTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTaa                                                                                >  1:3355655/1‑93 (MQ=255)

          tGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCCGGGCTCTCTGCTTAAt                                                                               >  2:980630/1‑93 (MQ=255)

          tGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAAt                                                                               >  1:8889807/1‑93 (MQ=255)

           ggggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATa                                                                              <  2:3801416/93‑1 (MQ=255)

            gggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAg                                                                             <  1:6435841/93‑1 (MQ=255)

              gATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCg                                                                           <  1:7480137/93‑1 (MQ=255)

               aTTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCgg                                                                          <  2:7736844/93‑1 (MQ=255)

                ttGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGa                                                                         >  2:419935/1‑93 (MQ=255)

                  ggCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGaaa                                                                       <  1:6952631/93‑1 (MQ=255)

                   gCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAAt                                                                      >  1:808577/1‑93 (MQ=255)

                   gCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAAt                                                                      <  2:5129454/93‑1 (MQ=255)

                   gCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAAt                                                                      >  2:4486761/1‑93 (MQ=255)

                     gTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTa                                                                    >  2:6834711/1‑93 (MQ=255)

                      ttttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTAc                                                                   >  2:4964472/1‑93 (MQ=255)

                         tGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCt                                                                <  2:4087881/93‑1 (MQ=255)

                          gCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCtt                                                               >  1:5983963/1‑93 (MQ=255)

                           gcGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTc                                                              <  1:6262026/92‑1 (MQ=255)

                           ccGGAGGGGAAATTGATTATTACAGAGGCCCAAGAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTc                                                              <  1:9705518/93‑1 (MQ=255)

                           ccGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTc                                                              >  2:9642752/1‑93 (MQ=255)

                            cGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCa                                                             <  2:9429481/93‑1 (MQ=255)

                            cGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCAGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCa                                                             >  2:275394/1‑93 (MQ=255)

                             ggAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCaa                                                            <  1:6354944/93‑1 (MQ=255)

                             ggAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCaa                                                            >  1:1097136/1‑93 (MQ=255)

                             ggAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCaa                                                            <  2:1892754/93‑1 (MQ=255)

                              gAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAAt                                                           <  2:8651072/93‑1 (MQ=255)

                              gAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAAt                                                           <  1:656859/93‑1 (MQ=255)

                               aGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATg                                                          >  1:2302826/1‑93 (MQ=255)

                                 gggAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgc                                                        <  2:3394640/93‑1 (MQ=255)

                                 gggAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgc                                                        >  1:4412460/1‑93 (MQ=255)

                                  ggAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgca                                                       <  1:607341/93‑1 (MQ=255)

                                    aaaTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcagc                                                     >  2:4600733/1‑93 (MQ=255)

                                      aTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCAc                                                   <  2:8211449/93‑1 (MQ=255)

                                      aTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCAc                                                   <  2:6639310/93‑1 (MQ=255)

                                      aTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGCCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCAc                                                   <  2:645564/93‑1 (MQ=255)

                                       ttGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACg                                                  >  1:9246823/1‑93 (MQ=255)

                                          attattACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCaa                                               <  2:6288950/93‑1 (MQ=255)

                                           ttattaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAt                                              >  2:2465280/1‑93 (MQ=255)

                                           ttattaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAt                                              <  2:8815475/93‑1 (MQ=255)

                                            tattaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGCGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCaagt                                             >  2:5739650/1‑91 (MQ=255)

                                             attaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAttt                                            <  2:1587987/93‑1 (MQ=255)

                                              ttaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtttt                                           >  1:9126197/1‑93 (MQ=255)

                                              ttaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtttt                                           <  1:7107558/93‑1 (MQ=255)

                                              ttaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtttt                                           <  1:4724507/93‑1 (MQ=255)

                                                 cAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTc                                        >  1:9504129/1‑93 (MQ=255)

                                                 cAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTc                                        <  1:3533163/93‑1 (MQ=255)

                                                  agagGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCa                                       <  1:1963002/93‑1 (MQ=255)

                                                   ggggCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGACAGCACGCAATTTTTTCAt                                      >  2:7026325/3‑93 (MQ=255)

                                                   gagGCCCAAAAAACAAAAACCCCGTACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCAt                                      <  2:2606994/93‑1 (MQ=255)

                                                   gagGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCAt                                      <  1:3008886/93‑1 (MQ=255)

                                                   gagGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCAt                                      <  2:1299895/93‑1 (MQ=255)

                                                       cccAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcgcg                                  >  1:5156237/1‑93 (MQ=255)

                                                         cAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGtt                                <  2:2904240/93‑1 (MQ=255)

                                                          aaaaaaCAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTc                               <  2:3122171/93‑1 (MQ=255)

                                                          aaaaaaCAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTc                               >  1:8357674/1‑93 (MQ=255)

                                                           aaaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCt                              <  1:1868813/93‑1 (MQ=255)

                                                            aaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCtt                             >  2:481346/1‑93 (MQ=255)

                                                            aaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCtt                             <  1:1646944/93‑1 (MQ=255)

                                                            aaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCtt                             >  1:2835764/1‑93 (MQ=255)

                                                            aaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCtt                             <  2:8112835/93‑1 (MQ=255)

                                                             aaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCttt                            <  2:3614672/93‑1 (MQ=255)

                                                               acaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTc                          <  1:9855635/93‑1 (MQ=255)

                                                                caaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTcc                         >  1:4465451/1‑93 (MQ=255)

                                                                  aaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAg                       >  1:4212147/1‑93 (MQ=255)

                                                                  aaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAg                       <  2:7989109/93‑1 (MQ=255)

                                                                  aaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAg                       <  2:2304574/93‑1 (MQ=255)

                                                                    aacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCt                     <  2:6071988/93‑1 (MQ=255)

                                                                    aacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCt                     >  2:3760748/1‑93 (MQ=255)

                                                                    aacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCt                     <  2:1244886/93‑1 (MQ=255)

                                                                    aacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCt                     <  2:1667923/93‑1 (MQ=255)

                                                                     acCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTg                    <  1:1243760/93‑1 (MQ=255)

                                                                      ccccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTgg                   <  2:8345635/93‑1 (MQ=255)

                                                                       cccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGc                  >  1:10019604/1‑93 (MQ=255)

                                                                       cccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGc                  >  1:9390196/1‑93 (MQ=255)

                                                                       cccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGc                  >  1:8586374/1‑93 (MQ=255)

                                                                       cccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGc                  >  1:2480347/1‑93 (MQ=255)

                                                                         cGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGt                >  2:7225004/1‑93 (MQ=255)

                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               >  1:3662649/1‑93 (MQ=255)

                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               <  1:4572344/93‑1 (MQ=255)

                                                                           gACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTac              >  1:327701/1‑93 (MQ=255)

                                                                           gACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTac              <  1:5605203/93‑1 (MQ=255)

                                                                            aCTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTacc             >  1:7979160/1‑92 (MQ=255)

                                                                                tcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCt         >  1:6777253/1‑93 (MQ=255)

                                                                                 cATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTc        <  2:5452380/93‑1 (MQ=255)

                                                                                  aTCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCa       <  1:2222402/93‑1 (MQ=255)

                                                                                   tCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAg      <  1:255678/93‑1 (MQ=255)

                                                                                   tCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAg      >  1:7300484/1‑93 (MQ=255)

                                                                                   tCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAg      >  2:37229/1‑93 (MQ=255)

                                                                                    ccAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGc     >  1:4905735/1‑93 (MQ=255)

                                                                                    ccAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGc     <  2:9160590/93‑1 (MQ=255)

                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    <  1:3793914/93‑1 (MQ=255)

                                                                                       gggCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGa  >  1:5324360/1‑93 (MQ=255)

ACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGA  >  NC_000913/3597449‑3597628

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