BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	247,207	A→G	T166A (ACG→GCG)	yafL →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	247,207	0	A	G	91.3%	758.9 / 47.8	289	T166A (ACG→GCG)	yafL	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base A (12/13); new base G (136/128); total (148/141)
Fisher's exact test for biased strand distribution p-value = 8.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 310 total aligned reads displayed.

TCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCACGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCAC  >  NC_000913/247115‑247296

gcTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCa                                                                                           <  1:2835143/92‑1 (MQ=255)

 cTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCAc                                                                                          <  2:3431792/93‑1 (MQ=255)

   ggggCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCgg                                                                                        >  1:8886823/1‑93 (MQ=255)

    gggCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGc                                                                                       <  1:8606807/93‑1 (MQ=255)

       cGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCaa                                                                                    >  1:3643306/1‑93 (MQ=255)

       cGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCACGGCCaa                                                                                    <  2:1313231/93‑1 (MQ=255)

         gTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAATATCCTTGAGGCTATGCTCCCGCGCGCGGCCAATg                                                                                  >  2:3790014/1‑93 (MQ=255)

          tACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATga                                                                                 >  2:8629246/1‑93 (MQ=255)

          tACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATga                                                                                 <  2:4128893/93‑1 (MQ=255)

          tACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATga                                                                                 <  2:4964799/93‑1 (MQ=255)

            cgcgTCCTGATAAAGGCTTTGACTGTAGCGGGTTGGGTTTTTATTCCTACAACAAGATCCTTGATTCTAAGCACCCTCGCGCGGCCGATgagg                                                                               >  2:7556499/1‑92 (MQ=255)

            cgcgGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATgaga                                                                               <  1:2129666/93‑1 (MQ=255)

            cgcgGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATgaga                                                                               <  2:7604155/93‑1 (MQ=255)

            cgcgGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATgaga                                                                               <  2:5999229/93‑1 (MQ=255)

             gcgGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGGGGCTAAGCTCCCGCGCGCGGCCAATGAGAt                                                                              >  1:4143365/1‑93 (MQ=255)

             gcgGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGAt                                                                              >  1:7797363/1‑93 (MQ=255)

             gcgGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCACGGCCAATGAGAt                                                                              <  1:8590039/93‑1 (MQ=255)

              cgGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATg                                                                             <  1:10014373/93‑1 (MQ=255)

              cgGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATg                                                                             <  2:4245658/93‑1 (MQ=255)

                 ccTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTAc                                                                          <  1:6913240/93‑1 (MQ=255)

                   tGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCa                                                                        >  1:7596933/1‑93 (MQ=255)

                   tGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCa                                                                        >  1:4549550/1‑93 (MQ=255)

                   tGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCa                                                                        >  2:7967245/1‑93 (MQ=255)

                    gATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCAc                                                                       >  2:7714602/1‑93 (MQ=255)

                    gATAAAGGCTTTGACTGTAGCGGGGTGGGTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCACGGCCCATGAGATGTAcccc                                                                       >  2:8091140/1‑91 (MQ=255)

                     aTAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACt                                                                      >  1:367843/1‑93 (MQ=255)

                      tAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTa                                                                     >  1:9092395/1‑93 (MQ=255)

                       aaaGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGCGATGTACCACTAt                                                                    >  1:5983954/1‑93 (MQ=255)

                        aaGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATc                                                                   <  1:7361670/93‑1 (MQ=255)

                        aaGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATc                                                                   <  1:795127/93‑1 (MQ=255)

                         aGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCg                                                                  >  1:1886300/1‑93 (MQ=255)

                          ggCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGc                                                                 >  2:6387576/1‑93 (MQ=255)

                           gCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGcc                                                                >  1:578359/1‑93 (MQ=255)

                                gACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCa                                                           <  1:9640264/93‑1 (MQ=255)

                                  cTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAAc                                                         <  1:5575428/93‑1 (MQ=255)

                                    gTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGa                                                       >  1:5085273/1‑93 (MQ=255)

                                     tAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGAt                                                      <  2:8919519/93‑1 (MQ=255)

                                      aGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGAtt                                                     <  2:587242/93‑1 (MQ=255)

                                       gCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATtg                                                    >  1:9396448/1‑93 (MQ=255)

                                       gCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATtg                                                    <  1:7713025/93‑1 (MQ=255)

                                         gggTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATtgtg                                                  >  1:5384024/1‑93 (MQ=255)

                                         gggTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTACCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATtgtg                                                  >  2:2914834/1‑93 (MQ=255)

                                          ggttggttTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTgg                                                 >  1:6184013/1‑93 (MQ=255)

                                          ggttggttTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTgg                                                 >  1:6826522/1‑93 (MQ=255)

                                          ggttggttTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTgg                                                 <  2:2235373/93‑1 (MQ=255)

                                          ggttggttTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCACGGCCAATGAGATGTACCACTATCGCCGGGCACCGATTGTgg                                                 >  2:1505853/1‑93 (MQ=255)

                                            ttggttTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCg                                               <  2:364013/93‑1 (MQ=255)

                                               gttTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGaac                                            <  1:1017681/93‑1 (MQ=255)

                                                 tttttATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGaacaa                                          >  1:13037/1‑93 (MQ=255)

                                                 ttgttATGCCTACACCAAGAACCTTGAGGCTAAGCCCCCGCGCGCGGCCAATGAGATGTACCCCTATCGCCGGGCAACGATTGTGGCGaacaa                                          <  2:6233181/90‑1 (MQ=255)

                                                   tttATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACg                                        <  1:5270805/93‑1 (MQ=255)

                                                   tttATGCCTACAACAAGATCCTTAAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGAGAACAACg                                        >  1:342631/1‑93 (MQ=255)

                                                    ttATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGa                                       <  1:6952312/93‑1 (MQ=255)

                                                      aTGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGAcc                                     <  1:4948153/93‑1 (MQ=255)

                                                      aTGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGAcc                                     <  1:3306561/93‑1 (MQ=255)

                                                        gCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTg                                   <  1:2376465/93‑1 (MQ=255)

                                                        gCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTg                                   <  1:2646402/93‑1 (MQ=255)

                                                         ccTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTgc                                  >  1:2955830/1‑93 (MQ=255)

                                                          cTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTgcg                                 >  1:7123049/1‑93 (MQ=255)

                                                          cTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTgcg                                 >  1:6849428/1‑93 (MQ=255)

                                                           tACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTgcgc                                >  1:8108247/1‑93 (MQ=255)

                                                           tACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTgcgc                                >  2:488463/1‑93 (MQ=255)

                                                           tACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTgcgc                                >  2:2057734/1‑93 (MQ=255)

                                                           tACAACAAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTgcgc                                >  2:6509975/1‑93 (MQ=255)

                                                           tACAACAAGATCCTTGAGGCTAAGCTCCCGCGCACGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTgcgc                                >  1:9355398/1‑93 (MQ=255)

                                                            acaacaAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGcc                               >  1:3185172/1‑93 (MQ=255)

                                                            acaacaAGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGcc                               <  2:2607819/93‑1 (MQ=255)

                                                             caacaaGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCg                              >  1:9191955/1‑93 (MQ=255)

                                                              aacaaGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCgg                             >  2:2446490/1‑93 (MQ=255)

                                                              aacaaGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCgg                             <  1:1365881/93‑1 (MQ=255)

                                                               acaaGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCggg                            <  1:3554449/93‑1 (MQ=255)

                                                               acaaGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCggg                            <  1:8438885/93‑1 (MQ=255)

                                                                caaGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCgggg                           <  1:1777769/93‑1 (MQ=255)

                                                                caaGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCgggg                           <  1:9775164/93‑1 (MQ=255)

                                                                 aaGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGga                          >  1:4329446/1‑93 (MQ=255)

                                                                  aGATCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGCGGCGAACAACGACCTGCGCCGGGgag                         >  1:6675127/1‑93 (MQ=255)

                                                                    aTCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGAt                       >  2:10033039/1‑93 (MQ=255)

                                                                     tCCTTGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACTATTGTGGCGAACAACGACCTGCGCCGGGGAGAtt                      >  2:2207970/1‑93 (MQ=255)

                                                                     accTTGAGGCTAAGCTCCCGCGCACGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGAtt                      <  1:1519824/92‑1 (MQ=255)

                                                                        ttGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGc                   <  2:1849700/93‑1 (MQ=255)

                                                                        ttGAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGc                   >  2:8322540/1‑93 (MQ=255)

                                                                          gAGGCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTg                 <  2:7155547/93‑1 (MQ=255)

                                                                           aGGCTAAGCTCCCGCGCACGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGt                <  2:3161483/93‑1 (MQ=255)

                                                                            ggCTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGtt               <  1:4831691/93‑1 (MQ=255)

                                                                              cTAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGtttt             <  2:4744176/93‑1 (MQ=255)

                                                                               tAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGttttt            >  1:309773/1‑93 (MQ=255)

                                                                               tAAGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGttttt            <  2:884071/93‑1 (MQ=255)

                                                                                aaGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTc           <  2:3081869/93‑1 (MQ=255)

                                                                                 aGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTcc          >  1:5236763/1‑93 (MQ=255)

                                                                                 aGCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTcc          <  2:6729180/93‑1 (MQ=255)

                                                                                  gCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCa         >  2:1623457/1‑93 (MQ=255)

                                                                                  gCTCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCa         <  2:927981/93‑1 (MQ=255)

                                                                                    tCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCata       <  2:9510682/93‑1 (MQ=255)

                                                                                    tCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCata       >  2:546431/1‑93 (MQ=255)

                                                                                    tCCCGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCata       <  2:4411559/93‑1 (MQ=255)

                                                                                    tCCCACGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCata       <  1:6541787/93‑1 (MQ=255)

                                                                                      ccGCGCGCGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATc     >  1:987189/1‑93 (MQ=255)

                                                                                       cgcgcgcgGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGACCAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATcc    <  2:7950004/93‑1 (MQ=255)

                                                                                         cgcgcgGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTTCATATCcac  >  1:383375/1‑93 (MQ=255)

                                                                                         cgcgcgGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCcac  <  2:9129311/93‑1 (MQ=255)

TCTGGGGCGGTACGCGGCCTGATAAAGGCTTTGACTGTAGCGGGTTGGTTTTTTATGCCTACAACAAGATCCTTGAGGCTAAGCTCCCGCGCACGGCCAATGAGATGTACCACTATCGCCGGGCAACGATTGTGGCGAACAACGACCTGCGCCGGGGAGATTTGCTGTTTTTCCATATCCAC  >  NC_000913/247115‑247296

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