BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	580,771	0	C	T	48.0%	104.2 / inf	1099	G65G (GGC→GGT)	ybcX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (328/241); new base T (242/284); total (574/525)
Fisher's exact test for biased strand distribution p-value = 1.35e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 1155 total aligned reads displayed.

GCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAA  >  NC_000913/580676‑580866

  ggATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCTGGCCATCATGAATGCGATGGGCAg                                                                                               >  2:683030/1‑96 (MQ=255)

  ggATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAg                                                                                               <  2:13611593/96‑1 (MQ=255)

   gATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGc                                                                                              <  1:10842187/96‑1 (MQ=255)

    aTGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCg                                                                                             >  1:2062849/1‑96 (MQ=255)

     cgatAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGa                                                                                            >  2:765887/4‑96 (MQ=255)

     cgatAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGAGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGa                                                                                            >  2:10293032/4‑96 (MQ=255)

      gCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGAc                                                                                           <  2:9667638/96‑1 (MQ=255)

       cTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACt                                                                                          <  1:4906907/96‑1 (MQ=255)

       atAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGatt                                                                                          <  1:9398603/95‑3 (MQ=255)

       atAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGatt                                                                                          <  2:570247/95‑3 (MQ=255)

        tAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTa                                                                                         <  2:9816188/96‑1 (MQ=255)

        tAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTa                                                                                         <  2:2566611/96‑1 (MQ=255)

          agtACTCTCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACTCTGCCCTTTCACCGGGCCTTCATGAATGCGATGGGCAGCGACTACa                                                                                       <  1:2460593/94‑1 (MQ=255)

          aTTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCCGCCCATCATGAATGCGATGGGCAGCGACTACa                                                                                       >  2:11644614/1‑96 (MQ=255)

          aTTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACa                                                                                       >  2:4917460/1‑96 (MQ=255)

           ttACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACAt                                                                                      >  2:2014455/1‑96 (MQ=255)

           ttACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAt                                                                                      >  1:4481676/1‑96 (MQ=255)

           ttACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAt                                                                                      <  1:7515264/96‑1 (MQ=255)

            tACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAta                                                                                     >  1:3586879/1‑95 (MQ=255)

            tACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAta                                                                                     <  2:11999994/96‑2 (MQ=255)

            tACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACAta                                                                                     <  2:2122456/96‑2 (MQ=255)

             aCTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATcc                                                                                    <  2:7340514/96‑1 (MQ=255)

             aCTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATcc                                                                                    <  2:8033290/96‑1 (MQ=255)

              cTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTAAATCCg                                                                                   >  1:3095988/1‑96 (MQ=255)

              cTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACg                                                                                   <  2:7114063/96‑1 (MQ=255)

               tATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGt                                                                                  >  1:8958202/1‑96 (MQ=255)

                aTCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTg                                                                                 >  2:9215128/1‑96 (MQ=255)

                aTCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTg                                                                                 <  2:5298989/96‑1 (MQ=255)

                  ctcCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAg                                                                               <  1:13821104/96‑1 (MQ=255)

                   tcCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAgg                                                                              >  2:9053380/1‑96 (MQ=255)

                     ccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTg                                                                            >  1:2887005/1‑96 (MQ=255)

                     ccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTg                                                                            >  1:5994242/1‑96 (MQ=255)

                       gAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGaa                                                                          >  1:3405770/1‑96 (MQ=255)

                        aaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCaac                                                                         <  1:655416/96‑2 (MQ=255)

                        aaaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAt                                                                         >  1:5137787/1‑96 (MQ=255)

                        aaaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAt                                                                         <  1:10001121/96‑1 (MQ=255)

                          aGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAtgt                                                                       <  2:13356265/96‑1 (MQ=255)

                            aaTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTGCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGGGAGGTGAATgtgg                                                                     >  1:731266/1‑96 (MQ=255)

                             aTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATgtggt                                                                    >  2:4552930/1‑96 (MQ=255)

                                ggcctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTaa                                                                 >  2:555123/5‑96 (MQ=255)

                                  cctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTaatt                                                               <  1:10564612/94‑3 (MQ=255)

                                  cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGt                                                               >  1:8678199/1‑96 (MQ=255)

                                   ctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTc                                                              >  1:1880085/2‑96 (MQ=255)

                                    tACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGGGAGGTGAATGTGGTGAGGTCt                                                             >  1:6921317/1‑96 (MQ=255)

                                       cAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGcc                                                          >  1:845570/1‑96 (MQ=255)

                                       cAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGcc                                                          <  1:14066094/96‑1 (MQ=255)

                                        aGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGcac                                                         <  1:11517192/96‑3 (MQ=255)

                                        aGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCAGTCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGccc                                                         <  1:14283961/96‑1 (MQ=255)

                                         ggAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCg                                                        >  2:1420594/1‑96 (MQ=255)

                                          gAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgt                                                       >  1:12116568/1‑96 (MQ=255)

                                          gAAGGGCGCTGGAAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgt                                                       <  2:1523550/96‑1 (MQ=255)

                                           aaGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtg                                                      >  1:3812285/1‑96 (MQ=255)

                                             ggggGCTGGTAAACCCTGCCCTTTCAGCGGGTCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCAGTGTc                                                    <  2:3435454/96‑1 (MQ=255)

                                              ggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCg                                                   <  2:11359266/96‑1 (MQ=255)

                                              ggCGCTGGAAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCg                                                   <  2:10750777/96‑1 (MQ=255)

                                               gcgcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCgg                                                  >  1:5651314/1‑96 (MQ=255)

                                                 gcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGtt                                                >  2:6164220/1‑96 (MQ=255)

                                                 gcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGtt                                                >  1:5378473/1‑96 (MQ=255)

                                                  cTGGGAAACACTGCCCTGTCAGCGGGCCATCATTAATGCGATGGGCAGCTACTACATCCGGGATGTGATATTAGTGAAGTCTGCCAGTGTCgggtg                                               >  1:2662772/1‑93 (MQ=255)

                                                   tGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTaa                                              <  1:11334791/96‑2 (MQ=255)

                                                   tGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTAt                                              >  1:4650698/1‑96 (MQ=255)

                                                   tGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGCCTGCCCGTGTCGGTTAt                                              >  2:6279300/1‑96 (MQ=255)

                                                   tGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGAGGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCAGTGTAGGtttg                                              >  1:4117255/1‑94 (MQ=255)

                                                     ngAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGGCGGGTAtgc                                            >  1:6212245/2‑94 (MQ=255)

                                                     ggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTc                                            <  2:12221955/96‑1 (MQ=255)

                                                        aaCACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaa                                         <  2:13756998/96‑1 (MQ=255)

                                                         acgttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCaaaa                                        <  1:8150130/92‑1 (MQ=255)

                                                         acacTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaa                                        >  2:12666211/1‑96 (MQ=255)

                                                           gttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAt                                      >  2:8932445/3‑96 (MQ=255)

                                                            ttGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtg                                     >  2:11142681/2‑96 (MQ=255)

                                                               cccTTTCAGCGGTCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAGTGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgctg                                  >  1:12954934/1‑96 (MQ=255)

                                                               cccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgctg                                  >  2:2320567/1‑96 (MQ=255)

                                                                  tttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGttgctg                               <  1:12335106/96‑1 (MQ=255)

                                                                   ttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTgg                              >  1:2246044/1‑96 (MQ=255)

                                                                    tCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTggg                             <  2:8258114/96‑1 (MQ=255)

                                                                    tCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTggg                             <  2:428043/96‑1 (MQ=255)

                                                                       gCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtgt                          >  2:406243/1‑96 (MQ=255)

                                                                        cGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGtt                         >  2:1741139/1‑96 (MQ=255)

                                                                        cGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGtt                         >  1:12667993/1‑96 (MQ=255)

                                                                         gggCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGttt                        >  1:993959/1‑96 (MQ=255)

                                                                         gggCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGttt                        >  1:7736595/1‑96 (MQ=255)

                                                                          ggCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTa                       <  2:9435676/96‑1 (MQ=255)

                                                                            cgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATg                     <  1:13308199/94‑1 (MQ=255)

                                                                            cgatcatGAATGCGATGGGTAACGATTACAGACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATg                     <  1:10293032/94‑1 (MQ=255)

                                                                            ccATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATg                     >  2:3823896/1‑96 (MQ=255)

                                                                             catcatGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGc                    >  2:3060927/1‑96 (MQ=255)

                                                                                catGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTa                 <  2:10461273/96‑1 (MQ=255)

                                                                                 atGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTat                >  1:9140795/1‑95 (MQ=255)

                                                                                   gAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACtt              <  2:6483709/96‑1 (MQ=255)

                                                                                     aTGCGATGGGTAACGATTACATACGTGGGCTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTta            <  1:1594418/96‑1 (MQ=255)

                                                                                       gCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtatt          <  1:9616420/96‑2 (MQ=255)

                                                                                       gCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtata          <  2:2839462/96‑1 (MQ=255)

                                                                                         gATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattca        >  2:13138902/1‑93 (MQ=255)

                                                                                         gATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattca        >  1:5751726/1‑93 (MQ=255)

                                                                                         gATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATaga        >  1:13555389/1‑96 (MQ=255)

                                                                                         gATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATaga        >  1:11829222/1‑96 (MQ=255)

                                                                                          aTGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATagag       >  2:8640322/1‑96 (MQ=255)

                                                                                           tGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGc      >  1:7085877/1‑96 (MQ=255)

                                                                                            gggCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCa     <  2:14023494/96‑1 (MQ=255)

                                                                                               cAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATaa  >  2:3209410/1‑96 (MQ=255)

GCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAA  >  NC_000913/580676‑580866

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