BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,181,174	G→A	G252S (GGC→AGC)	ygiC →	ATP‑Grasp family ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,181,174	0	G	A	99.1%	736.5 / ‑2.9	214	G252S (GGC→AGC)	ygiC	ATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (97/115); total (98/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

Only 100 of 214 total aligned reads displayed.

ACGGATTTACAGGATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTG  >  NC_000913/3181078‑3181269

 cGGATTTACAGGATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCg                                                                                              <  2:3304329/99‑1 (MQ=255)

  ggATTTACAGGATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGt                                                                                             <  2:1294910/99‑1 (MQ=255)

    aTTTACAGGATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTAc                                                                                           >  1:5321614/1‑99 (MQ=255)

     tttACAGGATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACg                                                                                          >  1:1579839/1‑99 (MQ=255)

        aCAGGATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACgctg                                                                                       >  1:2193070/1‑99 (MQ=255)

          aGGATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACgctggc                                                                                     <  2:3470909/99‑1 (MQ=255)

           ggATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACgctggct                                                                                    <  1:5559464/99‑1 (MQ=255)

           ggATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACgctggct                                                                                    <  1:1239773/99‑1 (MQ=255)

            gATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGACTTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACgctggctg                                                                                   >  1:727756/1‑99 (MQ=255)

            gATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACgctggctg                                                                                   >  1:3910232/1‑99 (MQ=255)

            gATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACgctggctg                                                                                   <  1:2646191/99‑1 (MQ=255)

            gATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACgctggctg                                                                                   >  1:4575719/1‑99 (MQ=255)

            gATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACgctggctg                                                                                   >  1:4583261/1‑99 (MQ=255)

                 ggTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAAcc                                                                              <  1:4520343/99‑1 (MQ=255)

                  gTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCg                                                                             >  1:851778/1‑99 (MQ=255)

                  gTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCg                                                                             <  1:3916491/99‑1 (MQ=255)

                  gTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCg                                                                             >  1:4480382/1‑99 (MQ=255)

                   tAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCgg                                                                            >  2:4123621/1‑99 (MQ=255)

                   tAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATACAAGCGTACGCTGGCTGGAACCgg                                                                            >  2:5039366/1‑99 (MQ=255)

                    aaTTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGc                                                                           <  2:3354894/99‑1 (MQ=255)

                      tttCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGt                                                                         <  2:4498178/99‑1 (MQ=255)

                       ttCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTg                                                                        <  2:636491/99‑1 (MQ=255)

                        tCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTgg                                                                       >  1:4830221/1‑99 (MQ=255)

                           aaCCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAag                                                                    >  2:4821612/1‑99 (MQ=255)

                           aaCCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAag                                                                    >  2:4354227/1‑99 (MQ=255)

                            aCCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAaga                                                                   <  1:2990447/99‑1 (MQ=255)

                            aCCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAaga                                                                   <  1:2707398/99‑1 (MQ=255)

                             ccTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAagag                                                                  <  1:38851/99‑1 (MQ=255)

                             ccTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAagag                                                                  >  2:1312668/1‑99 (MQ=255)

                              cTGTTCAAACTGGATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGc                                                                 <  1:1116927/99‑1 (MQ=255)

                              atcttCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGc                                                                 <  2:2449279/96‑1 (MQ=255)

                               tGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCa                                                                >  2:2258115/1‑99 (MQ=255)

                                gTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCAt                                                               >  2:4224584/1‑99 (MQ=255)

                                gTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCAt                                                               >  1:2846584/1‑99 (MQ=255)

                                 ttCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCAtt                                                              <  1:3448723/99‑1 (MQ=255)

                                 ttCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCAtt                                                              <  1:3071649/99‑1 (MQ=255)

                                    aaaCTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTAtc                                                           <  1:3721757/99‑1 (MQ=255)

                                    aaaCTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTAtc                                                           >  2:1964763/1‑99 (MQ=255)

                                    aaaCTGTATCCGAGGGGATTTATGATGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTAtc                                                           >  2:1069181/1‑99 (MQ=255)

                                     aaCTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTAtct                                                          >  1:3812996/1‑99 (MQ=255)

                                      aCTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTAtctc                                                         >  2:4634949/1‑99 (MQ=255)

                                      aCTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTAtctc                                                         <  1:3062169/99‑1 (MQ=255)

                                      aCTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTAtctc                                                         <  1:1342749/99‑1 (MQ=255)

                                      aCTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTAtctc                                                         >  1:1512044/1‑99 (MQ=255)

                                         gTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCcaa                                                      <  1:5188497/99‑1 (MQ=255)

                                         gTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCcaa                                                      <  2:542057/99‑1 (MQ=255)

                                          tATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCcaac                                                     >  1:5636766/1‑99 (MQ=255)

                                          tATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCcaac                                                     >  1:968488/1‑99 (MQ=255)

                                             ccGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAg                                                  <  2:3569501/99‑1 (MQ=255)

                                             ccGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAg                                                  <  2:5194588/99‑1 (MQ=255)

                                             ccGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAg                                                  >  1:3549573/1‑99 (MQ=255)

                                              cttGGGAATTTAGGTTGCGTGAGAGGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAgg                                                 <  2:3879979/97‑1 (MQ=255)

                                              cGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAgg                                                 >  2:4453298/1‑99 (MQ=255)

                                              cGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAgg                                                 <  2:2078510/99‑1 (MQ=255)

                                              cGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAgg                                                 <  1:2310390/99‑1 (MQ=255)

                                                tGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCa                                               >  1:1731777/1‑99 (MQ=255)

                                                tGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCa                                               >  2:748247/1‑99 (MQ=255)

                                                tGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCa                                               >  1:4546454/1‑99 (MQ=255)

                                                 gggAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCAc                                              <  1:5207798/99‑1 (MQ=255)

                                                 gggAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCAc                                              <  1:2371240/99‑1 (MQ=255)

                                                   gAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAGGAGCATTATCTCCAACAAGGCACtt                                            >  2:3979840/1‑99 (MQ=255)

                                                          tGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGc                                     <  1:5574528/99‑1 (MQ=255)

                                                            ttGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTa                                   >  2:3340425/1‑99 (MQ=255)

                                                            ttGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTa                                   >  1:5467095/1‑99 (MQ=255)

                                                                gTGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACtgt                               <  1:2827070/99‑1 (MQ=255)

                                                                 tGAGATGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACtgtg                              >  1:1991919/1‑99 (MQ=255)

                                                                  gagaTGTTCTCAACCAAGTTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTgg                             <  1:99723/99‑1 (MQ=255)

                                                                  gagaTGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTgg                             <  1:3394705/99‑1 (MQ=255)

                                                                    gaTGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGga                           <  2:3045073/99‑1 (MQ=255)

                                                                     aTGTTCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGgag                          >  1:1388997/1‑99 (MQ=255)

                                                                        ttCTCAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATg                       >  1:3843663/1‑99 (MQ=255)

                                                                          ctcAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGtt                     >  2:1441518/1‑99 (MQ=255)

                                                                           tcAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCGTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTc                    >  1:816911/1‑99 (MQ=255)

                                                                           tcAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTc                    <  1:2251276/99‑1 (MQ=255)

                                                                            cAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTcc                   >  2:5222908/1‑99 (MQ=255)

                                                                            cAACCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTcc                   <  1:4187907/99‑1 (MQ=255)

                                                                             aaCCAAGCTGGAGGATGCAAGCGTACGGTGGCTGGACACGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTccc                  <  1:2449047/99‑1 (MQ=255)

                                                                             aaCCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTccc                  >  1:1021845/1‑99 (MQ=255)

                                                                             aaCCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTccc                  <  2:2378027/99‑1 (MQ=255)

                                                                             aaCCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTccc                  <  2:3655962/99‑1 (MQ=255)

                                                                              aCCAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCg                 <  1:4327379/99‑1 (MQ=255)

                                                                               ccAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGa                >  1:3536984/1‑99 (MQ=255)

                                                                               ccAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGa                >  2:4800705/1‑99 (MQ=255)

                                                                               ccAAGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGa                <  1:553276/99‑1 (MQ=255)

                                                                                cAAGGTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGGGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGaa               <  2:3288955/99‑1 (MQ=255)

                                                                                 aaGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGAGAGGAAGAGCAGTATATCCAACAAGGCTCTTCTACCGCTACTGTGGGAAATGTTCCCGAAt              >  2:783461/1‑99 (MQ=255)

                                                                                  aGCTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATc             >  1:2573732/1‑99 (MQ=255)

                                                                                    cTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAGTCAc           <  2:4849612/99‑1 (MQ=255)

                                                                                    cTGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCAc           >  1:524424/1‑99 (MQ=255)

                                                                                     tGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGCGCATTATCTTCAACATGGCCCTTCTACCGCTACGGGGGGAGATGTTCGCTAATCAcc          <  2:1910241/99‑1 (MQ=255)

                                                                                     tGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCAcc          <  1:4105882/99‑1 (MQ=255)

                                                                                     tGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCAcc          <  1:3316929/99‑1 (MQ=255)

                                                                                     tGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCAcc          <  1:3454848/99‑1 (MQ=255)

                                                                                     tGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCAcc          <  1:4729379/99‑1 (MQ=255)

                                                                                     tGGAGGATGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCAcc          <  2:1593745/99‑1 (MQ=255)

                                                                                       gaggaTGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCAccag        >  2:3169852/1‑97 (MQ=255)

                                                                                       gaggaTGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCg        >  1:3575792/1‑99 (MQ=255)

                                                                                       gaggaTGCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCg        >  1:4926705/1‑99 (MQ=255)

                                                                                            tGCAAGCGTACGCTGGCTGGAACCGGCATGGAAGAGCATTATGTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACct   >  2:1176105/1‑99 (MQ=255)

                                                                                             gCAAGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACctg  <  2:4466583/99‑1 (MQ=255)

ACGGATTTACAGGATCAGGTAATTTCCAACCTGTTCAAACTGTATCCGTGGGAATTTATGTTGCGTGAGATGTTCTCAACCAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTG  >  NC_000913/3181078‑3181269

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