BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,005,601	+C	coding (237/282 nt)	b4568 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,005,600	1	.	C	99.5%	708.3 / NA	219	G79A (GGC→GCC)	b4568	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); major base C (110/108); minor base A (1/0); total (111/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

CGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑C‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTG  >  NC_000913/4005508‑4005687

 gTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑c                                                                                         >  1:1341951/1‑93 (MQ=255)

  tAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cc                                                                                        >  1:634656/1‑92 (MQ=255)

  tAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cc                                                                                        >  1:442894/1‑92 (MQ=255)

  tAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cc                                                                                        >  1:3660452/1‑92 (MQ=255)

   aaTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cct                                                                                       >  2:5748121/1‑91 (MQ=255)

   aaTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cct                                                                                       <  1:216604/93‑3 (MQ=255)

   aaTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cct                                                                                       >  1:5540846/1‑91 (MQ=255)

    aTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cctg                                                                                      >  2:1863312/1‑90 (MQ=255)

    aTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cctg                                                                                      >  1:4325436/1‑90 (MQ=255)

       gTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTAt                                                                                   <  1:1847460/93‑1 (MQ=255)

       gTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTAt                                                                                   <  2:2618919/93‑1 (MQ=255)

        tttGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTAtc                                                                                  <  2:1330187/93‑1 (MQ=255)

          tGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTAtctc                                                                                <  1:3492736/93‑1 (MQ=255)

           gggCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATctct                                                                               >  1:5265329/1‑93 (MQ=255)

           gggCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATctct                                                                               <  2:5155917/93‑1 (MQ=255)

            ggCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTctg                                                                              >  2:1140930/1‑93 (MQ=255)

               tGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTctgctg                                                                           <  2:6229343/93‑1 (MQ=255)

               tGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTctgctg                                                                           <  1:238602/93‑1 (MQ=255)

                gggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTgg                                                                          >  1:1529381/1‑93 (MQ=255)

                gggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTgg                                                                          >  2:2435372/1‑93 (MQ=255)

                 ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGa                                                                         <  1:4009842/93‑1 (MQ=255)

                 ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGa                                                                         <  2:1073930/93‑1 (MQ=255)

                 ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGa                                                                         <  1:4559867/93‑1 (MQ=255)

                   aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATg                                                                       >  1:6166805/1‑93 (MQ=255)

                   aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATg                                                                       <  2:756596/93‑1 (MQ=255)

                    ttGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATgg                                                                      <  2:3324744/93‑1 (MQ=255)

                    ttGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATgg                                                                      <  2:5689104/93‑1 (MQ=255)

                     tGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATggg                                                                     <  2:1790861/93‑1 (MQ=255)

                      gcgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGt                                                                    <  2:1967522/93‑1 (MQ=255)

                       cgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGtt                                                                   <  1:935065/93‑1 (MQ=255)

                        gctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTa                                                                  >  1:3806642/1‑93 (MQ=255)

                        gctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTa                                                                  <  2:5282353/93‑1 (MQ=255)

                          tgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTAcc                                                                <  1:831595/93‑1 (MQ=255)

                           gctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCa                                                               <  1:2287912/93‑1 (MQ=255)

                            ctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAg                                                              >  2:2845812/1‑93 (MQ=255)

                             ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGa                                                             >  2:1670209/1‑93 (MQ=255)

                              tgcccTCCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTAATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGAt                                                            <  2:3939060/90‑1 (MQ=255)

                              tGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGAt                                                            <  2:42954/93‑1 (MQ=255)

                                gCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGc                                                          >  1:5577719/1‑93 (MQ=255)

                                gCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGc                                                          >  1:1190606/1‑93 (MQ=255)

                                 ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGCGGGAGGTGGCC‑TGTATATCTGCTGGATGGGTTACCAGATGCt                                                         >  1:4826556/1‑93 (MQ=255)

                                  cTGCATTTGATTATCGACAAAATGGCCTGGCTGCATCCGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCtc                                                        <  1:2806720/93‑2 (MQ=255)

                                  cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTa                                                        <  1:1426958/93‑1 (MQ=255)

                                   tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                       >  2:446659/1‑93 (MQ=255)

                                    gCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACg                                                      >  2:214576/1‑93 (MQ=255)

                                    gCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGGTACCAGATGCTACg                                                      >  1:1988394/1‑93 (MQ=255)

                                     cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACgt                                                     >  1:1072749/1‑93 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                    >  1:1026227/1‑93 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                    >  1:3987734/1‑93 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                    >  2:5809087/1‑93 (MQ=255)

                                       tttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACgtgg                                                   >  1:4531708/1‑93 (MQ=255)

                                       tttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACgtgg                                                   >  2:587370/1‑93 (MQ=255)

                                        ttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggt                                                  <  1:2592549/93‑1 (MQ=255)

                                         tGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggtg                                                 >  2:6267960/1‑93 (MQ=255)

                                         tGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggtg                                                 <  1:4726465/93‑1 (MQ=255)

                                         tGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggtg                                                 >  2:4016817/1‑93 (MQ=255)

                                          gATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGc                                                >  1:628060/1‑93 (MQ=255)

                                           aTTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCa                                               <  1:5650227/93‑1 (MQ=255)

                                             tATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACt                                             >  1:5367102/1‑93 (MQ=255)

                                             tATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACt                                             >  2:3519478/1‑93 (MQ=255)

                                               tCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGa                                           >  1:2925339/1‑93 (MQ=255)

                                               tCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGa                                           >  2:265231/1‑93 (MQ=255)

                                                  aaaaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACAAGATGCTACGTGGTGCACTGaaaa                                        >  1:1666682/1‑93 (MQ=255)

                                                         ggCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGc                                 <  1:5009499/93‑1 (MQ=255)

                                                         ggCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGc                                 <  1:4676475/93‑1 (MQ=255)

                                                           ccTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCgg                               >  1:3912080/1‑93 (MQ=255)

                                                            cTGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGt                              >  1:5474315/1‑93 (MQ=255)

                                                             tGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGtt                             <  1:338267/93‑1 (MQ=255)

                                                             tGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGtt                             >  1:5755250/1‑93 (MQ=255)

                                                             tGGCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGtt                             >  1:4108728/1‑93 (MQ=255)

                                                               gCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTc                           >  2:5718454/1‑93 (MQ=255)

                                                               gCTGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTc                           <  2:3362624/93‑1 (MQ=255)

                                                               gCTGCATACGCTGATTATGGTGGGCGGTGG‑CGTGTATATCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACGGAAAAAAGAGGCGGTTTc                           <  2:512355/93‑1 (MQ=255)

                                                                 tGCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGGTTCTg                         >  2:2552065/1‑93 (MQ=255)

                                                                  gCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGc                        <  1:5053530/93‑1 (MQ=255)

                                                                  gCATACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGc                        >  1:593168/1‑93 (MQ=255)

                                                                   cATACGCTGA‑TATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAc                      <  2:407029/93‑1 (MQ=255)

                                                                    aTACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAc                      >  2:5590420/1‑93 (MQ=255)

                                                                     tACGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGACGGTTTCTGCAcc                     <  1:3201364/93‑1 (MQ=255)

                                                                     tACGCTGATTATGGTGGGCGGTGGAC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                     >  2:908798/1‑93 (MQ=255)

                                                                      aCGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCt                    <  1:5807296/93‑1 (MQ=255)

                                                                       cGCTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTg                   >  2:3055017/1‑93 (MQ=255)

                                                                         cTGATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgcg                 >  1:5619764/1‑93 (MQ=255)

                                                                           gATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGcc               >  1:6089398/1‑93 (MQ=255)

                                                                           gATTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGcc               >  2:3603292/1‑93 (MQ=255)

                                                                            aTTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCca              >  2:3776431/1‑93 (MQ=255)

                                                                            aTTATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCca              >  2:4244717/1‑93 (MQ=255)

                                                                              tATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcaca            <  2:2506464/93‑1 (MQ=255)

                                                                              tATGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcaca            <  2:4878539/93‑1 (MQ=255)

                                                                               aTGGTGGGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAg           >  1:43975/1‑93 (MQ=255)

                                                                                 ggtggGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGt         <  1:927509/93‑1 (MQ=255)

                                                                                  gtggGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTc        <  2:6085691/93‑1 (MQ=255)

                                                                                   tggGCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCg       >  1:6207584/1‑93 (MQ=255)

                                                                                    gggCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGa      >  2:5380798/1‑93 (MQ=255)

                                                                                    gggCGGTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGa      <  2:1392668/93‑1 (MQ=255)

                                                                                      gCGGTGGCC‑TGTATCTCTGCTGGATGGGTTCCCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGc    <  2:2297471/93‑1 (MQ=255)

                                                                                        ggTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTg  >  2:2861649/1‑93 (MQ=255)

                                                                                        ggTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTg  >  1:4276350/1‑93 (MQ=255)

                                                                                        ggTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTg  >  1:252211/1‑93 (MQ=255)

                                                                                        ggTGGCC‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTg  >  2:495453/1‑93 (MQ=255)

CGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑C‑TGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTG  >  NC_000913/4005508‑4005687

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: `‑`