BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,490,327	T→C	P41P (CCT→CCC)	yhfC →	putative transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,490,327	0	T	C	99.4%	1537.4 / ‑2.5	479	P41P (CCT→CCC)	yhfC	putative transport
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (187/288); total (189/290)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

Only 100 of 489 total aligned reads displayed.

GATTAGCTTTCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCC  >  NC_000913/3490237‑3490416

gATTAGCTTTCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGt                                                                                         >  1:3962870/1‑93 (MQ=255)

gATTAGCTTTCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGt                                                                                         >  1:6344445/1‑93 (MQ=255)

   tAGCTTTCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTc                                                                                      >  2:209247/1‑93 (MQ=255)

     gCTTTCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCa                                                                                    <  2:4511769/93‑1 (MQ=255)

     gCTTTCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCa                                                                                    <  2:3471405/93‑1 (MQ=255)

       tttCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGt                                                                                  <  1:6668840/93‑1 (MQ=255)

        ttCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTa                                                                                 <  2:5155551/93‑1 (MQ=255)

        ttCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTa                                                                                 >  2:2625421/1‑93 (MQ=255)

        ttCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTa                                                                                 <  1:5679305/93‑1 (MQ=255)

         tctctcCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGGTGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCTCGTAt                                                                                >  1:172698/1‑93 (MQ=255)

              ccTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTa                                                                           <  1:7674548/93‑1 (MQ=255)

               cTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTaa                                                                          <  2:7606909/93‑1 (MQ=255)

               cTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTaa                                                                          <  1:1516464/93‑1 (MQ=255)

               cTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTaa                                                                          <  2:3869077/93‑1 (MQ=255)

               cTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTaa                                                                          >  2:1068619/1‑93 (MQ=255)

                tACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAAt                                                                         <  1:9136950/93‑1 (MQ=255)

                  cGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATac                                                                       >  2:5789882/1‑93 (MQ=255)

                  cGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATac                                                                       <  2:1576546/93‑1 (MQ=255)

                   gCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATacc                                                                      <  2:6817533/93‑1 (MQ=255)

                   gCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATacc                                                                      >  1:4244500/1‑93 (MQ=255)

                   gCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATacc                                                                      <  1:3770530/93‑1 (MQ=255)

                    cACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATacct                                                                     >  1:2660619/1‑93 (MQ=255)

                      ctgactgGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATaccttc                                                                   >  1:9591093/1‑93 (MQ=255)

                       tgactgGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATaccttca                                                                  <  2:3083547/93‑1 (MQ=255)

                       tgactgGTGCGTTGGTTATTGTCACCGGGATGGGGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATaccttca                                                                  <  1:9732084/93‑1 (MQ=255)

                        gactgGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATaccttcac                                                                 >  1:3544488/1‑93 (MQ=255)

                         actgGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATaccttcacc                                                                <  2:7771304/93‑1 (MQ=255)

                         actgGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATaccttcacc                                                                <  2:6172796/93‑1 (MQ=255)

                         actgGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATaccttcacc                                                                <  1:7313182/93‑1 (MQ=255)

                          cttgTGCGTTTGGTATTGGCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATaccttcacct                                                               <  1:5634915/90‑1 (MQ=255)

                            ggTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATaccttcaccttc                                                             <  2:4388184/93‑1 (MQ=255)

                            ggTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATaccttcaccttc                                                             <  1:7022867/93‑1 (MQ=255)

                             gTGCGTTGGTGATTGTCACCGGGATGGTGATGGGAAATAGCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTcc                                                            <  1:926245/93‑1 (MQ=255)

                               gCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTc                                                          <  2:8405791/93‑1 (MQ=255)

                                cGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCa                                                         <  1:8002053/93‑1 (MQ=255)

                                 gTTGGTTATTGTCACCGGGATGGTTATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCaa                                                        <  2:5225255/93‑1 (MQ=255)

                                 gTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCaa                                                        >  2:5890948/1‑93 (MQ=255)

                                  ttagTTATTGTCACCGGGATAGTGATGGAAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTACTACCTTCAACTTCCTCaaa                                                       <  1:8943254/90‑2 (MQ=255)

                                  ttGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAAc                                                       >  2:2338949/1‑93 (MQ=255)

                                  ttGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAAc                                                       <  2:9830659/93‑1 (MQ=255)

                                  ttGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAAc                                                       >  1:7537790/1‑93 (MQ=255)

                                  ttGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAAc                                                       <  1:2938881/93‑1 (MQ=255)

                                  gtGGTTATTGTCGCCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAAc                                                       <  1:5367384/92‑1 (MQ=255)

                                   tGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACg                                                      <  2:226487/93‑1 (MQ=255)

                                       tATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCgg                                                  <  1:6242376/93‑1 (MQ=255)

                                        aTTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGc                                                 >  2:2798040/1‑93 (MQ=255)

                                        aTTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGc                                                 >  2:8158507/1‑93 (MQ=255)

                                         ttGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCa                                                <  1:7373892/93‑1 (MQ=255)

                                         ttGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCa                                                >  2:9346870/1‑93 (MQ=255)

                                          tGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCAt                                               >  1:8513919/1‑93 (MQ=255)

                                            tCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCAttt                                             >  1:2086237/1‑93 (MQ=255)

                                            tCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCAttt                                             <  2:6551409/93‑1 (MQ=255)

                                               ccGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTaa                                          <  1:5039807/93‑1 (MQ=255)

                                                 gggATTGTTATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAAtc                                        <  1:3226822/93‑1 (MQ=255)

                                                 gggATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAAtc                                        >  2:8029323/1‑93 (MQ=255)

                                                 gggATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAAtc                                        <  2:8080122/93‑1 (MQ=255)

                                                  ggATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAAtct                                       <  1:336853/93‑1 (MQ=255)

                                                   gATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAAtctc                                      <  1:3691946/93‑1 (MQ=255)

                                                    aTGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATctct                                     <  1:5213631/93‑1 (MQ=255)

                                                     tGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTa                                    >  2:4280179/1‑93 (MQ=255)

                                                     tGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTa                                    <  1:2521730/93‑1 (MQ=255)

                                                     tGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTa                                    >  2:490492/1‑93 (MQ=255)

                                                     tGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTa                                    <  2:3400137/93‑1 (MQ=255)

                                                       gTGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATc                                  >  1:6605175/1‑93 (MQ=255)

                                                        tGATGGGAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCt                                 >  1:4853806/1‑93 (MQ=255)

                                                           agggAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTcc                              <  2:8801045/92‑1 (MQ=255)

                                                            gggAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCt                             <  1:4149117/93‑1 (MQ=255)

                                                             ggAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTc                            <  2:6952435/93‑1 (MQ=255)

                                                             ggAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTc                            >  1:10005473/1‑93 (MQ=255)

                                                             ggAAATATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTc                            >  1:1862390/1‑93 (MQ=255)

                                                               aaaTATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCaa                          <  2:406548/93‑1 (MQ=255)

                                                               aaaTATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCaa                          <  1:607285/93‑1 (MQ=255)

                                                                aaTATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAAc                         <  1:4570667/93‑1 (MQ=255)

                                                                aaTATCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAAc                         >  2:7380698/1‑93 (MQ=255)

                                                                 atatCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACg                        >  2:8220733/1‑93 (MQ=255)

                                                                 atatCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACg                        <  2:533624/93‑1 (MQ=255)

                                                                  tatCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGc                       >  1:9854454/1‑93 (MQ=255)

                                                                  tatCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGc                       <  1:9941704/93‑1 (MQ=255)

                                                                   atCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGcc                      <  2:9029327/93‑1 (MQ=255)

                                                                    tCGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCt                     <  1:4499706/93‑1 (MQ=255)

                                                                     cGCCGATTATTTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTg                    <  2:555855/93‑1 (MQ=255)

                                                                          attattTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGa               <  2:6397461/93‑1 (MQ=255)

                                                                          attattTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGa               <  2:2015111/93‑1 (MQ=255)

                                                                          attattTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGa               >  1:510009/1‑93 (MQ=255)

                                                                          attattTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGa               <  1:821725/93‑1 (MQ=255)

                                                                            tattTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATg             <  2:4686345/93‑1 (MQ=255)

                                                                            tattTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATg             <  1:5830711/93‑1 (MQ=255)

                                                                            tattTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATg             >  1:8795635/1‑93 (MQ=255)

                                                                             attTCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATgg            >  2:3776861/1‑93 (MQ=255)

                                                                              tttCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGa           <  1:219447/93‑1 (MQ=255)

                                                                                tCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGaaa         <  2:1690486/93‑1 (MQ=255)

                                                                                tCAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGaaa         >  1:2318427/1‑93 (MQ=255)

                                                                                 cAATCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAAt        <  2:9713271/93‑1 (MQ=255)

                                                                                   aTCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCg      <  2:7108862/93‑1 (MQ=255)

                                                                                   aTCTGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCg      <  1:7855709/93‑1 (MQ=255)

                                                                                      tGCCCGTTTCCAGTATGAGTAATCCCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTCCCTCAACGCCTGGCTGATGGAAATCGTcc   <  1:1113078/93‑1 (MQ=255)

                                                                                      tGCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTcc   <  2:1344902/93‑1 (MQ=255)

                                                                                       gCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTccc  <  2:324164/93‑1 (MQ=255)

                                                                                       gCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTccc  >  1:4298566/1‑93 (MQ=255)

                                                                                       gCCCGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTccc  <  2:8622382/93‑1 (MQ=255)

GATTAGCTTTCTCTCCTACGCACTGACTGGTGCGTTGGTTATTGTCACCGGGATGGTGATGGGAAATATCGCCGATTATTTCAATCTGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCC  >  NC_000913/3490237‑3490416

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