BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,182,712	G→A	A55T (GCC→ACC)	zupT →	zinc transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,182,712	0	G	A	97.7%	442.7 / ‑0.2	129	A55T (GCC→ACC)	zupT	zinc transporter
Reads supporting (aligned to +/- strand): ref base G (1/2); new base A (56/69); total (58/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 129 total aligned reads displayed.

TTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCAC  >  NC_000913/3182617‑3182808

ttCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGc                                                                                              <  1:1305039/100‑1 (MQ=255)

 tctcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCa                                                                                             >  1:1023889/1‑100 (MQ=255)

  ctcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCAc                                                                                            >  1:808555/1‑100 (MQ=255)

  ctcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCAc                                                                                            >  2:3114015/1‑100 (MQ=255)

   tcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACt                                                                                           <  2:1031765/100‑1 (MQ=255)

   tcGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACt                                                                                           <  1:2600836/100‑1 (MQ=255)

     gggCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTag                                                                                         >  1:3162142/1‑100 (MQ=255)

      ggCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagc                                                                                        <  1:1593851/100‑1 (MQ=255)

      ggCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagc                                                                                        >  2:2732402/1‑100 (MQ=255)

      ggCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagc                                                                                        >  2:1040449/1‑100 (MQ=255)

      ggCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagc                                                                                        <  1:2673648/100‑1 (MQ=255)

      ggCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagc                                                                                        <  2:2623334/100‑1 (MQ=255)

       gCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagca                                                                                       <  2:3374312/100‑1 (MQ=255)

        cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagcag                                                                                      >  1:1185793/1‑100 (MQ=255)

        cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagcag                                                                                      >  1:385816/1‑100 (MQ=255)

        cAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGACGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagcag                                                                                      <  1:2848613/100‑1 (MQ=255)

          aaaaCCCTCGAACCTCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATTATGTTGCTCATCTCAGTAATGGAAATGCTTCCTACCGCACTagcagc                                                                                     <  1:2883323/99‑1 (MQ=255)

              cccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAg                                                                                >  1:2392223/1‑100 (MQ=255)

              cccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTagatcggaaga                                                                                >  1:3277877/1‑91 (MQ=255)

              cccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAg                                                                                <  1:3322037/100‑1 (MQ=255)

              cccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAg                                                                                >  1:2399264/1‑100 (MQ=255)

               ccTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAgg                                                                               <  2:2458871/100‑1 (MQ=255)

                cTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGa                                                                              <  2:385815/100‑1 (MQ=255)

                cTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGa                                                                              >  1:1961838/1‑100 (MQ=255)

                 tCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGaa                                                                             <  1:2274321/100‑1 (MQ=255)

                  cGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAAt                                                                            <  2:2712770/100‑1 (MQ=255)

                  cGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAAt                                                                            >  2:3197836/1‑100 (MQ=255)

                  cGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAAt                                                                            <  1:2438726/100‑1 (MQ=255)

                     aCCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTc                                                                         <  2:808555/100‑1 (MQ=255)

                      ccGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCg                                                                        <  1:2019317/100‑1 (MQ=255)

                         cTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCt                                                                     <  2:3307459/100‑1 (MQ=255)

                          ttACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCtg                                                                    <  2:1228223/100‑1 (MQ=255)

                            aCTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCtgtg                                                                  <  1:658274/100‑1 (MQ=255)

                            aCTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCtgtg                                                                  <  2:2932418/100‑1 (MQ=255)

                               ggCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTg                                                               >  2:493828/1‑100 (MQ=255)

                                gCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTgg                                                              <  2:1286350/100‑1 (MQ=255)

                                gCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTgg                                                              <  2:634500/100‑1 (MQ=255)

                                 cGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTggg                                                             >  2:3342890/1‑100 (MQ=255)

                                  gTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGt                                                            >  2:2961489/1‑100 (MQ=255)

                                  gTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGt                                                            <  2:1763510/100‑1 (MQ=255)

                                  gTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGt                                                            <  2:2221720/100‑1 (MQ=255)

                                      tCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATg                                                        <  1:2466821/100‑1 (MQ=255)

                                       cGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATgg                                                       <  1:2248659/100‑1 (MQ=255)

                                        gCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGa                                                      >  1:272978/1‑100 (MQ=255)

                                        gCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGa                                                      >  1:565541/1‑100 (MQ=255)

                                         cTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGaa                                                     >  1:1371227/1‑100 (MQ=255)

                                          tAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAAt                                                    <  1:2376240/100‑1 (MQ=255)

                                           aGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATg                                                   >  1:3351019/1‑100 (MQ=255)

                                            ggTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGt                                                  >  1:1666142/1‑100 (MQ=255)

                                            ggTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGt                                                  >  1:3145522/1‑100 (MQ=255)

                                             gTTTTGCGGCGAGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGtt                                                 >  2:1790717/1‑100 (MQ=255)

                                              ttttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGttt                                                >  1:437084/1‑100 (MQ=255)

                                               tttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGtttn                                               <  2:433454/100‑2 (MQ=255)

                                               gttGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTa                                               >  2:2142274/2‑100 (MQ=255)

                                                   cggcggTGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCtt                                           >  2:1746001/1‑100 (MQ=255)

                                                      tccgatctaTGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCgg                                        <  1:1948394/92‑1 (MQ=255)

                                                       ggggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGt                                       <  2:1023888/100‑1 (MQ=255)

                                                        gggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTc                                      <  1:496157/100‑1 (MQ=255)

                                                         ggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCt                                     >  2:1741427/1‑100 (MQ=255)

                                                         ggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCt                                     >  2:494825/1‑100 (MQ=255)

                                                          gATCATGTTGCTCATCTCATTAATGGAAATGCTGCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCgg                                        <  1:1299196/96‑1 (MQ=255)

                                                              aTGTTGCTCATCTCATTAATGGAAATGCTTCCTCCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCggagatcgga                                >  2:1948390/1‑92 (MQ=255)

                                                               tGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAGTGTTTATCTTCGGTCTGCTTgg                               >  1:1993678/1‑100 (MQ=255)

                                                               tGTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTgg                               <  2:731583/100‑1 (MQ=255)

                                                               tGTTGCTCATCTCATTAATGGAAATGATTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTgg                               >  1:2073026/1‑100 (MQ=255)

                                                                gTTGCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGc                              <  1:3114027/100‑1 (MQ=255)

                                                                   gCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAGTGTTTATCTTCGGTCTGCTTGGCTAt                           <  2:1993674/100‑1 (MQ=255)

                                                                   gCTCATCTCATTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAt                           >  1:417148/1‑100 (MQ=255)

                                                                      catctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAtttt                        <  1:3197849/100‑1 (MQ=255)

                                                                      catctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAtttt                        >  1:1934359/1‑100 (MQ=255)

                                                                      catctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAtttt                        >  1:2005163/1‑100 (MQ=255)

                                                                       gatctatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTg                       <  2:2992751/95‑1 (MQ=255)

                                                                       atctcatTAATGGAAATGCTTCCTACCACACTAGCAGCTGAAGAAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTg                       <  1:1006873/100‑1 (MQ=255)

                                                                        tctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTgg                      <  2:949392/100‑1 (MQ=255)

                                                                        tctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTgg                      <  2:3162130/100‑1 (MQ=255)

                                                                        tctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTgg                      >  2:2076657/1‑100 (MQ=255)

                                                                        tctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTgg                      <  1:655742/100‑1 (MQ=255)

                                                                        tctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTgg                      <  1:655961/100‑1 (MQ=255)

                                                                        tctcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTgg                      <  1:2961501/100‑1 (MQ=255)

                                                                          tcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGcc                    <  2:2452878/100‑1 (MQ=255)

                                                                          tcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGcc                    <  1:466414/100‑1 (MQ=255)

                                                                          tcatTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGcc                    <  2:1371224/100‑1 (MQ=255)

                                                                           natTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTgagat                   >  2:1006872/2‑96 (MQ=255)

                                                                           catTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCt                   >  2:1343683/1‑100 (MQ=255)

                                                                           catTAATGGAAATGCTTCCTACCGCACTAGCAACTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCt                   >  1:3228547/1‑100 (MQ=255)

                                                                            atTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTgagatc                  >  1:2992763/1‑95 (MQ=255)

                                                                            atTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTg                  >  2:2662486/1‑100 (MQ=255)

                                                                            atTAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTCTGGCCTg                  <  1:1790721/100‑1 (MQ=255)

                                                                             ttAATGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTg                  >  2:629472/1‑99 (MQ=255)

                                                                                aTGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGAcc              <  2:1325438/100‑1 (MQ=255)

                                                                                aTGGAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGAcc              <  2:1955971/100‑1 (MQ=255)

                                                                                  ggAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGc            >  2:1434402/1‑100 (MQ=255)

                                                                                   gAAATGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCa           >  1:1094358/1‑100 (MQ=255)

                                                                                     aaTGCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATg         <  2:1185790/100‑1 (MQ=255)

                                                                                        gCTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTg      <  1:1507200/100‑1 (MQ=255)

                                                                                         cTTCCTACCGCACTCGCCGCCGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATGTTGGCCTGTACCGCATGTTGc     >  2:1416789/1‑100 (MQ=255)

                                                                                         cTTCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGc     <  2:640999/100‑1 (MQ=255)

                                                                                           tCCTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCca   >  2:1000177/1‑100 (MQ=255)

                                                                                            ccTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGAGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCcac  <  2:2392216/100‑1 (MQ=255)

                                                                                            ccTACCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCcac  >  2:31078/1‑100 (MQ=255)

TTCTCGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCAC  >  NC_000913/3182617‑3182808

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