BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	171,718	C→T	P86L (CCG→CTG)	fhuB →	FhuB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	171,718	0	C	T	84.8%	783.8 / 91.7	289	P86L (CCG→CTG)	fhuB	FhuB
Reads supporting (aligned to +/- strand): ref base C (25/19); new base T (114/131); total (139/150)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 294 total aligned reads displayed.

TACAGCTTGTTGCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCCGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTGCGATGGC  >  NC_000913/171627‑171808

tACAGCTTGTTGCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCCg                                                                                           >  1:1078862/1‑93 (MQ=255)

 aCAGCTTGTTGCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGc                                                                                          >  2:3810691/1‑93 (MQ=255)

 aCAGCTTGTTGCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGc                                                                                          <  1:5925896/93‑1 (MQ=255)

  cAGCTTGTTGCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCt                                                                                         <  2:4562000/93‑1 (MQ=255)

   aGCTTGTTGCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACctgctg                                                                                        >  1:6850795/1‑93 (MQ=255)

   aGCTTGTTGCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACctgctg                                                                                        <  1:7018071/93‑1 (MQ=255)

   aGCTTGTTGCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCCGCTg                                                                                        <  2:6246467/93‑1 (MQ=255)

    gCTTGTTGCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTgg                                                                                       <  1:1641075/93‑1 (MQ=255)

     cTTGTTGCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGc                                                                                      <  2:3948366/93‑1 (MQ=255)

     cTTGTTGCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGc                                                                                      <  1:3035671/93‑1 (MQ=255)

       tgttgCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCgg                                                                                    <  2:554664/93‑1 (MQ=255)

        gttgCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGa                                                                                   <  2:5994835/93‑1 (MQ=255)

            ccGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCCGCTGGCGGAGCcg                                                                               <  2:7519025/93‑1 (MQ=255)

             cgcgTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCcga                                                                              <  2:3152982/93‑1 (MQ=255)

              gcgTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCcgac                                                                             <  2:3403364/93‑1 (MQ=255)

              gcgTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCcgac                                                                             <  1:1858380/93‑1 (MQ=255)

              gcgTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCcgac                                                                             >  2:6599055/1‑93 (MQ=255)

              gcgTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCcgac                                                                             >  1:6303332/1‑93 (MQ=255)

                gTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCcgacga                                                                           >  2:918187/1‑93 (MQ=255)

                   tGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCCGCTGGCGGAGCCGACGACGc                                                                        >  2:1195274/1‑93 (MQ=255)

                    ggCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCt                                                                       <  2:1512972/93‑1 (MQ=255)

                       gATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTgg                                                                    >  1:2775372/1‑93 (MQ=255)

                       gATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAACTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTgg                                                                    <  2:2605/93‑1 (MQ=255)

                        aTTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGc                                                                   >  2:6268230/1‑93 (MQ=255)

                        aTTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGc                                                                   >  1:3052920/1‑93 (MQ=255)

                        aTTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGc                                                                   <  2:3575680/93‑1 (MQ=255)

                        aTTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGc                                                                   >  1:438214/1‑93 (MQ=255)

                        aTTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTCCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGc                                                                   >  1:539157/1‑93 (MQ=255)

                          ttCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGt                                                                 <  2:4765660/93‑1 (MQ=255)

                          ttCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGt                                                                 <  2:4902024/93‑1 (MQ=255)

                              ctgctgGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCt                                                             >  2:437678/1‑93 (MQ=255)

                                gctgGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACCCTTGGCGTTGCTAc                                                           <  2:5189708/93‑1 (MQ=255)

                                 ctgGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACa                                                          >  1:5617561/1‑93 (MQ=255)

                                 ctgGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACa                                                          <  1:3400448/93‑1 (MQ=255)

                                  tggtggGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAg                                                         >  2:3322921/1‑93 (MQ=255)

                                  tggtggGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAg                                                         >  1:7055518/1‑93 (MQ=255)

                                    gtggGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGgc                                                       <  2:5846888/93‑1 (MQ=255)

                                     tggGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGgcg                                                      <  1:6457137/93‑1 (MQ=255)

                                       ggCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGgcgcg                                                    >  1:2070430/1‑93 (MQ=255)

                                           cgGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAAc                                                <  1:958967/93‑1 (MQ=255)

                                             ggTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTg                                              <  1:162152/93‑1 (MQ=255)

                                              gTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTgg                                             >  2:4290838/1‑93 (MQ=255)

                                              gTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTgg                                             <  1:35104/93‑1 (MQ=255)

                                               tCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTggg                                            >  1:5783279/1‑93 (MQ=255)

                                               tCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTggg                                            <  1:7454482/93‑1 (MQ=255)

                                                cTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTgggg                                           <  1:1229443/93‑1 (MQ=255)

                                                 tGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGa                                          <  2:129745/93‑1 (MQ=255)

                                                  ggggCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCCGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGAt                                         <  2:5429551/93‑1 (MQ=255)

                                                  ggggCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCCGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGAt                                         <  2:5817934/93‑1 (MQ=255)

                                                   gggCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGAtt                                        >  1:4586192/1‑93 (MQ=255)

                                                    ggCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTa                                       <  1:1932762/93‑1 (MQ=255)

                                                    ggCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTa                                       <  2:3719486/93‑1 (MQ=255)

                                                     gCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTAc                                      >  2:7169189/1‑93 (MQ=255)

                                                       tggtggGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCg                                    >  2:6897911/1‑93 (MQ=255)

                                                       tggtggGCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCg                                    >  1:6619349/1‑93 (MQ=255)

                                                           gggCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCAc                                <  1:5435513/93‑1 (MQ=255)

                                                            ggCGTGCTGTTTCCGCCAGGTGTGCGTAACCTGCTGGCGGGGCCGACGACGCTTTGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACt                               <  2:6034131/93‑1 (MQ=255)

                                                            ggCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACt                               >  1:6524815/1‑93 (MQ=255)

                                                             gCGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTa                              <  2:1450548/93‑1 (MQ=255)

                                                              cGTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTAc                             <  1:856936/93‑1 (MQ=255)

                                                              cGTGCTGTTTCAGCAAGTGCTGCGTAACCCGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTAc                             <  2:2824596/93‑1 (MQ=255)

                                                               gTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACg                            >  2:1391015/1‑93 (MQ=255)

                                                               gTGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACg                            <  1:2620674/93‑1 (MQ=255)

                                                                tGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGc                           >  2:4247062/1‑93 (MQ=255)

                                                                tGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGc                           >  1:5225996/1‑93 (MQ=255)

                                                                tGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGc                           >  2:2225867/1‑93 (MQ=255)

                                                                tGCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGc                           >  1:5579098/1‑93 (MQ=255)

                                                                 gCTGTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGct                          >  1:2837638/1‑93 (MQ=255)

                                                                    gTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTg                       <  2:5311420/93‑1 (MQ=255)

                                                                    gTTTCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTg                       <  1:1409263/93‑1 (MQ=255)

                                                                    gTTTCAGCAAGTGCTGCGTAACCCGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTg                       >  1:6701269/1‑93 (MQ=255)

                                                                     tttCAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTgg                      >  2:2677201/1‑93 (MQ=255)

                                                                        cAGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCg                   >  2:1628622/1‑93 (MQ=255)

                                                                         aGCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGa                  >  1:385086/1‑93 (MQ=255)

                                                                          gCAAGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGAt                 >  1:4037822/1‑93 (MQ=255)

                                                                            aaGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATcc               <  1:3305648/93‑1 (MQ=255)

                                                                             aGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATccc              <  2:1805811/93‑1 (MQ=255)

                                                                             aGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATccc              <  2:2186048/93‑1 (MQ=255)

                                                                             aGTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATccc              >  2:6195495/1‑93 (MQ=255)

                                                                              gTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCt             >  1:5669879/1‑93 (MQ=255)

                                                                              gTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCt             <  1:7184709/93‑1 (MQ=255)

                                                                              gTGCTGCGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCt             >  1:2275626/1‑93 (MQ=255)

                                                                               tgctgcGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTg            <  1:2934111/93‑1 (MQ=255)

                                                                               tgctgcGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTg            <  1:4154633/93‑1 (MQ=255)

                                                                                 ctgcGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAGCTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGt          <  2:2037350/93‑1 (MQ=255)

                                                                                 ctgcGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGt          <  1:1503106/93‑1 (MQ=255)

                                                                                 ctgcGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGt          >  2:2475558/1‑93 (MQ=255)

                                                                                 ctgcGTAACCCGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGt          <  1:4244283/93‑1 (MQ=255)

                                                                                  tgcGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTg         >  2:1853172/1‑93 (MQ=255)

                                                                                  tgcGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTg         <  1:390413/93‑1 (MQ=255)

                                                                                   gcGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTGc        >  1:991899/1‑93 (MQ=255)

                                                                                   gcGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTGc        >  1:3223590/1‑93 (MQ=255)

                                                                                    cGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTGCg       >  2:3075970/1‑93 (MQ=255)

                                                                                    cGTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTGCg       >  1:5643640/1‑93 (MQ=255)

                                                                                     gTAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTGCGa      >  2:6747729/1‑93 (MQ=255)

                                                                                      tAACCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTGCGAt     >  2:4729833/1‑93 (MQ=255)

                                                                                        aCCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTGCGATgg   <  2:2741892/93‑1 (MQ=255)

                                                                                        aCCTGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTGCGATgg   <  1:4912415/93‑1 (MQ=255)

                                                                                        aCCCGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTGCGATgg   <  1:2195783/93‑1 (MQ=255)

                                                                                         cccGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTGCGATGGc  >  1:3720002/1‑93 (MQ=255)

TACAGCTTGTTGCCGCGTCTGGCGATTTCGCTGCTGGTGGGCGCGGGTCTGGGGCTGGTGGGCGTGCTGTTTCAGCAAGTGCTGCGTAACCCGCTGGCGGAGCCGACGACGCTTGGCGTTGCTACAGGCGCGCAACTGGGGATTACCGTCACTACGCTCTGGGCGATCCCTGGTGCGATGGC  >  NC_000913/171627‑171808

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