BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,949,356	0	T	C	48.2%	76.4 / inf	487	S178S (AGT→AGC)	yecD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (124/127); new base C (119/115); total (245/242)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-04
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 494 total aligned reads displayed.

TCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTC  >  NC_000913/1949264‑1949447

  cACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGTCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTAt                                                                                           >  1:4316749/1‑93 (MQ=255)

  cACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCAt                                                                                           >  1:5013223/1‑93 (MQ=255)

  cACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCAt                                                                                           >  1:4792164/1‑93 (MQ=255)

   aCCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATc                                                                                          <  1:3377525/93‑1 (MQ=255)

      gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAAc                                                                                       >  2:3143029/1‑93 (MQ=255)

       cccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAAcc                                                                                      <  2:4047460/93‑1 (MQ=255)

       cccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTaatc                                                                                      >  2:4586671/1‑91 (MQ=255)

        ccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCa                                                                                     >  2:671309/1‑93 (MQ=255)

         cgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCat                                                                                    <  1:5751208/93‑1 (MQ=255)

          gcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCata                                                                                   <  1:6211383/93‑1 (MQ=255)

           cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCatat                                                                                  <  1:2289145/93‑1 (MQ=255)

            aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATc                                                                                 >  1:128371/1‑93 (MQ=255)

            aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATc                                                                                 <  1:1392906/93‑1 (MQ=255)

            aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACATAATAACAGCATTAATCATATc                                                                                 <  1:2675941/93‑1 (MQ=255)

             aTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCt                                                                                <  2:5522635/93‑1 (MQ=255)

              tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTa                                                                               >  1:1971483/1‑93 (MQ=255)

              tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTa                                                                               <  1:4922276/93‑1 (MQ=255)

              tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTa                                                                               >  2:5772076/1‑93 (MQ=255)

                ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAcc                                                                             <  1:3955133/93‑1 (MQ=255)

                ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAcc                                                                             <  2:6511249/93‑1 (MQ=255)

                  tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcg                                                                           >  2:5893925/1‑93 (MQ=255)

                  tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                           >  1:1043871/1‑93 (MQ=255)

                   gggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcgc                                                                          >  1:1774234/1‑93 (MQ=255)

                   gggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcgc                                                                          >  1:3023696/1‑93 (MQ=255)

                    ggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcgcg                                                                         <  1:6103016/93‑1 (MQ=255)

                    ggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcgcg                                                                         >  2:2781513/1‑93 (MQ=255)

                    ggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcgcg                                                                         >  1:3684338/1‑93 (MQ=255)

                     gAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCgcgc                                                                        <  2:1810870/93‑1 (MQ=255)

                     gAACTTGGCTTTAATCTGGCGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCgcgc                                                                        >  1:44921/1‑93 (MQ=255)

                      aaCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCa                                                                       >  2:5113761/1‑93 (MQ=255)

                      aaCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCa                                                                       <  2:4213755/93‑1 (MQ=255)

                      aaCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCa                                                                       >  1:1502398/1‑93 (MQ=255)

                        cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATc                                                                     <  2:6075134/93‑1 (MQ=255)

                         ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCg                                                                    <  1:235735/93‑1 (MQ=255)

                         ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCg                                                                    >  1:1251019/1‑93 (MQ=255)

                         ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCg                                                                    >  2:563420/1‑93 (MQ=255)

                            gCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGccc                                                                 <  2:4530572/93‑1 (MQ=255)

                             cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCg                                                                <  2:4965522/93‑1 (MQ=255)

                              tttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgt                                                               <  1:1209076/93‑1 (MQ=255)

                               ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCgtg                                                              <  2:2606474/93‑1 (MQ=255)

                                tAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCgtgt                                                             <  1:6346581/93‑1 (MQ=255)

                                 aaTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGtgtg                                                            <  2:822744/93‑1 (MQ=255)

                                  aTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGc                                                           >  1:2425050/1‑93 (MQ=255)

                                   tCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCg                                                          <  1:4716856/93‑1 (MQ=255)

                                   tCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCg                                                          <  2:6369998/93‑1 (MQ=255)

                                     tGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTa                                                        >  2:379884/1‑93 (MQ=255)

                                      ggTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                                       >  2:234924/1‑93 (MQ=255)

                                      ggTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                                       >  2:4375473/1‑93 (MQ=255)

                                       gTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGc                                                      <  2:5986664/93‑1 (MQ=255)

                                         gATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGt                                                    >  2:3662814/1‑93 (MQ=255)

                                          aTTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTg                                                   >  1:6287379/1‑93 (MQ=255)

                                            tGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGa                                                 <  2:720608/93‑1 (MQ=255)

                                             gCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGaa                                                <  2:2909985/93‑1 (MQ=255)

                                               cGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAaga                                              >  1:3575183/1‑93 (MQ=255)

                                                gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAagag                                             <  1:5284256/93‑1 (MQ=255)

                                                 aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                            <  1:5458010/93‑1 (MQ=255)

                                                  aGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGCGGAAGAGAt                                           >  2:199683/1‑93 (MQ=255)

                                                  aGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGAt                                           >  1:4557030/1‑93 (MQ=255)

                                                    aCGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                                         >  1:2935633/1‑93 (MQ=255)

                                                    aCGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                                         >  2:3482393/1‑93 (MQ=255)

                                                    aCGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                                         <  1:5914785/93‑1 (MQ=255)

                                                       ccTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCa                                      >  1:1623553/1‑93 (MQ=255)

                                                        cTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCaa                                     >  1:3191451/1‑93 (MQ=255)

                                                        cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttca                                     <  2:3286125/93‑4 (MQ=255)

                                                         tGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAc                                    <  2:3935204/93‑1 (MQ=255)

                                                         tGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttcac                                    >  2:1564939/1‑89 (MQ=255)

                                                          gTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcg                                   >  2:3104767/1‑93 (MQ=255)

                                                           tAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcgc                                  <  2:1312270/93‑1 (MQ=255)

                                                           tAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCAcgc                                  <  2:4651465/93‑1 (MQ=255)

                                                            ggTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTCCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCAcgcg                                 <  2:1154167/92‑1 (MQ=255)

                                                            aGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcgcg                                 <  1:2185607/93‑1 (MQ=255)

                                                              tGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGtt                               <  1:3290998/93‑1 (MQ=255)

                                                               gccgccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTa                              <  2:5182011/93‑1 (MQ=255)

                                                                ccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTAt                             <  1:4158477/93‑1 (MQ=255)

                                                                   ccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAt                          <  1:34705/93‑1 (MQ=255)

                                                                   ccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAt                          <  1:2784036/93‑1 (MQ=255)

                                                                    cAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAtt                         >  2:2129486/1‑93 (MQ=255)

                                                                     aGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAttt                        <  1:890607/93‑1 (MQ=255)

                                                                     aGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAttt                        <  1:6569936/93‑1 (MQ=255)

                                                                     aGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAttt                        >  2:3876856/1‑93 (MQ=255)

                                                                       cgcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTAc                      <  2:1078861/93‑1 (MQ=255)

                                                                       cgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTAc                      <  1:723979/93‑1 (MQ=255)

                                                                        gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACa                     <  2:6118017/93‑1 (MQ=255)

                                                                           gAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCg                  >  2:173431/1‑93 (MQ=255)

                                                                           gAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCg                  <  2:6284663/93‑1 (MQ=255)

                                                                             gcagcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGt                >  1:1101579/1‑93 (MQ=255)

                                                                             gcaacaCAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGt                >  1:809724/1‑93 (MQ=255)

                                                                              cagcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTc               <  1:261182/93‑1 (MQ=255)

                                                                              caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTc               <  1:5141731/90‑1 (MQ=255)

                                                                               aacacaATAACAGCATTAATCATATCTACCCGCGCATCTCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCt              <  2:4464224/91‑1 (MQ=255)

                                                                               aacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCt              <  2:4559921/91‑1 (MQ=255)

                                                                                 cacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtgc            >  2:3417449/1‑91 (MQ=255)

                                                                                 cacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtgc            <  1:2764338/93‑3 (MQ=255)

                                                                                  acaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGcc           >  2:5342446/1‑93 (MQ=255)

                                                                                  acaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGcc           >  1:2539233/1‑93 (MQ=255)

                                                                                   caATAACAGTATCACCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCg          <  2:6187657/93‑1 (MQ=255)

                                                                                   caATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTACGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCg          >  2:3557517/1‑93 (MQ=255)

                                                                                     aTAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCa        >  1:987856/1‑93 (MQ=255)

                                                                                       aaCAGCATTAAGCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAAt      <  1:3839086/93‑1 (MQ=255)

                                                                                         cAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATgg    >  2:2430603/1‑93 (MQ=255)

TCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTC  >  NC_000913/1949264‑1949447

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