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.4%	80.6 / inf	533	S167S (AGT→AGC)	yecD	isochorismatase family protein
Reads supporting (aligned to +/- strand): ref base T (125/149); new base C (125/132); total (251/282)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-07
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 549 total aligned reads displayed.

TGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAA  >  NC_000913/1949257‑1949456

   tGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCa                                                                                                  <  1:7036550/101‑1 (MQ=255)

     aaTCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTaat                                                                                                <  2:4499846/101‑2 (MQ=255)

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

       tCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCa                                                                                              >  1:1785952/1‑101 (MQ=255)

       tCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCa                                                                                              <  1:858594/101‑1 (MQ=255)

         cACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCata                                                                                            <  1:6731879/101‑1 (MQ=255)

         cACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCata                                                                                            >  1:1718188/1‑101 (MQ=255)

          aCCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCatat                                                                                           <  2:5428167/101‑1 (MQ=255)

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

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

             gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTa                                                                                        >  2:5594720/1‑101 (MQ=255)

             gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTa                                                                                        <  2:1833112/101‑1 (MQ=255)

               ccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAcc                                                                                      >  1:5606440/1‑101 (MQ=255)

               ccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAcc                                                                                      >  1:5801827/1‑101 (MQ=255)

               ccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAcc                                                                                      >  2:4753686/1‑101 (MQ=255)

                cgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAccc                                                                                     >  1:4062072/1‑101 (MQ=255)

                cgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAccc                                                                                     >  2:5857366/1‑101 (MQ=255)

                 gcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGAATCAACCATATCTACCcg                                                                                    <  2:3160116/101‑1 (MQ=255)

                   aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcgcg                                                                                  >  2:6472145/1‑101 (MQ=255)

                    aTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCgcgc                                                                                 >  2:3019527/1‑101 (MQ=255)

                     tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCa                                                                                >  1:4962771/1‑101 (MQ=255)

                      gCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTCCCCGCGCAt                                                                               <  2:2968743/101‑1 (MQ=255)

                      gCCTGGGAACTTGGCTTTAATCTGGGGATTGCCGAAGGCGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCAt                                                                               <  1:993277/101‑1 (MQ=255)

                       ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATc                                                                              >  1:6418209/1‑101 (MQ=255)

                       ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACATAATAACAGCATTAATCATATCTACCCGCGCATc                                                                              <  2:1036256/101‑1 (MQ=255)

                        cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCg                                                                             <  2:388076/101‑1 (MQ=255)

                          gggAACTTTGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGcc                                                                           <  2:4392271/101‑1 (MQ=255)

                          gggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGcc                                                                           <  2:4899973/101‑1 (MQ=255)

                           ggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGccc                                                                          <  1:1633969/101‑1 (MQ=255)

                           ggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGccc                                                                          <  2:1306825/101‑1 (MQ=255)

                            gAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCg                                                                         >  2:4016801/1‑101 (MQ=255)

                              aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgtg                                                                       >  1:6819722/1‑101 (MQ=255)

                               cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgtgt                                                                      <  2:5212956/101‑1 (MQ=255)

                               cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgtgt                                                                      >  1:347324/1‑101 (MQ=255)

                               cTTGGCTTTAATCTGGTGATTGCCGAAAACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATGTACCCGCGCATCGCCCgtgg                                                                      >  2:2009737/1‑100 (MQ=255)

                                 tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGc                                                                    >  2:6718794/1‑101 (MQ=255)

                                    cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTa                                                                 <  2:875547/101‑1 (MQ=255)

                                    cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTa                                                                 >  1:2192645/1‑101 (MQ=255)

                                    cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTa                                                                 <  1:6961274/101‑1 (MQ=255)

                                    cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTa                                                                 <  2:5819894/101‑1 (MQ=255)

                                    cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTa                                                                 >  1:4325958/1‑101 (MQ=255)

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

                                     tttAATCTGGTGATTGCCGAAGACGCCTGGAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                                                <  2:3802405/101‑1 (MQ=255)

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

                                          tCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTgg                                                           <  1:4724230/101‑1 (MQ=255)

                                          tCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCCCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTgg                                                           >  2:2489673/1‑101 (MQ=255)

                                            tGGTGTTTTCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCCACAGAATAACAGCATTAATCATCCCTACCCGCGGATCGCCCGTGTGCGTAGCGttaaa                                                         <  2:2187251/101‑5 (MQ=255)

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

                                             ggTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAag                                                        <  1:2887588/101‑1 (MQ=255)

                                             ggTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAag                                                        >  1:2287861/1‑101 (MQ=255)

                                             ggTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAag                                                        <  1:3811339/101‑1 (MQ=255)

                                             ggTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAag                                                        >  1:2451239/1‑101 (MQ=255)

                                               tGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAagag                                                      >  1:1312584/1‑101 (MQ=255)

                                                   tGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                                                  <  1:3633754/101‑1 (MQ=255)

                                                   tGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                                                  >  1:6575914/1‑101 (MQ=255)

                                                      cGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCa                                               <  1:2122631/101‑1 (MQ=255)

                                                      cGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCa                                               >  2:4181650/1‑101 (MQ=255)

                                                      cGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttc                                               >  1:1023671/1‑99 (MQ=255)

                                                        aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAc                                             <  1:1399122/101‑1 (MQ=255)

                                                         aGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttcacg                                            <  2:1751649/101‑6 (MQ=255)

                                                         aGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttcacg                                            >  2:328486/1‑96 (MQ=255)

                                                          gACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCAcgc                                           <  2:1900461/101‑1 (MQ=255)

                                                          caCGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcgc                                           <  2:7023059/100‑1 (MQ=255)

                                                            cGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGt                                         <  1:5861434/101‑1 (MQ=255)

                                                             gCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGtt                                        >  1:4191725/1‑101 (MQ=255)

                                                              ccTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTa                                       >  2:271184/1‑101 (MQ=255)

                                                              ccTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTa                                       <  1:2984360/101‑1 (MQ=255)

                                                               cTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTAt                                      >  1:1189540/1‑101 (MQ=255)

                                                               cTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTAt                                      >  2:2848256/1‑101 (MQ=255)

                                                               cTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTAt                                      <  1:4370598/101‑1 (MQ=255)

                                                                tGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATg                                     >  2:5048728/1‑101 (MQ=255)

                                                                     tGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTa                                <  2:1493895/101‑1 (MQ=255)

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

                                                                         gccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATc                            <  1:3592686/101‑1 (MQ=255)

                                                                         gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATc                            <  1:5313377/101‑1 (MQ=255)

                                                                            aGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGt                         >  1:2770165/1‑101 (MQ=255)

                                                                             gggcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTc                        >  1:3436850/3‑101 (MQ=255)

                                                                             gcgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTc                        <  1:502501/101‑1 (MQ=255)

                                                                               gcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTa                      <  2:1162241/101‑1 (MQ=255)

                                                                               gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtg                      >  2:945253/1‑100 (MQ=255)

                                                                                 cGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTAcc                    >  2:4271613/1‑101 (MQ=255)

                                                                                 cGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTAcc                    <  1:5151260/101‑1 (MQ=255)

                                                                                 cGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGcc                    <  2:6557692/101‑1 (MQ=255)

                                                                                  gAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCg                   >  2:5493594/1‑101 (MQ=255)

                                                                                     cagcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCaa                >  2:5232533/1‑101 (MQ=255)

                                                                                     cagcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCaa                <  2:6022725/101‑1 (MQ=255)

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

                                                                                      agcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAAt               <  1:4553901/101‑1 (MQ=255)

                                                                                      aacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAAt               >  1:3243962/3‑101 (MQ=255)

                                                                                        cacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATgg             >  1:6289448/1‑101 (MQ=255)

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

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

                                                                                          caATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTc           <  2:2426235/101‑1 (MQ=255)

                                                                                          caATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTc           <  2:4567726/101‑1 (MQ=255)

                                                                                               aCAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCCACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATc      >  1:6220599/1‑101 (MQ=255)

                                                                                               aCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTTGCATc      <  1:5956870/101‑1 (MQ=255)

                                                                                                 cgTATCAACCATATCTACCCGCGCACCGCCCGTGTGCGGAGCGTGGAAGAGATCCTCAACGCGTTATGATTGACATCGGTCTACCGCAATGGTCGCATCCt    <  1:4852425/100‑1 (MQ=255)

                                                                                                 aGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCt    <  2:6173637/101‑1 (MQ=255)

                                                                                                  gTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTa   <  1:1510368/101‑1 (MQ=255)

                                                                                                   tATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTaa  <  1:5406985/101‑1 (MQ=255)

TGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAA  >  NC_000913/1949257‑1949456

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