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.4%	100.4 / inf	1074	I168I (ATC→ATT)	yecD	isochorismatase family protein
Reads supporting (aligned to +/- strand): ref base C (279/285); new base T (246/262); total (526/548)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-06
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 1123 total aligned reads displayed.

TTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATG  >  NC_000913/1949259‑1949459

    aTCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAAcc                                                                                                  >  1:14370178/1‑101 (MQ=255)

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

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

           nnncgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGTCAGCGCCGAGCAGCACAATACCAGTATCAACCATCTCTa                                                                                           <  2:13006328/98‑1 (MQ=255)

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

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

             ccgcAATGCCTGGGAACTTGGCTTTAACCTGGTGATTCCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAcc                                                                                         <  1:632000/101‑1 (MQ=255)

              nnnaaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAccc                                                                                        <  2:12061752/98‑1 (MQ=255)

              cgcAATGCCTGGTAACTTGGCTTTAATCTGGTGATTGCCGACGCCGCCTGTAGTGCCGCCAGCGCCCAGCAGCACAATAACAGTATCAACCATATCTAccc                                                                                        <  1:511117/101‑1 (MQ=255)

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

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

                  aTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCgcgc                                                                                    >  1:8168780/1‑101 (MQ=255)

                  aTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCgcgc                                                                                    >  1:3013309/1‑101 (MQ=255)

                  aTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCgcgc                                                                                    <  2:1437242/101‑1 (MQ=255)

                     nnngggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATc                                                                                 <  2:7954401/98‑1 (MQ=255)

                     ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATc                                                                                 <  2:3566201/101‑1 (MQ=255)

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

                      cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCg                                                                                <  2:5564868/101‑1 (MQ=255)

                       taggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCAtnnn                                                                               >  2:11037527/3‑98 (MQ=255)

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

                       tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGc                                                                               >  1:12711720/1‑101 (MQ=255)

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

                           nnnttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCgt                                                                           <  2:7705593/98‑1 (MQ=255)

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

                               tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGc                                                                       >  1:371233/1‑101 (MQ=255)

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

                                  nnntAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTa                                                                    <  2:8512839/98‑1 (MQ=255)

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

                                    ttAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCGGCGCATCGCCCGTGGGCGTGGc                                                                  >  1:1638259/1‑101 (MQ=255)

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

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

                                    ttAATCTGGTGATTGCCGAACACGCCTGTAGTGCCGCCAGCGCAGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGc                                                                  <  1:8993785/101‑1 (MQ=255)

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

                                      aaTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAgnnn                                                                >  2:12671805/1‑98 (MQ=255)

                                       aTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTg                                                               >  1:8758225/1‑101 (MQ=255)

                                         cTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGa                                                             >  2:1593893/1‑101 (MQ=255)

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

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

                                            gTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCGCAATAACAGTATCAACCATATCTCCCCGCGCATCGCCCGTGTGCGTAGCGTGGAaga                                                          <  2:3322156/101‑1 (MQ=255)

                                            gTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAaga                                                          <  2:6325193/101‑1 (MQ=255)

                                            gTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAaga                                                          <  2:3817273/101‑1 (MQ=255)

                                                 tGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAgagannn                                                     >  2:8394048/1‑98 (MQ=255)

                                                 tGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                                                     <  1:5380866/101‑1 (MQ=255)

                                                 nnncGATGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTATGCGTAGCGTGGAAGAGATcc                                                     <  2:10705071/98‑1 (MQ=255)

                                                  gCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATccc                                                    >  2:4727650/1‑100 (MQ=255)

                                                  gCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGGAGCGTGGAAGAGATCCt                                                    >  2:1030533/1‑101 (MQ=255)

                                                    nnnaGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCa                                                  <  2:11216578/98‑1 (MQ=255)

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

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

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

                                                      aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTcnnn                                                >  2:9708300/1‑98 (MQ=255)

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

                                                      aaGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttcac                                                >  1:4650651/1‑97 (MQ=255)

                                                        gACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGGGCGTAGCGTGGAAGAGATCCTTCAcgc                                              >  1:9575937/1‑101 (MQ=255)

                                                         aCGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcgcg                                             <  1:4290166/101‑1 (MQ=255)

                                                         aCGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcgcg                                             <  1:9819517/101‑1 (MQ=255)

                                                          nnntcgTAGTGCCGCCAGCGGCGCGCAGCACAACACTAGTATCAACCATATCTACCCGCGAATGGCCCGTGTGCGTAGCGTGGAAGACGTCCTCAACGCGt                                            <  2:14641556/96‑1 (MQ=255)

                                                          cGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGt                                            >  1:1973695/1‑101 (MQ=255)

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

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

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

                                                            nnngTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTCCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTa                                          <  2:13906897/98‑1 (MQ=255)

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

                                                             nnntAGTGCCGCCAGCGCCGATCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTAt                                         <  2:13180838/98‑1 (MQ=255)

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

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

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

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

                                                              nnnaGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATg                                        <  2:7398625/98‑1 (MQ=255)

                                                                 aGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCGACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCTTGGAAGAGATCCTCAACGCGTTATgnnn                                     >  2:13473811/1‑98 (MQ=255)

                                                                  gTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAttt                                    <  2:822644/101‑1 (MQ=255)

                                                                   tGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTa                                   >  1:6134674/1‑101 (MQ=255)

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

                                                                    gccgccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTAc                                  <  1:12154695/101‑1 (MQ=255)

                                                                    gccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTAc                                  <  2:1817162/101‑1 (MQ=255)

                                                                    gccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGGGCGTAGCGTGGAGGAGATCCTTCACGCGTTATGATTTAc                                  >  1:413754/1‑101 (MQ=255)

                                                                      cgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTannn                                >  2:9016918/1‑98 (MQ=255)

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

                                                                        ccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCg                              <  2:2763551/101‑1 (MQ=255)

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

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

                                                                             nnngAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCtc                         <  2:11644144/98‑2 (MQ=255)

                                                                              nnnagcaacACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtgc                        <  2:9361613/98‑3 (MQ=255)

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

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

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

                                                                                gAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCg                      <  1:665567/101‑1 (MQ=255)

                                                                                  gcaacaCAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCa                    >  1:6950812/1‑101 (MQ=255)

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

                                                                                   caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCaa                   >  1:12391394/4‑101 (MQ=255)

                                                                                    nnnacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAAt                  <  2:11880521/98‑1 (MQ=255)

                                                                                     gcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCannn                 >  2:8779177/1‑98 (MQ=255)

                                                                                     gcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATg                 >  1:9311529/1‑101 (MQ=255)

                                                                                      cacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATgg                <  1:2322860/101‑1 (MQ=255)

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

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

                                                                                          aTAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGc            <  1:8394048/101‑1 (MQ=255)

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

                                                                                            aaCAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCac          >  1:5959814/1‑100 (MQ=255)

                                                                                                  aTCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTaaa    >  2:6124599/1‑101 (MQ=255)

TTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATG  >  NC_000913/1949259‑1949459

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