BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,949,362	0	C	T	50.5%	88.4 / inf	406	N180N (AAC→AAT)	yecD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (109/92); new base T (100/105); total (209/197)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 416 total aligned reads displayed.

GCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCT  >  NC_000913/1949270‑1949454

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

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

  ccgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCa                                                                                            >  1:5598301/1‑93 (MQ=255)

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

   cgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCat                                                                                           >  1:5559000/1‑93 (MQ=255)

     cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCATCGCCGAGCAACACAATAACAGCATTAATCatat                                                                                         <  1:6142919/93‑1 (MQ=255)

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

      aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATc                                                                                        <  1:371456/93‑1 (MQ=255)

       aTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCt                                                                                       >  2:614976/1‑93 (MQ=255)

       aTGCCTGGGAACTTGGCTTTAATCTGGTGAGTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACCGTATCAACCGTATCt                                                                                       >  1:2289789/1‑93 (MQ=255)

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

         gCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAc                                                                                     <  1:2195741/93‑1 (MQ=255)

         gCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAc                                                                                     >  1:2527950/1‑93 (MQ=255)

          ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAcc                                                                                    <  2:5720411/93‑1 (MQ=255)

           cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAccc                                                                                   >  1:3544928/1‑93 (MQ=255)

           cTGGGAACTTGGATTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAccc                                                                                   <  2:5965553/93‑1 (MQ=255)

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

             gggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcgc                                                                                 <  1:5122342/93‑1 (MQ=255)

               gAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACACgcgc                                                                               >  1:2878886/1‑93 (MQ=255)

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

                aaCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCa                                                                              <  2:8461/93‑1 (MQ=255)

                 aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCAt                                                                             >  2:442391/1‑93 (MQ=255)

                 aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCAt                                                                             >  1:3466350/1‑93 (MQ=255)

                 aCTTGGCTTTAATCTGGCGATTGCCGAACACGCCTGTAGTGCCGACAGCGCCGAGCAGCAGAATAACAGTATCAACCATATCTACCCGCGCAt                                                                             >  1:622068/1‑93 (MQ=255)

                  cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATc                                                                            >  1:3624200/1‑93 (MQ=255)

                   ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATCACAGTATCAACCATATCTCCCCGGGCATCg                                                                           >  2:5550194/1‑93 (MQ=255)

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

                     ggCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGcc                                                                         <  1:5229274/93‑1 (MQ=255)

                     ggCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGcc                                                                         >  2:5251397/1‑93 (MQ=255)

                     ggCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGcc                                                                         <  1:1536782/93‑1 (MQ=255)

                       cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCg                                                                       >  1:2051711/1‑93 (MQ=255)

                         ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCgtg                                                                     >  1:4808502/1‑93 (MQ=255)

                         ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgtg                                                                     <  1:3654317/93‑1 (MQ=255)

                         ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgtg                                                                     <  2:5333102/93‑1 (MQ=255)

                           aaTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGtgtg                                                                   >  2:5091198/1‑93 (MQ=255)

                            aTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGc                                                                  >  1:3253624/1‑93 (MQ=255)

                             tCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCg                                                                 >  1:1421602/1‑93 (MQ=255)

                                ggTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                                              <  2:6077152/93‑1 (MQ=255)

                                 gTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGc                                                             <  1:4250679/93‑1 (MQ=255)

                                  tGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCg                                                            >  2:4005633/1‑93 (MQ=255)

                                   gATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGt                                                           <  2:2655693/93‑1 (MQ=255)

                                   gATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGt                                                           >  1:4660509/1‑93 (MQ=255)

                                    aTTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTg                                                          <  2:4258601/93‑1 (MQ=255)

                                    aTTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTg                                                          <  2:6248723/93‑1 (MQ=255)

                                    aTTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTg                                                          >  2:5408382/1‑93 (MQ=255)

                                       gCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGaa                                                       >  2:3189452/1‑93 (MQ=255)

                                       gCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGaa                                                       >  2:4195841/1‑93 (MQ=255)

                                        ccGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAag                                                      <  2:516862/93‑1 (MQ=255)

                                        ccGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAag                                                      <  1:4845841/93‑1 (MQ=255)

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

                                          gAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAagag                                                    >  1:724012/1‑93 (MQ=255)

                                           aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                                   >  1:2967166/1‑93 (MQ=255)

                                           aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgaga                                                   >  2:1710518/1‑93 (MQ=255)

                                            aGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGAt                                                  <  2:2854870/93‑1 (MQ=255)

                                             gACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATc                                                 >  2:3557701/1‑93 (MQ=255)

                                               cGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCt                                               <  2:5831779/93‑1 (MQ=255)

                                               cGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCt                                               <  1:1810297/93‑1 (MQ=255)

                                                  cTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCaa                                            <  1:948818/93‑1 (MQ=255)

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

                                                       gTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGt                                       >  1:3092318/1‑93 (MQ=255)

                                                          ccgccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTAt                                    <  1:4414876/93‑1 (MQ=255)

                                                             ccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAt                                 >  1:2506423/1‑93 (MQ=255)

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

                                                              cAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAtt                                <  2:3441736/93‑1 (MQ=255)

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

                                                                gcgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTa                              >  2:485232/1‑93 (MQ=255)

                                                                gcgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCACCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTa                              >  2:4631786/1‑93 (MQ=255)

                                                                   ccGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACAt                           <  1:1020908/93‑1 (MQ=255)

                                                                    cGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATc                          >  2:2392428/1‑93 (MQ=255)

                                                                    cGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATc                          <  2:3016768/93‑1 (MQ=255)

                                                                    cGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATc                          <  1:4188196/93‑1 (MQ=255)

                                                                     gAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCg                         >  1:3919769/1‑93 (MQ=255)

                                                                      agcagcACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCgg                        >  1:5723462/1‑93 (MQ=255)

                                                                      agcaacACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCgg                        <  1:569397/93‑1 (MQ=255)

                                                                      agcaacACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGAGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCgg                        >  1:6628/1‑93 (MQ=255)

                                                                       gcaacaCAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGt                       <  1:1233045/93‑1 (MQ=255)

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

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

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

                                                                          gcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTa                    <  2:4518383/93‑1 (MQ=255)

                                                                          acacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtg                    >  2:1349828/2‑92 (MQ=255)

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

                                                                            acaATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTAcc                  <  2:4949649/93‑1 (MQ=255)

                                                                            acaATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTAcc                  >  1:5934928/1‑93 (MQ=255)

                                                                            acaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGcc                  <  2:5598249/93‑1 (MQ=255)

                                                                              aaTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGc                >  2:3677576/1‑93 (MQ=255)

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

                                                                                tAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCaa              >  1:1569291/1‑93 (MQ=255)

                                                                                tAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGTGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCaa              >  2:4113056/1‑93 (MQ=255)

                                                                                 aaCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAAt             <  2:3936914/93‑1 (MQ=255)

                                                                                 aaCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAAt             <  1:312197/93‑1 (MQ=255)

                                                                                  aCAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATGGGTCTACCGCAATg            >  2:4709685/1‑93 (MQ=255)

                                                                                  aCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATg            <  2:2464449/93‑1 (MQ=255)

                                                                                    aGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGt          <  1:3786143/93‑1 (MQ=255)

                                                                                      tATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCg        >  2:2200574/1‑93 (MQ=255)

                                                                                      tATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCg        <  2:6283682/93‑1 (MQ=255)

                                                                                         cAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCAt     >  2:17840/1‑93 (MQ=255)

                                                                                          aaCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATc    <  1:5022912/93‑1 (MQ=255)

                                                                                          aaCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGCGCGCATc    >  1:4502467/1‑93 (MQ=255)

                                                                                           aCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATcc   >  2:1212185/1‑93 (MQ=255)

GCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCT  >  NC_000913/1949270‑1949454

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