BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,463,677	C→T	G166G (GGC→GGT)	hofD →	leader peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,463,677	0	C	T	99.4%	1089.6 / NA	314	G166G (GGC→GGT)	hofD	leader peptidase
Reads supporting (aligned to +/- strand): ref base C (0/1); major base T (152/158); minor base G (2/0); total (155/159)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

GGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCTGCA  >  NC_000913/3463587‑3463766

ggTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGgtgc                                                                                         <  1:6396801/93‑1 (MQ=255)

  tGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAATTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCAt                                                                                       <  1:7084607/93‑1 (MQ=255)

   gCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATg                                                                                      <  1:4775073/93‑1 (MQ=255)

   gCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATg                                                                                      >  2:3425728/1‑93 (MQ=255)

     tgtgGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATggg                                                                                    >  2:7190575/1‑92 (MQ=255)

     tgtgGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGgt                                                                                    >  1:1817196/1‑93 (MQ=255)

        ggtgTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGtgtg                                                                                 >  2:5449407/1‑93 (MQ=255)

        ggtgTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGtgtg                                                                                 <  2:6330760/93‑1 (MQ=255)

          tgtgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggc                                                                               >  1:7649197/1‑93 (MQ=255)

          tgtgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggc                                                                               >  1:7502271/1‑93 (MQ=255)

          tgtgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggc                                                                               >  2:7337831/1‑93 (MQ=255)

            tgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggctg                                                                             <  1:5419880/93‑1 (MQ=255)

             gtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggctgg                                                                            >  2:4112983/1‑93 (MQ=255)

               ttACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCa                                                                          <  1:1239044/93‑1 (MQ=255)

                 aCTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCaaa                                                                        <  1:1570630/93‑1 (MQ=255)

                 aCTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCaaa                                                                        <  1:5106535/93‑1 (MQ=255)

                   tGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGGAAACg                                                                      >  2:4712123/1‑93 (MQ=255)

                   tGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACg                                                                      <  1:3392483/93‑1 (MQ=255)

                     gggCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGtt                                                                    <  1:2823291/93‑1 (MQ=255)

                     gggCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGtt                                                                    <  1:5655846/93‑1 (MQ=255)

                       gCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCGAACGTTGc                                                                  >  2:2282491/1‑93 (MQ=255)

                        cGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGcc                                                                 <  2:1459825/93‑1 (MQ=255)

                        cGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGcc                                                                 >  2:497390/1‑93 (MQ=255)

                          tCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAAAGTTGccca                                                               >  2:4615581/1‑91 (MQ=255)

                          tCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGGGCATGGTGTGGCTGGCAAACGTTGCCaa                                                               >  1:3176250/1‑93 (MQ=255)

                             ggTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGa                                                            >  1:498193/1‑93 (MQ=255)

                             ggTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGa                                                            <  1:2368871/93‑1 (MQ=255)

                               ttAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATa                                                          <  2:5688671/93‑1 (MQ=255)

                                tAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTTCATGGTGTGGCTGGCAAACGTTGCCAATGATAc                                                         <  2:3113631/93‑1 (MQ=255)

                                tAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATAc                                                         >  1:4265431/1‑93 (MQ=255)

                                tAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATAc                                                         >  1:3917751/1‑93 (MQ=255)

                                tAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATAc                                                         <  2:6102366/93‑1 (MQ=255)

                                tAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATAc                                                         >  1:7248558/1‑93 (MQ=255)

                                tAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATAc                                                         <  2:5119315/93‑1 (MQ=255)

                                 aGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGGTGCCCATGATAct                                                        >  2:1631633/1‑93 (MQ=255)

                                  gTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATActg                                                       <  1:5928621/93‑1 (MQ=255)

                                  gTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATActg                                                       >  1:4409645/1‑93 (MQ=255)

                                      gTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGa                                                   <  2:2104723/93‑1 (MQ=255)

                                      gTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGa                                                   <  1:5569351/93‑1 (MQ=255)

                                       tCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGAt                                                  >  2:3252123/1‑93 (MQ=255)

                                       tCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGGTGCCAATGATACTGCTGAt                                                  >  1:7395564/1‑93 (MQ=255)

                                        cacaAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGAtt                                                 <  1:3116096/93‑1 (MQ=255)

                                        cacaAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGAtt                                                 <  1:3609330/93‑1 (MQ=255)

                                        cacaAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGAtt                                                 <  2:371626/93‑1 (MQ=255)

                                         acaAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTg                                                >  1:1303482/1‑93 (MQ=255)

                                          caAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGc                                               <  1:3208647/93‑1 (MQ=255)

                                            aaGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCt                                             >  2:6954641/1‑93 (MQ=255)

                                            aaGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCt                                             >  1:372695/1‑93 (MQ=255)

                                             agagGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTc                                            >  2:128713/1‑93 (MQ=255)

                                              gagGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCg                                           >  2:5098852/1‑93 (MQ=255)

                                                  gATTGGGCTACGGTGATTTCAAGTTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGAAAACGTTGCCAATGATACTGCTGATTGCCTCGCTgg                                       <  1:701147/93‑1 (MQ=255)

                                                   aTTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCCCGCTggg                                      >  1:5716742/1‑93 (MQ=255)

                                                    ttGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCCATGATACTGCTGATTGCCTCGCTgggg                                     >  2:3557240/1‑92 (MQ=255)

                                                    ttGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGt                                     <  2:6663938/93‑1 (MQ=255)

                                                      gggCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTgggggg                                   >  2:852855/1‑93 (MQ=255)

                                                       ggCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGc                                  <  1:5299207/93‑1 (MQ=255)

                                                          tACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCAtt                               >  1:4121378/1‑93 (MQ=255)

                                                           aCGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGGAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGTATTg                              >  1:6316826/1‑93 (MQ=255)

                                                           aCGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTg                              <  1:6626723/93‑1 (MQ=255)

                                                           aCGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTg                              <  1:811694/93‑1 (MQ=255)

                                                           aCGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTg                              >  1:2897288/1‑93 (MQ=255)

                                                           aCGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTg                              <  2:1851203/93‑1 (MQ=255)

                                                           aCGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTg                              >  2:2974382/1‑93 (MQ=255)

                                                              gTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGtt                           <  2:4586651/93‑1 (MQ=255)

                                                               tGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTa                          >  1:5244850/1‑93 (MQ=255)

                                                               tGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTa                          <  2:7872033/93‑1 (MQ=255)

                                                                gATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTAc                         >  2:3304973/1‑93 (MQ=255)

                                                                gATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTAc                         <  2:6723107/93‑1 (MQ=255)

                                                                gATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTAc                         >  1:2917830/1‑93 (MQ=255)

                                                                 aTTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACg                        >  1:1349622/1‑93 (MQ=255)

                                                                   ttCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGcc                      <  2:7604147/93‑1 (MQ=255)

                                                                    tCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCa                     <  2:6526505/93‑1 (MQ=255)

                                                                    tCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCa                     >  2:638332/1‑93 (MQ=255)

                                                                    tCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCa                     >  2:767044/1‑93 (MQ=255)

                                                                     cAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCAt                    <  2:5988226/93‑1 (MQ=255)

                                                                     cAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCAt                    <  2:2373141/93‑1 (MQ=255)

                                                                      aaGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATc                   >  2:2602916/1‑93 (MQ=255)

                                                                      aaGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATc                   <  2:1314182/93‑1 (MQ=255)

                                                                      aaGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATc                   >  2:4301612/1‑93 (MQ=255)

                                                                      aaGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATc                   >  1:2290162/1‑93 (MQ=255)

                                                                        gCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGt                 <  2:5613234/93‑1 (MQ=255)

                                                                        acTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGt                 <  2:1799604/92‑1 (MQ=255)

                                                                          tACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGttt               <  1:948532/93‑1 (MQ=255)

                                                                           aCTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTc              >  2:2453405/1‑93 (MQ=255)

                                                                             tGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTcac            >  2:5942306/1‑93 (MQ=255)

                                                                             aggCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTcac            >  1:2744988/2‑93 (MQ=255)

                                                                              ggCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTcaca           <  2:1918188/93‑1 (MQ=255)

                                                                               gCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACaa          <  2:7437972/93‑1 (MQ=255)

                                                                               gCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACaa          <  1:3696192/93‑1 (MQ=255)

                                                                                ccGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAAc         >  1:5853688/1‑93 (MQ=255)

                                                                                  gCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACtt       <  2:1107938/93‑1 (MQ=255)

                                                                                   cTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTc      <  2:7741600/93‑1 (MQ=255)

                                                                                    tGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCt     <  2:5852218/93‑1 (MQ=255)

                                                                                     gCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGtt                <  2:6806393/81‑1 (MQ=255)

                                                                                     gCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCtc    >  1:6715864/1‑92 (MQ=255)

                                                                                     gCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCTg    >  1:1455666/1‑93 (MQ=255)

                                                                                      cAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCTGc   >  2:281329/1‑93 (MQ=255)

                                                                                      cAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCTGc   <  2:533237/93‑1 (MQ=255)

                                                                                       aGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCTGCa  <  2:5782592/93‑1 (MQ=255)

                                                                                       aGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCTGCa  >  2:4395599/1‑93 (MQ=255)

GGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCTGCA  >  NC_000913/3463587‑3463766

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