BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,158,202	+C	coding (961/975 nt)	ygiRQ ←	Fe‑S oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,158,201	1	.	C	98.3%	806.7 / 4.1	230	coding (962/975 nt)	ygiRQ	Fe‑S oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (0/4); new base C (109/117); total (109/121)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

CGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACA  >  NC_000913/3158110‑3158288

 gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cc                                                                                       <  1:4039892/93‑2 (MQ=255)

 gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cc                                                                                       <  1:673722/93‑2 (MQ=255)

 gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cc                                                                                       >  1:18632/1‑92 (MQ=255)

 gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cc                                                                                       <  1:2115189/93‑2 (MQ=255)

  cTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cct                                                                                      <  2:2484430/93‑3 (MQ=255)

       cATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggtt                                                                                  <  2:2030422/93‑1 (MQ=255)

        aTCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttg                                                                                 <  2:1652330/93‑1 (MQ=255)

         tctcTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttgg                                                                                >  2:4786817/1‑93 (MQ=255)

         tctcTTCGGTAGAAAGCGGAATAGCAGGCGGGGTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGgggtgg                                                                                >  1:4436777/1‑93 (MQ=255)

          ctctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttggt                                                                               >  1:1994128/1‑93 (MQ=255)

          ctctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttggt                                                                               <  2:2906635/93‑1 (MQ=255)

           tctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttggtt                                                                              >  2:95713/1‑93 (MQ=255)

           tctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttggtt                                                                              <  1:877546/93‑1 (MQ=255)

             ttCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTc                                                                            >  2:4909706/1‑93 (MQ=255)

               cGGTAGAAAGCGGAATAGCAGGCGTGTTGATCCACACATAGCCGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgt                                                                          <  1:5209361/93‑1 (MQ=255)

               cGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgt                                                                          <  1:1015420/93‑1 (MQ=255)

                ggTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgtg                                                                         >  2:3441477/1‑93 (MQ=255)

                ggTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGGTGGTTTCgtg                                                                         >  1:4505604/1‑93 (MQ=255)

                 gTAGAAAGCGGAATAGCAGGCGGGTTGATCCATACATAGCGGTCGCCGTTTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgtgg                                                                        <  1:5384433/93‑1 (MQ=255)

                 gTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgtgg                                                                        <  2:5684917/93‑1 (MQ=255)

                  tAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgtggt                                                                       <  2:2880885/93‑1 (MQ=255)

                  tAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgtggt                                                                       >  1:935416/1‑93 (MQ=255)

                   aGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgtggtg                                                                      >  2:3259212/1‑93 (MQ=255)

                     aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCa                                                                    <  2:2667700/93‑1 (MQ=255)

                     aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCa                                                                    <  1:1723360/93‑1 (MQ=255)

                     aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCa                                                                    <  2:1268869/93‑1 (MQ=255)

                      aaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAg                                                                   <  1:4923560/93‑1 (MQ=255)

                        gCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGaa                                                                 <  1:5923646/93‑1 (MQ=255)

                        gCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGaa                                                                 >  1:5190952/1‑93 (MQ=255)

                        gCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGaa                                                                 >  1:753306/1‑93 (MQ=255)

                        gCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGaa                                                                 <  2:3037086/93‑1 (MQ=255)

                          ggAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATa                                                               <  1:1137827/93‑1 (MQ=255)

                          ggAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATa                                                               <  2:2168030/93‑1 (MQ=255)

                           gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAc                                                              >  1:3767150/1‑93 (MQ=255)

                             aTAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAcgc                                                            <  1:2273566/93‑1 (MQ=255)

                              tAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAcgcg                                                           >  2:4051712/1‑93 (MQ=255)

                              tAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAcgcg                                                           >  1:2476885/1‑93 (MQ=255)

                               aGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGa                                                          >  1:201888/1‑93 (MQ=255)

                                gCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGaa                                                         >  1:4830767/1‑93 (MQ=255)

                                gCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGaa                                                         >  2:5088126/1‑93 (MQ=255)

                                gCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGaa                                                         <  2:4246630/93‑1 (MQ=255)

                                gCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGaa                                                         <  1:5005681/93‑1 (MQ=255)

                                   ggCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                      >  1:5595142/1‑93 (MQ=255)

                                   ggCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                      <  2:746240/93‑1 (MQ=255)

                                   ggCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                      >  1:2616762/1‑93 (MQ=255)

                                   ggCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                      <  1:317677/93‑1 (MQ=255)

                                   ggCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                      <  1:966825/93‑1 (MQ=255)

                                    gCGGTTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCAt                                                     <  1:4729672/93‑1 (MQ=255)

                                    gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCAt                                                     <  1:5674279/93‑1 (MQ=255)

                                     cGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                                                    <  1:3774187/93‑1 (MQ=255)

                                     cGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                                                    <  1:4888443/93‑1 (MQ=255)

                                     cGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                                                    <  2:2610551/93‑1 (MQ=255)

                                     cGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                                                    <  2:1135114/93‑1 (MQ=255)

                                     cGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                                                    <  2:798793/93‑1 (MQ=255)

                                      gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                                                   >  1:3585879/1‑93 (MQ=255)

                                       ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                                  <  2:5282757/93‑1 (MQ=255)

                                       ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                                  >  1:399209/1‑93 (MQ=255)

                                        gTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGc                                                 >  1:5651985/1‑93 (MQ=255)

                                        gTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGc                                                 >  2:5995326/1‑93 (MQ=255)

                                        gTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGc                                                 >  2:5906059/1‑93 (MQ=255)

                                        gTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGc                                                 >  1:1370853/1‑93 (MQ=255)

                                         ttGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCg                                                >  1:551964/1‑93 (MQ=255)

                                         ttGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCg                                                >  2:2311269/1‑93 (MQ=255)

                                          tGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGt                                               <  1:3259126/93‑1 (MQ=255)

                                           gATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                                              >  2:3517745/1‑93 (MQ=255)

                                             tCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACa                                            >  2:2251973/1‑93 (MQ=255)

                                              ccACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAg                                           >  1:3546544/1‑93 (MQ=255)

                                               cacacaTAGCGTTCTCCGTGTTTTTGCATCAATGCGCGGGCACATC‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGc                                          <  1:3302621/93‑1 (MQ=255)

                                               cacacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGc                                          <  1:5341291/93‑1 (MQ=255)

                                                  acaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACt                                       >  1:2915389/1‑93 (MQ=255)

                                                    aTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACttt                                     >  2:5871699/1‑93 (MQ=255)

                                                    aTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACttt                                     <  1:2766392/93‑1 (MQ=255)

                                                     tAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTa                                    <  2:255310/93‑1 (MQ=255)

                                                      aGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTAt                                   <  1:5775500/93‑1 (MQ=255)

                                                           tCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGccc                              >  2:1342732/1‑93 (MQ=255)

                                                             gCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCtt                            >  2:1428542/1‑93 (MQ=255)

                                                              ccGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                           <  1:609103/93‑1 (MQ=255)

                                                               cGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCa                          >  1:1568277/1‑93 (MQ=255)

                                                               cGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGGTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTcc                          >  1:5085434/1‑92 (MQ=255)

                                                                 tgtTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACt                        >  1:4182940/1‑93 (MQ=255)

                                                                 tgtTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACt                        >  2:560974/1‑93 (MQ=255)

                                                                 tgtTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACt                        >  2:4999237/1‑93 (MQ=255)

                                                                    ttttGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACtttt                     >  2:452597/1‑92 (MQ=255)

                                                                      ttGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtctc                   <  2:3959582/93‑1 (MQ=255)

                                                                       tGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCg                  <  1:386737/93‑1 (MQ=255)

                                                                        gCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGa                 >  2:5793476/1‑93 (MQ=255)

                                                                             aaTGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa            >  2:4107518/1‑93 (MQ=255)

                                                                             aaTGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa            >  2:110698/1‑93 (MQ=255)

                                                                              aTGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaag           <  1:3212202/93‑1 (MQ=255)

                                                                              aTGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaag           <  1:6008534/93‑1 (MQ=255)

                                                                               tGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAgg          >  2:5210344/1‑93 (MQ=255)

                                                                                 cgcgGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgca        <  2:3661126/93‑1 (MQ=255)

                                                                                 cgcgGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgca        <  2:2995417/93‑1 (MQ=255)

                                                                                 cgcgGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgca        <  1:3039189/93‑1 (MQ=255)

                                                                                   ctggCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgcagc      <  2:1515013/91‑1 (MQ=255)

                                                                                   cgGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgcagc      <  2:3806599/93‑1 (MQ=255)

                                                                                    gggCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgcagca     >  2:1477847/1‑93 (MQ=255)

                                                                                     ggCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgcagcac    >  1:4556803/1‑92 (MQ=255)

                                                                                      gCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAac   <  2:4405574/93‑1 (MQ=255)

                                                                                       caca‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAaca  <  1:2476574/93‑1 (MQ=255)

CGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACA  >  NC_000913/3158110‑3158288

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