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%	93.7 / inf	1086	S167S (AGT→AGC)	yecD	isochorismatase family protein
Reads supporting (aligned to +/- strand): ref base T (272/287); new base C (250/270); total (529/557)
Fisher's exact test for biased strand distribution p-value = 8.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.15e-08
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 1131 total aligned reads displayed.

GTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAA  >  NC_000913/1949256‑1949456

  gtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGCCGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTAt                                                                                                    <  1:1950034/101‑1 (MQ=255)

  gtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTAt                                                                                                    <  2:5161871/101‑1 (MQ=255)

  gtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCAt                                                                                                    <  2:6020769/101‑1 (MQ=255)

    tGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTACTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCTAGCGCCGAGCAGCACAATAACAGTATCa                                                                                                  >  1:13006328/1‑101 (MQ=255)

     gAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCaa                                                                                                 >  1:12363185/1‑101 (MQ=255)

     gAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTaa                                                                                                 <  1:14371530/101‑1 (MQ=255)

      aaTCCACCGCCCGCAATGTCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATcnnn                                                                                                >  2:10452679/1‑98 (MQ=255)

      aaTCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAAc                                                                                                >  1:648863/1‑101 (MQ=255)

      aaTCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTaat                                                                                                >  1:5539884/1‑100 (MQ=255)

      aaTCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCAGCCAGCGCCGAGCAGCACAATAACAGTATCAAc                                                                                                >  2:5472292/1‑101 (MQ=255)

       aTCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAAcc                                                                                               <  1:5070673/101‑1 (MQ=255)

         ccACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCat                                                                                             <  1:4111735/101‑1 (MQ=255)

          nnncgcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCata                                                                                            <  2:8333485/98‑1 (MQ=255)

          cACCGCCCTCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCata                                                                                            >  1:2983917/1‑101 (MQ=255)

           cccgcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCatat                                                                                           <  2:1536403/100‑1 (MQ=255)

            ccgcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATc                                                                                          >  1:11288738/1‑101 (MQ=255)

             cgcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCt                                                                                         <  2:327236/101‑1 (MQ=255)

             cgcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCt                                                                                         <  2:2674/101‑1 (MQ=255)

              gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTa                                                                                        >  1:9909436/1‑101 (MQ=255)

               cccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAc                                                                                       <  1:9404880/101‑1 (MQ=255)

                nnncAATGCCTGGGGACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGCAGTGCCGCCAGCGCCGAGCACCACAATACCAGCATTAATCATATCTAcc                                                                                      <  2:11063401/98‑1 (MQ=255)

                ccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCACCACAATAACAGCATTAATCATATCTAcc                                                                                      >  2:4199600/1‑101 (MQ=255)

                cagcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCtnnn                                                                                      >  2:10833274/3‑98 (MQ=255)

                 cgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAccc                                                                                     <  2:6259381/101‑1 (MQ=255)

                 cgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAccc                                                                                     <  1:5895765/101‑1 (MQ=255)

                  gcAATGCCTGGGACCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                                    <  1:11868334/101‑1 (MQ=255)

                  gcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                                    >  1:5291356/1‑101 (MQ=255)

                   nnntGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcgc                                                                                   <  2:7653498/98‑1 (MQ=255)

                    aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcgcg                                                                                  >  1:8495616/1‑101 (MQ=255)

                       gCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCAt                                                                               <  2:4617375/101‑1 (MQ=255)

                       gCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATAAACCATATCTACCCGCGCAt                                                                               >  1:7707012/1‑101 (MQ=255)

                          tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGc                                                                            >  1:14608282/1‑101 (MQ=255)

                            nnnaCTTGTCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCTGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGccc                                                                          <  2:12464031/98‑1 (MQ=255)

                             gAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGcnnn                                                                         >  2:12866334/1‑98 (MQ=255)

                                cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCgtgt                                                                      >  1:1213756/1‑101 (MQ=255)

                                 ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGtgtg                                                                     <  2:5678960/101‑1 (MQ=255)

                                 ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGtgtg                                                                     <  1:1659639/101‑1 (MQ=255)

                                 ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGtgtg                                                                     <  2:2351319/101‑1 (MQ=255)

                                  tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGc                                                                    <  1:2075139/101‑1 (MQ=255)

                                    gCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATCACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGt                                                                  <  1:14004891/101‑1 (MQ=255)

                                      tttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                                                >  1:4268401/1‑101 (MQ=255)

                                      tttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                                                >  1:14482352/1‑101 (MQ=255)

                                      tttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                                                <  1:2531481/101‑1 (MQ=255)

                                      nnnaaTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                                                <  2:8168551/98‑1 (MQ=255)

                                       ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGc                                                               <  1:10578680/101‑1 (MQ=255)

                                        tAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACATGATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTannn                                                              >  2:8993785/1‑98 (MQ=255)

                                        nnntCTGGTGATTGCCGAAGACGCCTGTAGTTCCGCTAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCTGTGTGCGTAGCg                                                              <  2:9659631/98‑1 (MQ=255)

                                        nnntCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCg                                                              <  2:14527876/98‑1 (MQ=255)

                                         aaTCTGGTGATTGCCGAATACGCCTGTAGTTCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCATAGCGt                                                             <  1:11573840/101‑1 (MQ=255)

                                         aaTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGt                                                             <  1:8264881/101‑1 (MQ=255)

                                            nnngTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGa                                                          <  2:9342874/98‑1 (MQ=255)

                                            cTGATGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGa                                                          >  1:10406914/1‑101 (MQ=255)

                                             tGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGaa                                                         <  1:10000190/101‑1 (MQ=255)

                                             tGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGaa                                                         >  1:6216517/1‑101 (MQ=255)

                                             tGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGaa                                                         <  1:2764231/101‑1 (MQ=255)

                                                   ttGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATc                                                   >  1:3842841/1‑101 (MQ=255)

                                                     gCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGAtnnn                                                 >  2:8501054/1‑98 (MQ=255)

                                                        gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCaa                                              >  1:9727676/1‑101 (MQ=255)

                                                        gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCaa                                              >  1:13367104/1‑101 (MQ=255)

                                                        gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                                     <  1:6794392/94‑1 (MQ=255)

                                                        gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                                     >  2:6794392/1‑94 (MQ=255)

                                                          aGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcg                                            >  1:2280627/1‑101 (MQ=255)

                                                          aGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcg                                            <  2:785261/101‑1 (MQ=255)

                                                              gCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGtt                                        >  1:11032205/1‑101 (MQ=255)

                                                              gCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGtt                                        <  1:8858468/101‑1 (MQ=255)

                                                               ccTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGGTa                                       >  2:20325/1‑101 (MQ=255)

                                                               ccTGCAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTACCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTa                                       <  1:14778997/101‑1 (MQ=255)

                                                                cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTCCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTAt                                      <  2:6482977/101‑1 (MQ=255)

                                                                cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTAt                                      <  1:9087631/101‑1 (MQ=255)

                                                                  atAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTannn                                    >  2:11722236/2‑98 (MQ=255)

                                                                    aGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAtt                                  <  2:3784791/101‑1 (MQ=255)

                                                                     nnnccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAttt                                 <  2:8772138/98‑1 (MQ=255)

                                                                     gTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAttt                                 >  1:3429007/1‑101 (MQ=255)

                                                                       gccgccAGCGCCGTGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAttnnn                               >  2:14263241/1‑98 (MQ=255)

                                                                       gccgccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTAc                               >  2:4020073/1‑101 (MQ=255)

                                                                       gccgccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTAc                               >  2:3525141/1‑101 (MQ=255)

                                                                        nnnccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACa                              <  2:12271277/98‑1 (MQ=255)

                                                                        cccccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACa                              <  2:3155042/98‑1 (MQ=255)

                                                                          gccAGCGCCGAGCAGCACAATAACAGTATCAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATc                            <  1:719225/101‑1 (MQ=255)

                                                                           ccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCg                           <  2:7019576/101‑1 (MQ=255)

                                                                            cAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCgg                          <  1:14190434/101‑1 (MQ=255)

                                                                               caccGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCggnnn                       >  2:11375171/3‑98 (MQ=255)

                                                                                gcCGAGCAGCACAATAACAGTATCAACCATATCTACCAGCGAATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATCTACATCGGtnnn                      >  2:13346486/1‑98 (MQ=255)

                                                                                  cGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGCTATGATTTACATCGGTCTAcc                    >  2:4876715/1‑101 (MQ=255)

                                                                                  cGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGcc                    >  1:4516054/1‑101 (MQ=255)

                                                                                   gAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCg                   <  1:13729173/101‑1 (MQ=255)

                                                                                       agcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAAt               >  2:730213/1‑101 (MQ=255)

                                                                                         cacaATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATgg             >  1:8029573/1‑101 (MQ=255)

                                                                                          acaATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGt            >  1:4778162/1‑101 (MQ=255)

                                                                                          acaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGt            >  2:1708255/1‑101 (MQ=255)

                                                                                           caATAACAGTATCAACCATATCTCCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTc           <  1:1019749/101‑1 (MQ=255)

                                                                                           caATAACAGCATTGATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATgnnn           >  2:13831153/1‑98 (MQ=255)

                                                                                            aaTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATggnnn          >  2:11959107/1‑98 (MQ=255)

                                                                                            aaTAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCg          <  1:10833274/101‑1 (MQ=255)

                                                                                             aTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGc         >  1:12664320/1‑101 (MQ=255)

                                                                                              tAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCa        >  2:2947933/1‑101 (MQ=255)

                                                                                               aaCAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATAGGTCTACCGCAATGGTCgnnn       >  2:13146303/1‑98 (MQ=255)

                                                                                                aCAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATc      >  1:3157798/1‑101 (MQ=255)

                                                                                                aCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATc      >  1:3686161/1‑101 (MQ=255)

                                                                                                aCAGCATTAATCATATCTACCCGCGCATCGCCCGTGGGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATc      >  1:6045820/1‑101 (MQ=255)

GTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAA  >  NC_000913/1949256‑1949456

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 17 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`