BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	308,935	G→A	R106C (CGC→TGC)	yagY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	308,935	0	G	A	99.6%	861.1 / NA	246	R106C (CGC→TGC)	yagY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (114/131); minor base C (0/1); total (114/132)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 254 total aligned reads displayed.

GGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGC  >  NC_000913/308845‑309023

gggcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCAct                                                                                        <  1:806776/93‑1 (MQ=255)

 ggcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCActc                                                                                       >  2:2007311/1‑93 (MQ=255)

      aaGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCttttt                                                                                  >  1:707146/1‑93 (MQ=255)

      aaGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCttttt                                                                                  <  1:1510319/93‑1 (MQ=255)

      aaGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCttttt                                                                                  >  2:2293191/1‑93 (MQ=255)

      aaGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCttttt                                                                                  >  1:5135810/1‑93 (MQ=255)

       aGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTc                                                                                 >  2:2690939/1‑93 (MQ=255)

       aGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTc                                                                                 <  1:5764929/93‑1 (MQ=255)

       aGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTc                                                                                 <  2:2856944/93‑1 (MQ=255)

       aGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTc                                                                                 <  1:2206535/93‑1 (MQ=255)

        gCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCa                                                                                >  1:2533706/1‑93 (MQ=255)

        gCAGTCGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCa                                                                                <  1:3784704/93‑1 (MQ=255)

         cAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCAt                                                                               <  1:440225/93‑1 (MQ=255)

          aGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATc                                                                              >  2:440571/1‑93 (MQ=255)

           gTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCt                                                                             >  1:911292/1‑93 (MQ=255)

           gTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCt                                                                             <  2:339642/93‑1 (MQ=255)

           gTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCt                                                                             <  2:257058/93‑1 (MQ=255)

                 acaGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGt                                                                       >  2:745357/1‑93 (MQ=255)

                  caGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTc                                                                      <  1:6290751/93‑1 (MQ=255)

                  caGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTc                                                                      >  1:4930285/1‑93 (MQ=255)

                    gcgcTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTccc                                                                    <  2:3510304/93‑1 (MQ=255)

                    gcgcTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTccc                                                                    >  1:3486154/1‑93 (MQ=255)

                       cTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGAGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTta                                                                 <  1:2943577/93‑1 (MQ=255)

                        ttgggTTGGCATTATCGCGCTGCACATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtat                                                                <  2:360460/90‑1 (MQ=255)

                          cGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATag                                                              <  1:3301871/93‑1 (MQ=255)

                          cGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATag                                                              <  1:2107829/93‑1 (MQ=255)

                           ggTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATaga                                                             <  2:4370246/93‑1 (MQ=255)

                            gTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaa                                                            >  2:5409484/1‑93 (MQ=255)

                            gTTGGCATTATCGCGCTGCGCATCACTGAGAGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaa                                                            >  1:3400285/1‑93 (MQ=255)

                             gtGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCCCTCTTTTTCATCTGCGGGTCCCTTATagaag                                                           <  1:4077705/92‑1 (MQ=255)

                              tGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaaga                                                          <  1:4662973/93‑1 (MQ=255)

                               ggCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAgaagaa                                                         <  2:5892374/93‑1 (MQ=255)

                               ggCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAgaagaa                                                         >  2:1879029/1‑93 (MQ=255)

                                gCATTATCGCGCTGCGCATCACTGAGGGCCTCATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGaagaag                                                        >  1:5112029/1‑93 (MQ=255)

                                gCATTATCGCGCTGCGCATCACTGAGGGCCGATTCAAACCAAACAATGCGGTATTACAACCTTTTTTCATCTGCGGTCCCCTTATAAaagaag                                                        >  2:3940072/1‑93 (MQ=255)

                                  aTTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCg                                                      <  1:1873685/93‑1 (MQ=255)

                                   ttATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCgg                                                     <  2:5119569/93‑1 (MQ=255)

                                    tATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                    <  1:5032183/93‑1 (MQ=255)

                                    tATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                    >  2:3913737/1‑93 (MQ=255)

                                    tATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                    >  2:519029/1‑93 (MQ=255)

                                    tATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                    >  1:5551976/1‑93 (MQ=255)

                                    tATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                    <  2:4496760/93‑1 (MQ=255)

                                     aTCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGAt                                                   <  1:1256218/93‑1 (MQ=255)

                                      tCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATc                                                  >  1:2007713/1‑93 (MQ=255)

                                        gcgcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCAc                                                >  2:662133/1‑93 (MQ=255)

                                            tGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCg                                            >  2:6196187/1‑93 (MQ=255)

                                            tGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCg                                            <  2:2445363/93‑1 (MQ=255)

                                            tGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCg                                            >  2:1830477/1‑93 (MQ=255)

                                             gcgcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGc                                           >  2:4079240/1‑93 (MQ=255)

                                               gcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTg                                         <  2:6088501/93‑1 (MQ=255)

                                                cATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTgg                                        >  1:588688/1‑93 (MQ=255)

                                                cATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTgg                                        >  1:3135411/1‑93 (MQ=255)

                                                  tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                      >  2:6234535/1‑93 (MQ=255)

                                                  tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                      >  1:3673651/1‑93 (MQ=255)

                                                  tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                      >  1:1308615/1‑93 (MQ=255)

                                                    aCTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTg                                    >  2:5136408/1‑93 (MQ=255)

                                                     cgtaGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTACAGAAGAAGCGGATCACTTCGCTGGCTTgg                                   <  1:1699939/90‑1 (MQ=255)

                                                     cTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTgg                                   >  1:987102/1‑93 (MQ=255)

                                                     cTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTgg                                   <  2:550106/93‑1 (MQ=255)

                                                      tGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTggg                                  <  2:1332946/93‑1 (MQ=255)

                                                      tGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTggg                                  <  1:3049305/93‑1 (MQ=255)

                                                      tGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTggg                                  <  1:1797899/93‑1 (MQ=255)

                                                       gAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTCTAGAAGCAGCGGATCACTTCGCTGGCTTGGGc                                 >  1:825996/1‑93 (MQ=255)

                                                        aGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCg                                <  1:1099041/93‑1 (MQ=255)

                                                         gggCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCgg                               >  2:1084630/1‑93 (MQ=255)

                                                         gggCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCgg                               <  1:4727687/93‑1 (MQ=255)

                                                          ggCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCggg                              >  2:5366452/1‑93 (MQ=255)

                                                           gCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGt                             <  2:2094474/93‑1 (MQ=255)

                                                           gCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGt                             <  2:1445151/93‑1 (MQ=255)

                                                            ccTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTa                            >  2:2809846/1‑93 (MQ=255)

                                                            ccTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTa                            <  2:2972764/93‑1 (MQ=255)

                                                            ccTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTa                            >  2:4888307/1‑93 (MQ=255)

                                                            ccTGATCAAACCAGACAATGCGGTAGTAGCAATCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTa                            >  1:2887325/1‑93 (MQ=255)

                                                               gATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagca                         <  1:5169808/93‑1 (MQ=255)

                                                                aTCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagcag                        >  1:2345236/1‑93 (MQ=255)

                                                                aTCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagcag                        >  2:4207239/1‑93 (MQ=255)

                                                                 tCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagcagc                       <  1:2649597/93‑1 (MQ=255)

                                                                  cAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAgcagca                      >  2:1706161/1‑93 (MQ=255)

                                                                   aaaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCaa                     <  2:5582049/93‑1 (MQ=255)

                                                                   aaaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCaa                     >  2:3704517/1‑93 (MQ=255)

                                                                    aaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAg                    >  1:3202701/1‑93 (MQ=255)

                                                                       cAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCTGGTAGCAGCAAGCTg                 <  1:1934111/93‑1 (MQ=255)

                                                                       cAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTg                 >  1:3016745/1‑93 (MQ=255)

                                                                       cAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTg                 <  1:3966690/93‑1 (MQ=255)

                                                                       cAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTg                 <  1:547190/93‑1 (MQ=255)

                                                                        aGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTgg                >  1:479107/1‑93 (MQ=255)

                                                                        aGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTgg                <  2:3108961/93‑1 (MQ=255)

                                                                         gACGATGCGGTAGTAGTACTGTTTTACATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGc               <  1:2314310/93‑1 (MQ=255)

                                                                            aaTGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCggg            <  1:1995346/93‑1 (MQ=255)

                                                                              tGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgag          <  2:941807/93‑1 (MQ=255)

                                                                              tGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgag          <  2:4828385/93‑1 (MQ=255)

                                                                               gCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagt         <  2:6010881/93‑1 (MQ=255)

                                                                               gCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagt         <  2:2707587/93‑1 (MQ=255)

                                                                                 ggTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtga       <  1:609701/93‑1 (MQ=255)

                                                                                 ggTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtga       <  1:1179893/93‑1 (MQ=255)

                                                                                  gtagtaGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtgag      <  2:5701692/93‑1 (MQ=255)

                                                                                  gtagtaGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtgag      >  2:1313937/1‑93 (MQ=255)

                                                                                  gtagtaGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtgag      <  1:565797/93‑1 (MQ=255)

                                                                                   tagtagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtgagt     >  2:5645120/1‑93 (MQ=255)

                                                                                      tagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGCGAGTGagcagc  >  2:5246490/1‑93 (MQ=255)

GGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGC  >  NC_000913/308845‑309023

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