BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,531,073	T→C	M94V (ATG→GTG)	yjhQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,531,073	0	T	C	79.5%	574.6 / 164.9	302	M94V (ATG→GTG)	yjhQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (29/33); new base C (96/144); total (125/177)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

GGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAA  >  NC_000913/4530984‑4531164

ggTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCAcccc                                                                                          <  2:2068141/93‑1 (MQ=255)

 gtggcgtggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCt                                                                                         <  1:3353989/93‑1 (MQ=255)

  tggcgtggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCtt                                                                                        <  2:3118994/93‑1 (MQ=255)

      gtggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTa                                                                                    <  2:7298041/93‑1 (MQ=255)

       tggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTAt                                                                                   >  2:3718121/1‑93 (MQ=255)

        ggcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTAtt                                                                                  >  1:477578/1‑93 (MQ=255)

         gcCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTc                                                                                 <  1:44575/93‑1 (MQ=255)

           cGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCg                                                                               >  2:5124505/1‑93 (MQ=255)

            gagaACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCgg                                                                              >  1:4388264/1‑93 (MQ=255)

             agaACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGt                                                                             >  1:1213829/1‑93 (MQ=255)

             agaACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGt                                                                             <  1:1285696/93‑1 (MQ=255)

              gaACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTa                                                                            <  2:7522982/93‑1 (MQ=255)

              gaACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTa                                                                            >  1:1408244/1‑93 (MQ=255)

              gaACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTa                                                                            <  1:712485/93‑1 (MQ=255)

              gaACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTa                                                                            <  1:1941239/93‑1 (MQ=255)

              gaACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTa                                                                            <  2:6650874/93‑1 (MQ=255)

               aaCAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTAt                                                                           <  2:1696499/93‑1 (MQ=255)

                aCAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATg                                                                          >  1:5544914/1‑93 (MQ=255)

                aCAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATg                                                                          <  1:3060088/93‑1 (MQ=255)

                 cAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGa                                                                         <  2:6160885/93‑1 (MQ=255)

                 cAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGa                                                                         <  1:6965123/93‑1 (MQ=255)

                  aaaaaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGAc                                                                        <  2:2401965/93‑1 (MQ=255)

                  aaaaaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGAc                                                                        <  1:2795474/93‑1 (MQ=255)

                   aaaaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGAcc                                                                       <  1:1990851/93‑1 (MQ=255)

                   aaaaCCGTCTGGCAGCCCATCAACCGCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGAcc                                                                       <  1:4042338/93‑1 (MQ=255)

                    aaaCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCg                                                                      <  1:2655885/93‑1 (MQ=255)

                      aCCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGcc                                                                    <  1:5961963/93‑1 (MQ=255)

                       ccGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCa                                                                   <  2:5367405/93‑1 (MQ=255)

                       ccGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCa                                                                   <  1:6047700/93‑1 (MQ=255)

                        cGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAg                                                                  <  1:859572/93‑1 (MQ=255)

                        cGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAg                                                                  >  1:4350666/1‑93 (MQ=255)

                         gTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAgg                                                                 <  1:705801/93‑1 (MQ=255)

                          tCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAggg                                                                >  1:7146348/1‑93 (MQ=255)

                           cTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAgggg                                                               <  2:7218853/93‑1 (MQ=255)

                              gCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGc                                                            <  2:5183125/93‑1 (MQ=255)

                              gCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGc                                                            <  1:2919361/93‑1 (MQ=255)

                               cAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCg                                                           >  2:3538708/1‑93 (MQ=255)

                                aGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGa                                                          <  1:3746689/93‑1 (MQ=255)

                                aGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGa                                                          <  2:251358/93‑1 (MQ=255)

                                aGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGa                                                          <  1:6825395/93‑1 (MQ=255)

                                aGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGa                                                          <  2:6111690/93‑1 (MQ=255)

                                aGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGa                                                          >  2:1646339/1‑93 (MQ=255)

                                 gCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAg                                                         >  2:6271202/1‑93 (MQ=255)

                                 gCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAg                                                         <  1:1440381/93‑1 (MQ=255)

                                 gCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAg                                                         >  2:3347773/1‑93 (MQ=255)

                                  cccATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAgg                                                        <  2:1101322/93‑1 (MQ=255)

                                     aTCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGAtg                                                     >  1:3153091/1‑93 (MQ=255)

                                      tCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGAtgt                                                    >  2:2840738/1‑93 (MQ=255)

                                      tCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGAtgt                                                    <  2:6492675/93‑1 (MQ=255)

                                        aaCCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGc                                                  >  2:1773878/1‑93 (MQ=255)

                                             cAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCaa                                             >  1:3918238/1‑93 (MQ=255)

                                             cAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCaa                                             <  2:6030899/93‑1 (MQ=255)

                                               ggTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACg                                           >  2:3909495/1‑93 (MQ=255)

                                               ggTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACg                                           <  2:1221800/93‑1 (MQ=255)

                                                gTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACgg                                          <  2:3618742/93‑1 (MQ=255)

                                                 tGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGt                                         >  1:1607146/1‑93 (MQ=255)

                                                 tGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGt                                         <  1:5338537/93‑1 (MQ=255)

                                                 tGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGt                                         >  1:5878362/1‑93 (MQ=255)

                                                  gCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTACGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTg                                        <  2:4524689/93‑1 (MQ=255)

                                                     cGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAAt                                     <  1:178251/93‑1 (MQ=255)

                                                       aTGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATcc                                   >  2:6571586/1‑93 (MQ=255)

                                                        tGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCa                                  >  2:3445360/1‑93 (MQ=255)

                                                        tGCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCa                                  <  1:4961012/93‑1 (MQ=255)

                                                         gCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCAt                                 >  1:7458454/1‑93 (MQ=255)

                                                         gCCAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCAt                                 <  1:5629674/93‑1 (MQ=255)

                                                          ccAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCAtt                                >  1:3641225/1‑93 (MQ=255)

                                                          ccAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCAtt                                <  1:646809/93‑1 (MQ=255)

                                                          ccAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCAtt                                <  2:5967145/93‑1 (MQ=255)

                                                           cAGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCAttt                               >  1:6715989/1‑93 (MQ=255)

                                                            aGTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTtc                              <  2:1423816/93‑1 (MQ=255)

                                                             gTGCGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTtct                             >  1:4178532/1‑93 (MQ=255)

                                                                cGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCt                          >  2:1496335/1‑93 (MQ=255)

                                                                cGTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCt                          <  2:4042569/93‑1 (MQ=255)

                                                                 gTATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCtt                         <  2:87752/93‑1 (MQ=255)

                                                                  tATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCttt                        <  2:837481/93‑1 (MQ=255)

                                                                  tATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCttt                        <  1:2021544/93‑1 (MQ=255)

                                                                  tATCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCttt                        <  2:1999123/93‑1 (MQ=255)

                                                                   aTCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTg                       <  1:992632/93‑1 (MQ=255)

                                                                   aTCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTg                       >  2:4213104/1‑93 (MQ=255)

                                                                    tCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGa                      <  1:5094096/93‑1 (MQ=255)

                                                                    tCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGa                      >  2:6782553/1‑93 (MQ=255)

                                                                    tCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGa                      >  2:2756266/1‑93 (MQ=255)

                                                                    tCAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGa                      <  1:6647889/93‑1 (MQ=255)

                                                                     cAACCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGaa                     <  1:3440981/93‑1 (MQ=255)

                                                                      aaCCGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAg                    <  1:225370/93‑1 (MQ=255)

                                                                        ccGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGt                  <  2:5718815/93‑1 (MQ=255)

                                                                         cGTCCTCCCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGtt                 <  2:6881056/93‑1 (MQ=255)

                                                                           tcctccCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGc               >  1:4528760/1‑93 (MQ=255)

                                                                           tcctccCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGc               <  2:2179882/93‑1 (MQ=255)

                                                                              tccCACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCg            >  2:2639586/1‑93 (MQ=255)

                                                                               cccACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCgg           <  2:7004129/93‑1 (MQ=255)

                                                                                ccACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCggg          <  1:1150886/93‑1 (MQ=255)

                                                                                 cACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGt         >  2:5519717/1‑93 (MQ=255)

                                                                                 cACGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGt         >  1:4099759/1‑93 (MQ=255)

                                                                                  aCGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGta        <  1:631036/93‑1 (MQ=255)

                                                                                  aCGCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGta        >  1:570751/1‑93 (MQ=255)

                                                                                    gCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaa      <  2:102753/93‑1 (MQ=255)

                                                                                    gCCCACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaa      <  2:3028274/93‑1 (MQ=255)

                                                                                     cccACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaat     <  1:6221388/93‑1 (MQ=255)

                                                                                      ccACCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGtaaata    <  1:5058470/93‑1 (MQ=255)

GGTGGCGTGGCCGAGAACAAAAACCGTCTGGCAGCCCATCAACCTCAGGTGCTCGATGCCAGTGCGTATCAACCGTCCTCCCACGCCCATCCCTTGGTATTCCGGTATGACCGCCAGGGGAGCGAGGATGTGCATCAACGGTGAATCCATTTCTCCTTTGAAGGTTGCCCGGGTAAATAAA  >  NC_000913/4530984‑4531164

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