BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,066,222	T→C	Y212C (TAC→TGC)	yihQ ←	putative glycosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,066,222	0	T	C	81.9%	803.4 / 166.7	371	Y212C (TAC→TGC)	yihQ	putative glycosidase
Reads supporting (aligned to +/- strand): ref base T (32/35); new base C (154/150); total (186/185)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 372 total aligned reads displayed.

GCAGGGCGGTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTAGGCTGTGGGA  >  NC_000913/4066130‑4066312

  agGGCGGTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTAt                                                                                          >  2:3311801/1‑93 (MQ=255)

   gggCGGTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCAtt                                                                                         >  1:6891609/1‑93 (MQ=255)

   gggCGGTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCAtt                                                                                         >  2:3133261/1‑93 (MQ=255)

   gggCGGTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCAtt                                                                                         <  2:1005061/93‑1 (MQ=255)

    ggCGGTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTc                                                                                        <  2:1204417/93‑1 (MQ=255)

     gCGGTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTcc                                                                                       >  2:1252084/1‑93 (MQ=255)

      cGGTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCg                                                                                      >  2:6282034/1‑93 (MQ=255)

       ggTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCgg                                                                                     <  1:773273/93‑1 (MQ=255)

       ggTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCgg                                                                                     <  1:3546370/93‑1 (MQ=255)

       ggTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCgg                                                                                     >  1:2566639/1‑93 (MQ=255)

        gTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCggg                                                                                    >  2:2466406/1‑93 (MQ=255)

        gTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCggg                                                                                    <  1:5318575/93‑1 (MQ=255)

        gTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCggg                                                                                    >  1:1567003/1‑93 (MQ=255)

           aaTTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCa                                                                                 >  2:2712278/1‑93 (MQ=255)

             tttttCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACt                                                                               <  2:2572258/93‑1 (MQ=255)

              ttttCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTa                                                                              >  2:5678792/1‑93 (MQ=255)

                ttCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTaaa                                                                            <  1:5682071/93‑1 (MQ=255)

                   cagcagGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTc                                                                         >  2:2270069/1‑93 (MQ=255)

                    agcagGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCg                                                                        >  1:3848510/1‑93 (MQ=255)

                     gcagGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGa                                                                       >  2:1842146/1‑93 (MQ=255)

                      cagGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGaa                                                                      <  1:549975/93‑1 (MQ=255)

                       agGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAg                                                                     >  1:1747900/1‑93 (MQ=255)

                          gAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTc                                                                  >  1:4167863/1‑93 (MQ=255)

                           aaaTGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCa                                                                 <  1:5856088/93‑1 (MQ=255)

                            aaTGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCat                                                                >  2:6494643/1‑93 (MQ=255)

                             atgtatgtGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCata                                                               <  2:7272320/93‑1 (MQ=255)

                               gtatgtGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCAtata                                                             <  2:5027869/93‑1 (MQ=255)

                                tatgtGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAg                                                            <  1:7322696/93‑1 (MQ=255)

                                 atgtGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGc                                                           >  1:58080/1‑93 (MQ=255)

                                 atgtGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGACATTCCGGGGCACTAAAGTCGAAGTTCATATAGc                                                           <  1:4742583/93‑1 (MQ=255)

                                  tgtgTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCa                                                          >  2:6181506/1‑93 (MQ=255)

                                  tgtgTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCa                                                          <  1:1521320/93‑1 (MQ=255)

                                    tgtCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAAc                                                        <  1:4274415/93‑1 (MQ=255)

                                        aGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGtt                                                    >  1:7496154/1‑93 (MQ=255)

                                         gCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTa                                                   <  2:6640044/93‑1 (MQ=255)

                                         gCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTa                                                   >  1:949761/1‑93 (MQ=255)

                                         gCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTa                                                   >  2:423508/1‑93 (MQ=255)

                                         gCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTa                                                   >  1:2385512/1‑93 (MQ=255)

                                          cacaTTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTAt                                                  <  1:6949381/93‑1 (MQ=255)

                                          cacaTTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTAt                                                  <  1:4314780/93‑1 (MQ=255)

                                           acaTTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATc                                                 <  1:4955835/93‑1 (MQ=255)

                                            caTTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCa                                                <  2:6017207/93‑1 (MQ=255)

                                             aTTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCaa                                               >  1:7055683/1‑93 (MQ=255)

                                              ttCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAAc                                              >  2:4342518/1‑93 (MQ=255)

                                               tCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACa                                             <  2:3077823/93‑1 (MQ=255)

                                                cAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACAt                                            <  1:5465221/93‑1 (MQ=255)

                                                 aaaaCGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATg                                           <  2:3907507/93‑1 (MQ=255)

                                                 aaaaCGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTAGAAGTTCATATAGCAACTGTTATCAACATg                                           <  1:6562055/93‑1 (MQ=255)

                                                  aaaCGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATgg                                          <  2:4033658/93‑1 (MQ=255)

                                                   aaCGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGc                                         <  2:4453997/93‑1 (MQ=255)

                                                   aaCGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGc                                         <  2:6360765/93‑1 (MQ=255)

                                                    aCGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCa                                        <  1:3311801/93‑1 (MQ=255)

                                                    aCGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCa                                        <  1:1150087/93‑1 (MQ=255)

                                                    aCGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCa                                        >  2:3731531/1‑93 (MQ=255)

                                                     cGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAg                                       >  2:342951/1‑93 (MQ=255)

                                                        aGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTaa                                    <  1:1243527/93‑1 (MQ=255)

                                                         gCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAAt                                   >  1:6659799/1‑93 (MQ=255)

                                                         gCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAAt                                   >  2:1385947/1‑93 (MQ=255)

                                                              gCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTc                              <  2:2764588/93‑1 (MQ=255)

                                                              gCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTc                              <  2:3848473/93‑1 (MQ=255)

                                                              gCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTc                              >  2:4672843/1‑93 (MQ=255)

                                                              gCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTc                              >  2:469592/1‑93 (MQ=255)

                                                              gCTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTc                              >  1:2844253/1‑93 (MQ=255)

                                                               cTTTGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCt                             <  2:5811302/93‑1 (MQ=255)

                                                                tttGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTg                            >  2:5901604/1‑93 (MQ=255)

                                                                tttGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTg                            >  1:2055100/1‑93 (MQ=255)

                                                                tttGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTg                            >  1:5938423/1‑93 (MQ=255)

                                                                tttGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTg                            <  1:4460285/93‑1 (MQ=255)

                                                                 ttGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGc                           >  1:2373772/1‑93 (MQ=255)

                                                                 ttGTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGc                           <  1:1000273/93‑1 (MQ=255)

                                                                  tGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCg                          <  1:3824430/93‑1 (MQ=255)

                                                                   gTCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGt                         <  1:1029060/93‑1 (MQ=255)

                                                                    tCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTg                        >  2:3685011/1‑93 (MQ=255)

                                                                    tCTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTg                        >  2:528197/1‑93 (MQ=255)

                                                                     cTTCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGc                       >  2:892431/1‑93 (MQ=255)

                                                                      ttCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCt                      >  1:6790186/1‑93 (MQ=255)

                                                                      ttCCCACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCt                      >  1:3991778/1‑93 (MQ=255)

                                                                        cccACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGa                    <  1:1900929/93‑1 (MQ=255)

                                                                        cccACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGa                    <  2:5792941/93‑1 (MQ=255)

                                                                         ccACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGAc                   >  2:1068991/1‑93 (MQ=255)

                                                                         ccACAGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCCAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGAc                   >  1:6213405/1‑93 (MQ=255)

                                                                          cacaGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACa                  >  1:1851032/1‑93 (MQ=255)

                                                                           acaGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACaa                 >  1:4675415/1‑93 (MQ=255)

                                                                           acaGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACaa                 <  2:7022218/93‑1 (MQ=255)

                                                                            caGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACaaa                <  2:5517339/93‑1 (MQ=255)

                                                                             aGCGCCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAAc               <  2:3610837/93‑1 (MQ=255)

                                                                               cgcCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGt             <  1:6408429/93‑1 (MQ=255)

                                                                                gcCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTa            >  2:4128925/1‑93 (MQ=255)

                                                                                gcCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTa            >  2:5157137/1‑93 (MQ=255)

                                                                                gcCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTa            >  2:1682570/1‑93 (MQ=255)

                                                                                gcCAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTa            <  2:3560199/93‑1 (MQ=255)

                                                                                 ccAGTTCATGGCATTCCGGGGCACTAAAGTGGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTAg           <  2:7434032/93‑1 (MQ=255)

                                                                                 ccAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTAg           >  2:6918295/1‑93 (MQ=255)

                                                                                  cAGTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTAgg          <  1:7231243/93‑1 (MQ=255)

                                                                                    gTTCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTAGGCt        <  1:362394/93‑1 (MQ=255)

                                                                                     ttCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTAGGCtg       >  1:3398914/1‑93 (MQ=255)

                                                                                      tCATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTAGGCtgt      <  1:3105044/93‑1 (MQ=255)

                                                                                       cATGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTAGGCtgtg     <  1:3042051/93‑1 (MQ=255)

                                                                                        aTGGCATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTAGGCTGTgg    <  1:295406/93‑1 (MQ=255)

                                                                                         tGGCATTCCGGGGCACTAAAGTCGAAGATCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTAGGCTGTggg   >  2:6429520/1‑93 (MQ=255)

GCAGGGCGGTTAATTTTTCCAGCAGGGAAATGTATGTGTCAGCACATTCAAAACGCAGCGTTGCTTTGTCTTCCCACAGCGCCAGTTCATGGTATTCCGGGGCACTAAAGTCGAAGTTCATATAGCAACTGTTATCAACATGGCAGTAATACTTCTGCGTGCTGACAAACGTAGGCTGTGGGA  >  NC_000913/4066130‑4066312

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