BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,949,410	0	C	T	71.3%	1624.4 / inf	949	L196L (CTC→CTT)	yecD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (157/115); new base T (363/313); total (520/429)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 980 total aligned reads displayed.

GCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCC  >  NC_000913/1949318‑1949500

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

       gTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGt                                                                                     >  2:7022834/1‑93 (MQ=255)

        tGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGtt                                                                                    <  1:8471123/93‑1 (MQ=255)

        tGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAACCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGtt                                                                                    >  1:9682952/1‑93 (MQ=255)

         gccgccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTa                                                                                   <  1:3188174/93‑1 (MQ=255)

         gccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTa                                                                                   <  1:5881452/93‑1 (MQ=255)

         gccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTa                                                                                   >  2:1465451/1‑93 (MQ=255)

         gccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTa                                                                                   <  1:420213/93‑1 (MQ=255)

          ccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTAt                                                                                  <  2:978647/93‑1 (MQ=255)

          ccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTAt                                                                                  <  2:1761415/93‑1 (MQ=255)

          ccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTAt                                                                                  >  1:3554467/1‑93 (MQ=255)

            gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATgt                                                                                >  1:7841308/1‑92 (MQ=255)

            gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGa                                                                                >  1:6100935/1‑93 (MQ=255)

            gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGa                                                                                <  1:1765321/93‑1 (MQ=255)

            cccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGa                                                                                <  1:257221/92‑1 (MQ=255)

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

              cAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAtt                                                                              <  2:7911466/93‑1 (MQ=255)

              cAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAtt                                                                              <  2:842906/93‑1 (MQ=255)

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

                gcgcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTa                                                                            <  1:4805100/93‑1 (MQ=255)

                gcgcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTa                                                                            <  2:7340414/93‑1 (MQ=255)

                gcgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTa                                                                            <  1:8186817/93‑1 (MQ=255)

                gcgcCAAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTa                                                                            >  2:7478230/1‑93 (MQ=255)

                  gcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACa                                                                          >  2:3206646/1‑93 (MQ=255)

                  gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACa                                                                          <  2:9466422/93‑1 (MQ=255)

                     gAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCg                                                                       <  2:9597119/93‑1 (MQ=255)

                     gAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCg                                                                       <  1:1021150/93‑1 (MQ=255)

                         aacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCt                                                                   >  2:4401247/3‑93 (MQ=255)

                         aacacaATAACAGCATTAATCATATCCCCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCt                                                                   <  1:6598932/91‑1 (MQ=255)

                           cacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtgc                                                                 <  2:5750724/93‑3 (MQ=255)

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

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

                               aTAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCa                                                             <  1:3923287/93‑1 (MQ=255)

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

                                   cAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATgg                                                         <  1:905574/93‑1 (MQ=255)

                                     gTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTc                                                       <  2:394317/93‑1 (MQ=255)

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

                                       aTCACCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGc                                                     <  2:528171/93‑1 (MQ=255)

                                          aatcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATc                                                  >  2:5156911/4‑93 (MQ=255)

                                           atcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATcc                                                 <  2:4889206/91‑1 (MQ=255)

                                             cATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTa                                               <  1:6941683/93‑1 (MQ=255)

                                               tatCTTTCCGCTCATCGCCCGTGTGCTTAGCGTCTAAGAGTTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                             <  2:2866663/93‑1 (MQ=255)

                                               tatCTTCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTaaa                                             <  2:3387818/93‑1 (MQ=255)

                                               tatCTACCCGCGCATCGCCGGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTaaa                                             <  1:99438/93‑1 (MQ=255)

                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                             >  1:440922/1‑93 (MQ=255)

                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                             >  2:8556305/1‑93 (MQ=255)

                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                             <  1:8608340/93‑1 (MQ=255)

                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                             >  2:5623193/1‑93 (MQ=255)

                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                             >  2:5927650/1‑93 (MQ=255)

                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTaaa                                             >  2:6557576/1‑93 (MQ=255)

                                                 tcgaCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATg                                           <  1:3754943/90‑1 (MQ=255)

                                                     cccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTg                                       >  1:4560657/1‑93 (MQ=255)

                                                      ccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGc                                      <  1:9341636/93‑1 (MQ=255)

                                                      ccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGc                                      <  2:5466054/93‑1 (MQ=255)

                                                       cgcgCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCg                                     <  2:2465775/93‑1 (MQ=255)

                                                       cgcgCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCg                                     >  1:5643620/1‑93 (MQ=255)

                                                        gcgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCgg                                    >  2:4627264/1‑93 (MQ=255)

                                                        gcgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCgg                                    >  2:1183779/1‑93 (MQ=255)

                                                         cgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGt                                   <  1:7286136/93‑1 (MQ=255)

                                                          gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                  >  2:2395805/1‑93 (MQ=255)

                                                          gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                  <  1:9155574/93‑1 (MQ=255)

                                                            aTCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTgg                                <  2:6697119/93‑1 (MQ=255)

                                                            aTCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTgg                                <  1:6312211/93‑1 (MQ=255)

                                                            aTCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTgg                                >  1:1390780/1‑93 (MQ=255)

                                                             tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                               >  1:9517266/1‑93 (MQ=255)

                                                             tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                               >  2:8157221/1‑93 (MQ=255)

                                                             tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                               >  2:2141563/1‑93 (MQ=255)

                                                              cgcgcgTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTgggt                              <  1:294901/93‑2 (MQ=255)

                                                               gCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTGGCATCCTAAGTGGGTGCGGGTgggta                             >  1:8356537/1‑91 (MQ=255)

                                                                cccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTAt                            >  2:2688821/1‑93 (MQ=255)

                                                                cccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTAt                            <  1:4930027/93‑1 (MQ=255)

                                                                  cGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATcc                          >  2:8771839/1‑92 (MQ=255)

                                                                   gtgtGCGTAGCGCGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAc                         <  1:7503417/93‑1 (MQ=255)

                                                                    tgtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAcc                        >  1:7133565/1‑93 (MQ=255)

                                                                    tgtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAcc                        >  1:1523547/1‑93 (MQ=255)

                                                                     gtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCa                       <  1:3341916/93‑1 (MQ=255)

                                                                      tgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAg                      >  2:804523/1‑93 (MQ=255)

                                                                        cGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGcc                    <  2:364117/93‑1 (MQ=255)

                                                                        cGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGcc                    <  2:1994623/93‑1 (MQ=255)

                                                                          tAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCtt                  <  2:5048050/93‑1 (MQ=255)

                                                                            gCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGa                <  1:5167178/93‑1 (MQ=255)

                                                                             cGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGaa               <  1:5444997/93‑1 (MQ=255)

                                                                              gTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAag              <  2:641638/93‑1 (MQ=255)

                                                                              gTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAag              <  2:4021506/93‑1 (MQ=255)

                                                                              gTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAag              <  2:1017014/93‑1 (MQ=255)

                                                                              gTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAag              <  1:4678879/93‑1 (MQ=255)

                                                                               tGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAaga             <  1:6425271/93‑1 (MQ=255)

                                                                                ggAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAagag            <  1:8805383/93‑1 (MQ=255)

                                                                                  aaGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTa          <  1:203348/93‑1 (MQ=255)

                                                                                  aaGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTa          >  2:2996515/1‑93 (MQ=255)

                                                                                  aaGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTa          >  2:7616299/1‑93 (MQ=255)

                                                                                   agagATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTAt         <  1:1000143/93‑1 (MQ=255)

                                                                                     agaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGc       >  2:4240906/1‑93 (MQ=255)

                                                                                     agaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGc       >  1:7719937/1‑93 (MQ=255)

                                                                                     agaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGc       <  1:6754696/93‑1 (MQ=255)

                                                                                     agaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGCATGc       >  1:1878158/1‑93 (MQ=255)

                                                                                      gaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGcc      <  1:5281725/93‑1 (MQ=255)

                                                                                       aTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGccc     <  2:6822443/93‑1 (MQ=255)

                                                                                       aTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGccc     <  2:7244029/93‑1 (MQ=255)

                                                                                       aTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGccc     >  2:9286275/1‑93 (MQ=255)

GCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCC  >  NC_000913/1949318‑1949500

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