BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,949,359	0	C	T	47.2%	165.1 / inf	530	I168I (ATC→ATT)	yecD	isochorismatase family protein
Reads supporting (aligned to +/- strand): ref base C (129/151); new base T (121/129); total (250/280)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-11
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 550 total aligned reads displayed.

TTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATG  >  NC_000913/1949259‑1949459

ttGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATc                                                                                                      <  1:69036/101‑1 (MQ=255)

  gAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCaa                                                                                                    <  1:4360640/101‑1 (MQ=255)

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

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

        aCCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCatat                                                                                              <  1:722748/101‑1 (MQ=255)

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

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

         ccgcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATc                                                                                             <  2:4175946/101‑1 (MQ=255)

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

           gcccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGAATCAACCATATCTa                                                                                           >  2:4741300/1‑101 (MQ=255)

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

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

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

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

              cgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAccc                                                                                        <  2:2185423/101‑1 (MQ=255)

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

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

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

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

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

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

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

                     ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATc                                                                                 <  1:4146225/101‑1 (MQ=255)

                      cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCg                                                                                >  2:529972/1‑101 (MQ=255)

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

                      cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCg                                                                                >  1:6086073/1‑101 (MQ=255)

                        gggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGcc                                                                              >  2:1161841/1‑101 (MQ=255)

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

                           aaCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCgt                                                                           >  1:5402870/1‑101 (MQ=255)

                           aaCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgt                                                                           >  1:3589597/1‑101 (MQ=255)

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

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

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

                                ggCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCg                                                                      >  1:4961586/1‑101 (MQ=255)

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

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

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

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

                                    ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGc                                                                  >  1:2748603/1‑101 (MQ=255)

                                     tAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCg                                                                 <  2:5433902/101‑1 (MQ=255)

                                      aaTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGt                                                                >  1:2082116/1‑101 (MQ=255)

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

                                         cTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGa                                                             >  1:2562648/1‑101 (MQ=255)

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

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

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

                                              gATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                                        <  1:1547546/101‑1 (MQ=255)

                                              gATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                                        >  2:1692301/1‑101 (MQ=255)

                                                ttGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATc                                                      <  1:2947144/101‑1 (MQ=255)

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

                                                 tGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCg                                                                                <  2:2413502/74‑1 (MQ=255)

                                                  gCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCt                                                    <  2:3035069/101‑1 (MQ=255)

                                                  gCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCt                                                    <  1:980126/101‑1 (MQ=255)

                                                   ccGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCtt                                                   >  2:5500144/1‑100 (MQ=255)

                                                      aaGACGCTTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAc                                                <  2:5159910/101‑1 (MQ=255)

                                                      aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttcac                                                <  2:6298931/101‑5 (MQ=255)

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

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

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

                                                           gCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTTCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGCTCCTCAACGCGtt                                           <  2:14353/101‑1 (MQ=255)

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

                                                             cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGttcg                                         <  2:1019883/101‑3 (MQ=255)

                                                             cTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTAt                                         >  2:5602259/1‑101 (MQ=255)

                                                                 aGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAtt                                     >  1:3415005/1‑101 (MQ=255)

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

                                                                   tGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTa                                   <  1:227942/101‑1 (MQ=255)

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

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

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

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

                                                                             gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGTTCCTTCACGCGTTATGATTTACATCGGTCtg                         >  1:1352578/1‑100 (MQ=255)

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

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

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

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

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

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

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

                                                                                 agcagcACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGc                     <  1:5796083/101‑1 (MQ=255)

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

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

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

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

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

                                                                                         aaTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCg             >  1:5106992/1‑101 (MQ=255)

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

                                                                                           tAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCa           >  1:6973282/1‑101 (MQ=255)

                                                                                            aaCAGTATCAACCATATCTTCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAt          <  1:5871693/101‑1 (MQ=255)

                                                                                            aaCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAt          <  1:4024411/101‑1 (MQ=255)

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

                                                                                             aCAGTATCAACCATATCTACCCGCACATCGCCCGTGTGCGTAGCGTGGAACAGATCCTCAACGCGTTATGATTTACATCGGCCTAATGCAATGGTCGCATc         >  2:5222873/1‑101 (MQ=255)

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

                                                                                             aCAGCATTAATCATATCTACCCGCGCATCGCCCGGGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATc         <  2:1898580/101‑1 (MQ=255)

                                                                                              cAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCACAACGCGTTATGATTTACATCGGTCTACCGCAATGGTAGCATct        >  1:380928/1‑100 (MQ=255)

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

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

                                                                                               aGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCt       >  1:425290/1‑101 (MQ=255)

                                                                                                   tCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAAt   <  2:5858131/101‑1 (MQ=255)

                                                                                                   tCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAAt   <  1:1147034/101‑1 (MQ=255)

                                                                                                   tCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAAt   <  1:1004725/101‑1 (MQ=255)

TTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATG  >  NC_000913/1949259‑1949459

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