BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,555	(G)_5→6	coding (55/618 nt)	b3646 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,550	1	.	G	99.2%	1583.4 / 2.5	496	coding (50/618 nt)	b3646	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/4); new base G (166/325); total (167/329)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 573 total aligned reads displayed.

CACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGT  >  NC_000913/3816458‑3816637

cacaAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAc                                                                                            <  1:1128574/93‑1 (MQ=255)

cacaAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAc                                                                                            >  1:6612917/1‑93 (MQ=255)

 acaAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑‑g                                                                                         <  1:1648957/93‑1 (MQ=255)

   aaaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑‑ggg                                                                                       <  1:8104520/93‑1 (MQ=255)

   aaaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑‑ggg                                                                                       <  1:5591605/93‑1 (MQ=255)

    aaTCGCTGCCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑‑gggg                                                                                      <  1:946261/93‑1 (MQ=255)

    aaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑‑gggg                                                                                      <  2:672796/93‑1 (MQ=255)

    aaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑‑gggg                                                                                      >  1:2402607/1‑93 (MQ=255)

     aTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑‑ggggg                                                                                     <  2:9936011/93‑1 (MQ=255)

     aTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑‑ggggg                                                                                     >  1:2029142/1‑93 (MQ=255)

      tCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑‑gggggg                                                                                    <  2:2163260/93‑2 (MQ=255)

       cGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑‑ggggggc                                                                                   <  1:8430360/93‑3 (MQ=255)

           aCCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTg                                                                               <  1:7360799/93‑1 (MQ=255)

           aCCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTg                                                                               >  1:4277563/1‑93 (MQ=255)

             ccTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTggc                                                                             <  1:2197469/93‑1 (MQ=255)

             ccTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTggc                                                                             <  2:904074/93‑1 (MQ=255)

              cTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTggcg                                                                            >  1:5884496/1‑93 (MQ=255)

                gCCAGACAGATTTTTAGGGAGAGAACCATGCTTTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTggcggc                                                                          <  2:5650351/93‑1 (MQ=255)

                gCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTggcggc                                                                          <  2:6377600/93‑1 (MQ=255)

                gCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTggcggc                                                                          <  1:1776741/93‑1 (MQ=255)

                  cagacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCg                                                                        <  2:6568540/93‑1 (MQ=255)

                  cagacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCg                                                                        >  2:4226377/1‑93 (MQ=255)

                  cagacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCg                                                                        >  1:4422871/1‑93 (MQ=255)

                   agacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCgg                                                                       <  1:7959270/93‑1 (MQ=255)

                    gacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGga                                                                      <  2:7863787/93‑1 (MQ=255)

                    gacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGga                                                                      <  2:5376375/93‑1 (MQ=255)

                     acagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGgac                                                                     >  2:4431144/1‑93 (MQ=255)

                      cagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGgacg                                                                    <  2:3681368/93‑1 (MQ=255)

                       agaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGgacga                                                                   <  1:7842004/93‑1 (MQ=255)

                       agaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGgacga                                                                   <  1:288089/93‑1 (MQ=255)

                       agaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGgacga                                                                   <  2:1763268/93‑1 (MQ=255)

                         aTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCGGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGacgacg                                                                 <  2:4184821/93‑1 (MQ=255)

                         aTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGacgacg                                                                 <  1:7225479/93‑1 (MQ=255)

                          tttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACgg                                                                <  2:6157370/93‑1 (MQ=255)

                          tttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACgg                                                                <  2:1987739/93‑1 (MQ=255)

                          tttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACgg                                                                <  2:9294473/93‑1 (MQ=255)

                           ttttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGgc                                                               <  1:8258854/93‑1 (MQ=255)

                              tAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCa                                                            <  1:3124044/93‑1 (MQ=255)

                              tAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCa                                                            <  1:3900802/93‑1 (MQ=255)

                              tAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCa                                                            <  2:3276566/93‑1 (MQ=255)

                               aGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCAt                                                           <  2:3932158/93‑1 (MQ=255)

                               aGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCAt                                                           >  1:4792154/1‑93 (MQ=255)

                               aGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCAt                                                           <  1:255054/93‑1 (MQ=255)

                               aGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCAt                                                           <  2:4922929/93‑1 (MQ=255)

                                ggaagagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATg                                                          <  2:5855291/90‑1 (MQ=255)

                                 ggAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATgg                                                         <  1:4525572/93‑1 (MQ=255)

                                  gagagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGa                                                        <  2:3615185/93‑1 (MQ=255)

                                   agagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGAt                                                       <  2:7301235/93‑1 (MQ=255)

                                    gagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATa                                                      <  2:2298556/93‑1 (MQ=255)

                                    gagaACCATGCTGTTACACATTTTGAATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATa                                                      <  2:6288671/93‑1 (MQ=255)

                                      gaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACa                                                    >  1:6734805/1‑93 (MQ=255)

                                      gaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACa                                                    <  2:8837653/93‑1 (MQ=255)

                                      gaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACa                                                    <  1:536983/93‑1 (MQ=255)

                                        aCCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑AGCTGGCGGCCGGACGACGGCGCATGGATACAtt                                                  >  1:1472636/1‑93 (MQ=255)

                                         ccATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACAttt                                                 >  1:687531/1‑93 (MQ=255)

                                          cATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTg                                                >  1:2684669/1‑93 (MQ=255)

                                           aTGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTgg                                               >  2:1164732/1‑93 (MQ=255)

                                           aTGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTgg                                               >  2:3823965/1‑93 (MQ=255)

                                            tGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGc                                              >  1:3681177/1‑93 (MQ=255)

                                             gCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGc                                              <  1:8884058/93‑1 (MQ=255)

                                             gCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCAAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCg                                             >  2:6907002/1‑93 (MQ=255)

                                               tGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTa                                           >  2:8760493/1‑93 (MQ=255)

                                                  tACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAatt                                        >  2:7757593/1‑93 (MQ=255)

                                                     acaTTTTTTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatt                                     <  2:7611810/93‑1 (MQ=255)

                                                     acaTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatt                                     <  1:8651522/93‑1 (MQ=255)

                                                      caTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatta                                    <  2:1558613/93‑1 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   <  2:7140293/93‑1 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   >  2:9471129/1‑93 (MQ=255)

                                                        ttttGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattatt                                  <  1:6264699/93‑1 (MQ=255)

                                                        ttttGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattatt                                  >  2:4757047/1‑93 (MQ=255)

                                                         tttGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTg                                 <  1:503134/93‑1 (MQ=255)

                                                           tGTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCg                               <  2:2804677/93‑1 (MQ=255)

                                                            gTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGa                              <  1:9653657/93‑1 (MQ=255)

                                                            gTATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGa                              <  1:5111305/93‑1 (MQ=255)

                                                             tATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             <  1:1007438/93‑1 (MQ=255)

                                                             tATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             <  2:9869597/93‑1 (MQ=255)

                                                             tATTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             <  2:514129/93‑1 (MQ=255)

                                                             gaTTTGGTTGGCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             <  1:2096975/92‑1 (MQ=255)

                                                                 tggttggCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGcc                         <  1:5538849/93‑1 (MQ=255)

                                                                    ttggCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCAcc                      <  2:3243102/93‑1 (MQ=255)

                                                                     tggCATTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCg                     >  2:8317477/1‑93 (MQ=255)

                                                                         aTTCCTGCCGAAGCCATGAC‑‑GGGGGCCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGCCCGCCACCGCAAt                 <  2:1320118/93‑1 (MQ=255)

                                                                         aTTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  2:8483386/93‑1 (MQ=255)

                                                                         aTTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  1:1570304/93‑1 (MQ=255)

                                                                         aTTACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  2:8102692/93‑1 (MQ=255)

                                                                          ttACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAtt                >  1:6048402/1‑93 (MQ=255)

                                                                          ttACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAtt                <  2:8038824/93‑1 (MQ=255)

                                                                           tACTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTg               <  2:9618863/93‑1 (MQ=255)

                                                                            aCTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTgg              <  2:3442500/93‑1 (MQ=255)

                                                                            aCTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTgg              <  1:356058/93‑1 (MQ=255)

                                                                            aCTGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTgg              <  2:4514835/93‑1 (MQ=255)

                                                                              tGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCg            <  1:9077338/93‑1 (MQ=255)

                                                                              tGCCGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCg            >  1:5764113/1‑93 (MQ=255)

                                                                                ccGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGa          <  1:401642/93‑1 (MQ=255)

                                                                                 cGAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAg         <  1:1261976/93‑1 (MQ=255)

                                                                                  gAAGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAgg        >  2:2622008/1‑93 (MQ=255)

                                                                                    aGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGt      <  1:9985971/93‑1 (MQ=255)

                                                                                    aGCCATGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGt      >  2:5484191/1‑93 (MQ=255)

                                                                                     gCCATGACGGGGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGt      <  1:9337208/93‑1 (MQ=255)

                                                                                        aTGACG‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGt  <  1:3461302/93‑1 (MQ=255)

CACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑‑GGGGG‑CGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGT  >  NC_000913/3816458‑3816637

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