BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,506,868	A→G	P193P (CCT→CCC)	glk ←	glucokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,506,868	0	A	G	99.3%	923.2 / ‑3.4	276	P193P (CCT→CCC)	glk	glucokinase
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (119/155); total (121/155)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

Only 100 of 305 total aligned reads displayed.

CAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGC  >  NC_000913/2506778‑2506957

cagcGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCggg                                                                                         <  1:5890391/93‑1 (MQ=255)

cagcGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCggg                                                                                         >  2:7711099/1‑93 (MQ=255)

 agcGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGTGCGCTACAAATTCACCAGCCCggga                                                                                        >  2:7940607/1‑92 (MQ=255)

  gcgcgcGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGc                                                                                       >  1:5822838/1‑93 (MQ=255)

   cgcgcgTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGcc                                                                                      >  1:549753/1‑93 (MQ=255)

   cgcgcgTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGcc                                                                                      <  2:2814941/93‑1 (MQ=255)

      gcgTTCGGTAATATCTTTTGGCTTGATATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGa                                                                                   >  1:2603377/1‑93 (MQ=255)

        gTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGaaa                                                                                 >  2:8458508/1‑93 (MQ=255)

           cGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCa                                                                              <  2:2865391/93‑1 (MQ=255)

           cGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCa                                                                              >  1:7562119/1‑93 (MQ=255)

            ggTAATATCTTTTGGCTGGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATGGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAc                                                                             <  1:8358207/93‑1 (MQ=255)

             gTATTATCTTTTGGCTTGAGATTTTCTGAAAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTGAACACCCCCGCGCCAGAAAGCAcg                                                                            >  2:6173306/1‑93 (MQ=255)

             gTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAcg                                                                            >  1:8171376/1‑93 (MQ=255)

             gTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAcg                                                                            >  2:7417313/1‑93 (MQ=255)

              tAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAcgc                                                                           >  2:8329539/1‑93 (MQ=255)

              tAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAcgc                                                                           >  1:8540497/1‑93 (MQ=255)

               aaTATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAcgcg                                                                          <  1:2099639/93‑1 (MQ=255)

               aaTATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAcgcg                                                                          <  2:7598935/93‑1 (MQ=255)

                atatCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACgcgc                                                                         <  1:7971894/93‑1 (MQ=255)

                atatCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACgcgc                                                                         >  1:2325221/1‑93 (MQ=255)

                atatCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACgcgc                                                                         >  2:4430208/1‑93 (MQ=255)

                 tatCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCAcgcgat                                                                        >  1:4064194/1‑91 (MQ=255)

                 tatCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCt                                                                        >  2:5389467/1‑93 (MQ=255)

                  atCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTc                                                                       <  1:7617822/93‑1 (MQ=255)

                  atCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTc                                                                       <  2:1560745/93‑1 (MQ=255)

                  atCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTc                                                                       <  1:2263004/93‑1 (MQ=255)

                     ttttGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTccgc                                                                    <  2:317185/93‑1 (MQ=255)

                      tttGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTccgcc                                                                   >  2:8734081/1‑93 (MQ=255)

                       ttGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTccgccg                                                                  <  2:1028294/93‑1 (MQ=255)

                        tGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGa                                                                 >  2:730763/1‑93 (MQ=255)

                        tGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGa                                                                 <  2:2484894/93‑1 (MQ=255)

                         ggCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTccgccgca                                                                >  2:3681163/1‑91 (MQ=255)

                         ggCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGaa                                                                >  1:4079055/1‑93 (MQ=255)

                         ggCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGaa                                                                >  1:6855761/1‑93 (MQ=255)

                          gcgtGATATTGTATGGCAGGCTGTTGCCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGaaa                                                               <  1:4386940/90‑1 (MQ=255)

                          gCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGaaa                                                               <  1:3430552/93‑1 (MQ=255)

                            ttGAGATTTTCTTGCAGGGGGTTGTCAGCTTTCGCAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACa                                                             <  2:6838098/93‑1 (MQ=255)

                            ttGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACa                                                             >  1:6453053/1‑93 (MQ=255)

                            ttGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACa                                                             >  1:6058540/1‑93 (MQ=255)

                              gagaTTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATg                                                           <  2:7084274/93‑1 (MQ=255)

                                 gttttCTGGCAGGCGGTTGTCAGATTTCACAATTGCGCGATACAAATTCACCAGCCCGGGTCCAGAAAGCACGCGCTCCGCCTAAACATGAcc                                                        <  2:5998170/92‑1 (MQ=255)

                                 aTTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGAcc                                                        <  1:7431927/93‑1 (MQ=255)

                                   tttCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCaa                                                      >  1:327794/1‑93 (MQ=255)

                                    ttCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGTCGAAACATGACCAAt                                                     <  1:1046611/93‑1 (MQ=255)

                                    ttCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAAt                                                     <  1:122160/93‑1 (MQ=255)

                                    ttCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAAt                                                     >  2:6427656/1‑93 (MQ=255)

                                     tCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAAtt                                                    <  1:4762473/93‑1 (MQ=255)

                                      cTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAAttt                                                   >  2:560940/1‑93 (MQ=255)

                                        ggCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTcc                                                 <  2:3857901/93‑1 (MQ=255)

                                        ggCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTcc                                                 >  2:309732/1‑93 (MQ=255)

                                          cAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCcgc                                               <  1:4689280/93‑1 (MQ=255)

                                           aGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCcgca                                              <  2:8193130/93‑1 (MQ=255)

                                           aGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCcgca                                              <  1:2192503/93‑1 (MQ=255)

                                           aGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCcgca                                              <  1:6215274/93‑1 (MQ=255)

                                             gCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCcgcacg                                            <  1:6706727/93‑1 (MQ=255)

                                                gTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCACTTTCCGCACGCaa                                         >  1:1100769/1‑93 (MQ=255)

                                                gTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCaa                                         >  2:1242382/1‑93 (MQ=255)

                                                 ttGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAat                                        >  1:3498203/1‑93 (MQ=255)

                                                 ttGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAat                                        <  2:6501172/93‑1 (MQ=255)

                                                   gccAGCTTTCACAATTGCTCGATACAAATTTACCAGGCCAGGGCCAGAAAGAACGCGCTCCGCCGAAACATCACAAATTTCCGCTAGCAatat                                      >  2:1534093/3‑93 (MQ=255)

                                                    tCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATAtt                                     <  2:2147478/93‑1 (MQ=255)

                                                    tCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATAtt                                     >  1:7106/1‑93 (MQ=255)

                                                      aGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTc                                   <  2:5304473/93‑1 (MQ=255)

                                                      aGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTc                                   <  2:3311893/93‑1 (MQ=255)

                                                      aGCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTc                                   <  1:5108870/93‑1 (MQ=255)

                                                       gCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCg                                  >  2:2692451/1‑93 (MQ=255)

                                                       gCTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGATCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCg                                  >  2:4027679/1‑93 (MQ=255)

                                                       gCTTTCACAATTGCGCGATACAAATTCACCAGACCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCg                                  >  1:170747/1‑93 (MQ=255)

                                                        cTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGa                                 <  1:5621380/93‑1 (MQ=255)

                                                        cTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGa                                 <  2:2204351/93‑1 (MQ=255)

                                                        cTTTCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGa                                 <  2:5475727/93‑1 (MQ=255)

                                                           tCACAATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGGTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGa                              <  1:3575875/93‑1 (MQ=255)

                                                            cacaATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGAt                             >  2:1701897/1‑93 (MQ=255)

                                                             acaATTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATa                            >  1:3067709/1‑93 (MQ=255)

                                                               aaTTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAAt                          <  1:8655830/93‑1 (MQ=255)

                                                                agtGCGCGATTCAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAACTTCCGCACGCAATACTTCGAGGATAATg                         >  1:7800435/3‑93 (MQ=255)

                                                                aTTGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATg                         <  2:6946720/93‑1 (MQ=255)

                                                                  tGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGc                       >  1:8799415/1‑93 (MQ=255)

                                                                  tGCGCGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGc                       >  1:7279609/1‑93 (MQ=255)

                                                                   gcgcGATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGcc                      <  1:3024240/93‑1 (MQ=255)

                                                                     gcgATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCctc                    >  2:4365353/1‑93 (MQ=255)

                                                                     gcgATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCctc                    <  1:4031182/93‑1 (MQ=255)

                                                                     gcgATACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCctc                    <  1:5214216/93‑1 (MQ=255)

                                                                      cgATACAAAGTCACCAGCCCGGGGCCAGAAAGCACGCGCCCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCctct                   <  1:3843022/93‑1 (MQ=255)

                                                                         tACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCtcttct                <  1:2886664/93‑1 (MQ=255)

                                                                         tACAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCtcttct                <  1:5305249/93‑1 (MQ=255)

                                                                          aCAAATTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTcttctt               <  1:6354524/93‑1 (MQ=255)

                                                                           cAAATTCACCAGCCCGGTGACAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCttcttc              <  2:5167692/93‑1 (MQ=255)

                                                                             aaTTCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCAc            >  2:1774580/1‑93 (MQ=255)

                                                                               ttCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACtt          >  2:6597753/1‑92 (MQ=255)

                                                                               ttCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTa          <  1:1572499/93‑1 (MQ=255)

                                                                               atCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTa          <  2:7250806/92‑1 (MQ=255)

                                                                                tCACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTAt         >  2:467277/1‑93 (MQ=255)

                                                                                 cACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTAtt        >  2:2332949/1‑93 (MQ=255)

                                                                                 cACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTAtt        <  2:8103698/93‑1 (MQ=255)

                                                                                 cACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTAtt        <  1:2731297/93‑1 (MQ=255)

                                                                                 cACCAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTAtt        >  1:2539170/1‑93 (MQ=255)

                                                                                  cccggCCCGGGGCCAGAAAGCAAGCGTTGCGCCGAAACATCACCAATTTCCGCACGCCATATCTCGAGGATAATGGCCTCTTCTTCACTATTc       <  2:919201/89‑1 (MQ=255)

                                                                                  cccAGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTc       <  2:5109026/92‑1 (MQ=255)

                                                                                     aGCCCGGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGgc    >  1:1314173/1‑93 (MQ=255)

CAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTTCACAATTGCGCGATACAAATTCACCAGCCCAGGGCCAGAAAGCACGCGCTCCGCCGAAACATGACCAATTTCCGCACGCAATATTTCGAGGATAATGGCCTCTTCTTCACTATTCGGCGC  >  NC_000913/2506778‑2506957

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