BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,158,827	G→T	G61V (GGC→GTC)	yijC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,158,827	0	G	T	99.6%	777.6 / NA	222	G61V (GGC→GTC)	yijC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base T (90/131); total (91/131)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 228 total aligned reads displayed.

GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCC  >  NC_000913/4158737‑4158916

gccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCa                                                                                         >  2:2157761/1‑93 (MQ=255)

 ccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCAt                                                                                        >  1:3386672/1‑93 (MQ=255)

 ccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCAt                                                                                        <  2:1354760/93‑1 (MQ=255)

   gTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTg                                                                                      <  1:4480048/93‑1 (MQ=255)

   gTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTg                                                                                      >  1:1278544/1‑93 (MQ=255)

   gTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTg                                                                                      >  1:956963/1‑93 (MQ=255)

    tgcGCTGGTGGATGCCGCATTGAGCCAATTAAGTGCTGAACGCAGCGTCGCCAGCCTGAGTTTGCGTGAGGTGGCGCGTGAAGCGGTCATTGc                                                                                     <  1:4363247/91‑1 (MQ=255)

      cGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTc                                                                                   >  2:2493152/1‑93 (MQ=255)

      cGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTc                                                                                   <  1:2240760/93‑1 (MQ=255)

        cTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTccc                                                                                 >  1:561713/1‑93 (MQ=255)

         tggtggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCa                                                                                <  1:2912542/93‑1 (MQ=255)

          ggtggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCAc                                                                               >  1:1199440/1‑93 (MQ=255)

          ggtggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCAc                                                                               >  1:1973652/1‑93 (MQ=255)

           gtggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCAcc                                                                              >  2:3437349/1‑93 (MQ=255)

            tggAAGCCTCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACct                                                                             <  2:130396/93‑1 (MQ=255)

                acctgCATTTAGCCAATTAAGTGCTGATCGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCttt                                                                         <  1:238442/89‑1 (MQ=255)

                aGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCttt                                                                         >  1:2615144/1‑93 (MQ=255)

                 gCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCtttt                                                                        <  1:4559205/93‑1 (MQ=255)

                 gCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCtttt                                                                        >  1:4126181/1‑93 (MQ=255)

                  ccGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCttttt                                                                       <  1:3077931/93‑1 (MQ=255)

                   cGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCGCCTCTTTTTa                                                                      <  2:2813519/93‑1 (MQ=255)

                   cGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTa                                                                      <  1:342959/93‑1 (MQ=255)

                    gCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTAt                                                                     >  1:604400/1‑93 (MQ=255)

                     cATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATc                                                                    >  1:4054559/1‑93 (MQ=255)

                     cATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATc                                                                    >  2:1769932/1‑93 (MQ=255)

                      aTTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCg                                                                   >  2:1912317/1‑93 (MQ=255)

                       tttAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCgc                                                                  >  2:1222371/1‑92 (MQ=255)

                       gttAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCgg                                                                  <  2:1674243/92‑1 (MQ=255)

                        ttAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGc                                                                 >  2:2263902/1‑93 (MQ=255)

                          aGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAt                                                               <  2:682978/93‑1 (MQ=255)

                          aGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAt                                                               <  2:4058253/93‑1 (MQ=255)

                          aGCCAATTAAGTGCGGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAt                                                               <  2:533509/93‑1 (MQ=255)

                           gCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAtt                                                              <  2:547696/93‑1 (MQ=255)

                           gCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAtt                                                              <  2:4704158/93‑1 (MQ=255)

                           gCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAtt                                                              <  2:3020838/93‑1 (MQ=255)

                             cAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTc                                                            >  1:802159/1‑93 (MQ=255)

                              aaTTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTcc                                                           <  2:3630298/93‑1 (MQ=255)

                                 tAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCcgcg                                                        >  2:1549213/1‑93 (MQ=255)

                                     tGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt                                                    <  2:3768048/93‑1 (MQ=255)

                                       cTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAg                                                  <  1:1077002/93‑1 (MQ=255)

                                       cTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAg                                                  >  2:2903264/1‑93 (MQ=255)

                                        tGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGa                                                 <  1:4290408/93‑1 (MQ=255)

                                         gAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGAc                                                >  1:1842317/1‑93 (MQ=255)

                                         gAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGAc                                                >  1:3150528/1‑93 (MQ=255)

                                             gCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAAc                                            >  2:4394718/1‑93 (MQ=255)

                                              cAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                           >  2:4586975/1‑93 (MQ=255)

                                              cAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                           <  1:3751964/93‑1 (MQ=255)

                                              cAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                           <  1:1830245/93‑1 (MQ=255)

                                               aGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCCCCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTg                                          <  1:2144649/93‑1 (MQ=255)

                                               aGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTg                                          <  2:3248370/93‑1 (MQ=255)

                                               aGCGTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTg                                          <  1:352327/93‑1 (MQ=255)

                                                gCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTgg                                         >  1:42585/1‑93 (MQ=255)

                                                 cTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTggg                                        <  2:2528251/93‑1 (MQ=255)

                                                 cTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTggg                                        >  1:4627665/1‑93 (MQ=255)

                                                 cTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTggg                                        <  1:892471/93‑1 (MQ=255)

                                                   tCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTc                                      <  1:3132353/93‑1 (MQ=255)

                                                   tCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTc                                      >  2:3855361/1‑93 (MQ=255)

                                                   tCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTc                                      >  2:3250449/1‑93 (MQ=255)

                                                      ccAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGa                                   <  1:1001861/93‑1 (MQ=255)

                                                      ccAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGGCTGa                                   >  1:2343186/1‑93 (MQ=255)

                                                        aGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGAcc                                 >  2:258789/1‑93 (MQ=255)

                                                        aGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGAcc                                 <  2:3973283/93‑1 (MQ=255)

                                                        aGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGAcc                                 >  2:2651383/1‑93 (MQ=255)

                                                         gCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCa                                >  1:4392995/1‑93 (MQ=255)

                                                         gCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTAGGTCTGACCa                                <  2:2269021/93‑1 (MQ=255)

                                                          ccTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCAt                               <  2:3556558/93‑1 (MQ=255)

                                                          ccTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCAt                               >  2:837784/1‑93 (MQ=255)

                                                          ccTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCAt                               <  1:3888047/93‑1 (MQ=255)

                                                           cTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATg                              <  2:1446120/93‑1 (MQ=255)

                                                           cTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATg                              <  2:3902718/93‑1 (MQ=255)

                                                            tGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATgg                             >  1:2769691/1‑93 (MQ=255)

                                                             gAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGt                            <  1:3506797/93‑1 (MQ=255)

                                                              aGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGtt                           <  1:1451434/93‑1 (MQ=255)

                                                                   gCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgatga                      >  2:631837/1‑93 (MQ=255)

                                                                     gTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgaga                    <  2:1382233/93‑1 (MQ=255)

                                                                     gTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgaga                    <  2:498609/93‑1 (MQ=255)

                                                                       gAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGc                  <  2:1618658/93‑1 (MQ=255)

                                                                       gAAGTGGCGCGTGAAGCGATCATTGCTCCCTCCTCTTTTTATCGGCCTTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGc                  <  1:686217/93‑1 (MQ=255)

                                                                        aaGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCg                 >  2:2635233/1‑93 (MQ=255)

                                                                         aGTGGCGCGTGAAGCGGTCATTGCTCCCATCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCgg                <  1:1027937/93‑1 (MQ=255)

                                                                         aGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCgg                <  1:1403493/93‑1 (MQ=255)

                                                                           tGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGtt              >  2:916537/1‑93 (MQ=255)

                                                                           tGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGtt              <  2:3317260/93‑1 (MQ=255)

                                                                            ggCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGttt             <  2:3518837/93‑1 (MQ=255)

                                                                              cgcgTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTaa           <  2:4230801/93‑1 (MQ=255)

                                                                              cgcgTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTaa           <  1:3565031/93‑1 (MQ=255)

                                                                                cgTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATg         <  2:3047310/93‑1 (MQ=255)

                                                                                cgTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATg         <  2:4266108/93‑1 (MQ=255)

                                                                                 gTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGc        <  1:1748697/93‑1 (MQ=255)

                                                                                   gAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTa      <  1:975721/93‑1 (MQ=255)

                                                                                   gAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTa      <  1:3647894/93‑1 (MQ=255)

                                                                                    aaGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTAc     <  1:4707542/93‑1 (MQ=255)

                                                                                    aaGCGGTCATTGCTCCCACCTCTTTTGATCGGCATTTCCGCGACGTAGACGAACAGGGTCTGACCAAGGTTGATGAGAGCGGTGGAATGCGAc     >  2:4219817/1‑93 (MQ=255)

                                                                                     aGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACg    >  1:650506/1‑93 (MQ=255)

                                                                                     aGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACg    >  2:4143958/1‑93 (MQ=255)

                                                                                      gCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGc   <  1:1726231/93‑1 (MQ=255)

                                                                                      gCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGc   <  2:4011650/93‑1 (MQ=255)

                                                                                      gCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGc   <  2:3485467/93‑1 (MQ=255)

                                                                                       cGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGcc  <  1:558160/93‑1 (MQ=255)

                                                                                       cGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGcc  <  1:3024773/93‑1 (MQ=255)

GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCC  >  NC_000913/4158737‑4158916

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