BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,679,306	T→C	S89G (AGC→GGC)	yhjJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,679,306	0	T	C	99.6%	727.7 / ‑5.0	236	S89G (AGC→GGC)	yhjJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (109/126); total (109/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

Only 100 of 247 total aligned reads displayed.

CATAAGAGACAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAA  >  NC_000913/3679216‑3679395

aaTAAGAGACAATAACAGGCGGCATCGGGCGATTAGAGCAAAACCCCTGATGCCACAATGAACGCGCCTGTGCAGCGTCAAGGCCACCGCctt                                                                                         <  2:4681126/92‑1 (MQ=255)

   aaGAGACAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCg                                                                                      >  1:2850831/1‑93 (MQ=255)

   aaGAGACAATTACCGGCGGCATCGGGCGTTTAGGGTCGAGCTCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCg                                                                                      >  2:19260/1‑93 (MQ=255)

    agagACAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGt                                                                                     >  2:1212535/1‑93 (MQ=255)

     gagaCAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTACTGCGTCAAGGACACCGCCTTGCGtt                                                                                    >  2:501931/1‑93 (MQ=255)

      agaCAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTa                                                                                   <  2:4046781/93‑1 (MQ=255)

      agaCAATTACCAGCGGCATCGGTCTTTTAGGATCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTa                                                                                   <  1:1952001/93‑1 (MQ=255)

        aCAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAgc                                                                                 <  2:1358901/93‑1 (MQ=255)

        aCAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAgc                                                                                 <  1:2495531/93‑1 (MQ=255)

           aTTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCa                                                                              >  1:1782979/1‑93 (MQ=255)

           aTTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCa                                                                              >  1:3514887/1‑93 (MQ=255)

            ttACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCaa                                                                             >  2:1707203/1‑93 (MQ=255)

            ttACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCaa                                                                             <  1:2082597/93‑1 (MQ=255)

            ttACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCaa                                                                             >  2:3813939/1‑93 (MQ=255)

             tACCGGCGGCATGGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAAt                                                                            <  2:3938934/93‑1 (MQ=255)

               ccGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATAc                                                                          <  2:1462557/93‑1 (MQ=255)

               ccGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATAc                                                                          >  2:2187551/1‑93 (MQ=255)

                cggcggCATCGGGCTTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACg                                                                         <  1:2463267/93‑1 (MQ=255)

                 ggcggcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGa                                                                        <  1:1196016/93‑1 (MQ=255)

                 ggcggcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGa                                                                        >  2:2463099/1‑93 (MQ=255)

                   cggcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAgg                                                                      <  1:565725/93‑1 (MQ=255)

                    ggcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAggg                                                                     >  2:2663840/1‑93 (MQ=255)

                    ggcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAggg                                                                     <  2:2355902/93‑1 (MQ=255)

                    ggcATCGGGCGTTTAGGGTCGATCCCCTGCTG‑CACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGa                                                                    <  1:3555716/93‑1 (MQ=255)

                     gcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGGGGGa                                                                    >  2:691318/1‑93 (MQ=255)

                     gcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGa                                                                    <  2:4652898/93‑1 (MQ=255)

                     gcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGa                                                                    <  1:168250/93‑1 (MQ=255)

                     gcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGa                                                                    <  1:3862329/93‑1 (MQ=255)

                     gcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGa                                                                    <  2:3517829/93‑1 (MQ=255)

                     gcATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGa                                                                    <  2:4207883/93‑1 (MQ=255)

                       aTCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATg                                                                  >  1:2428604/1‑93 (MQ=255)

                       aTCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATg                                                                  >  2:4425520/1‑93 (MQ=255)

                          gggCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCg                                                               >  2:3851185/1‑93 (MQ=255)

                            gCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTg                                                             <  1:753150/93‑1 (MQ=255)

                             cGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGa                                                            >  1:4430225/1‑93 (MQ=255)

                             cGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGa                                                            <  2:2565390/93‑1 (MQ=255)

                              gTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGAc                                                           <  2:1437606/93‑1 (MQ=255)

                              gTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGAc                                                           <  2:3598055/93‑1 (MQ=255)

                               tttAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACt                                                          <  1:4614393/93‑1 (MQ=255)

                               tttAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACt                                                          >  2:222907/1‑93 (MQ=255)

                               tttAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACt                                                          >  2:3230658/1‑93 (MQ=255)

                               tttAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACt                                                          >  2:3287595/1‑93 (MQ=255)

                                ttAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTg                                                         <  2:993605/93‑1 (MQ=255)

                                 tAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGt                                                        >  2:1139224/1‑93 (MQ=255)

                                   gggTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTaa                                                      <  2:362019/93‑1 (MQ=255)

                                    ggTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAAc                                                     <  1:1045430/93‑1 (MQ=255)

                                     gTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAAcc                                                    >  1:3200053/1‑93 (MQ=255)

                                     gTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAAcc                                                    >  2:710799/1‑93 (MQ=255)

                                      tCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGGGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGTCATGACTGTAACCg                                                   >  2:3125127/1‑93 (MQ=255)

                                      tCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCg                                                   >  1:1855898/1‑93 (MQ=255)

                                       cGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGc                                                  <  2:3256047/93‑1 (MQ=255)

                                       cGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCAACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGc                                                  >  1:2632553/1‑93 (MQ=255)

                                         aTCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctc                                                >  2:416551/1‑93 (MQ=255)

                                         aTCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctc                                                <  1:3971368/93‑1 (MQ=255)

                                          tCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGGGTGCCTGTAACCGctct                                               >  2:1054498/1‑93 (MQ=255)

                                          tCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctct                                               >  1:2369253/1‑93 (MQ=255)

                                          tCCCCGGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGctct                                               >  2:2968884/1‑93 (MQ=255)

                                            cccTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCtgt                                             >  1:2011610/1‑93 (MQ=255)

                                             ccTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCtgtt                                            >  1:1452686/1‑93 (MQ=255)

                                              ctgctgCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCtgttg                                           >  1:1174085/1‑93 (MQ=255)

                                               tgctgcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCtgttgt                                          <  1:674138/93‑1 (MQ=255)

                                                gctgcCCCAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTTCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTtgtg                                         <  1:2557609/93‑1 (MQ=255)

                                                gctgcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTtgtg                                         >  1:325642/1‑93 (MQ=255)

                                                gctgcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTtgtg                                         <  2:3327683/93‑1 (MQ=255)

                                                 ctgcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTgtgt                                        <  1:928793/93‑1 (MQ=255)

                                                 ctgcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTgtgt                                        >  1:1878202/1‑93 (MQ=255)

                                                   gcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTat                                      >  2:3970460/1‑92 (MQ=255)

                                                   gcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTAc                                      <  2:2396895/93‑1 (MQ=255)

                                                   gcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTAc                                      <  1:4382815/93‑1 (MQ=255)

                                                   gcCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTAc                                      <  1:1266098/93‑1 (MQ=255)

                                                    ccACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACt                                     <  1:2151391/93‑1 (MQ=255)

                                                     cacaATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACtt                                    <  2:95300/93‑1 (MQ=255)

                                                         aTGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCCATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCgg                                >  2:4469242/1‑93 (MQ=255)

                                                         aTGAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCgg                                >  1:3713509/1‑93 (MQ=255)

                                                           gAACGCGCCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGgcg                              <  1:228838/93‑1 (MQ=255)

                                                               gcgcCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGgcgagcg                          <  1:631160/93‑1 (MQ=255)

                                                                cgcCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGgcgagcga                         >  2:911607/1‑93 (MQ=255)

                                                                 gcCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGaa                        >  1:1076617/1‑93 (MQ=255)

                                                                 gcCTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGaa                        >  1:4224318/1‑93 (MQ=255)

                                                                   cTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAAcc                      <  1:4124769/93‑1 (MQ=255)

                                                                   cTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAAcc                      >  1:3414303/1‑93 (MQ=255)

                                                                   cTGTGCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAAcc                      >  1:3986097/1‑93 (MQ=255)

                                                                    tgtgCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCg                     <  1:4513361/93‑1 (MQ=255)

                                                                    tgtgCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCg                     >  1:4275815/1‑93 (MQ=255)

                                                                    tgtgCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCg                     <  1:3913184/93‑1 (MQ=255)

                                                                    tgtgCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCg                     >  2:3034099/1‑93 (MQ=255)

                                                                    tgtgCTGCGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCg                     <  2:3576794/93‑1 (MQ=255)

                                                                      tgctgcGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGt                   <  1:4141141/93‑1 (MQ=255)

                                                                      tgctgcGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGt                   <  1:2985475/93‑1 (MQ=255)

                                                                       gctgcGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTa                  >  2:3607539/1‑93 (MQ=255)

                                                                           cGTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTaa              <  2:539327/93‑1 (MQ=255)

                                                                            gTCAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAAc             <  2:1310698/93‑1 (MQ=255)

                                                                              cAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTACCCTCTCTGTTGTGTACTTTCGGCGAGCGAGCCGGTATTAACca           <  2:1863221/93‑1 (MQ=255)

                                                                              cAAGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACca           <  1:1508592/93‑1 (MQ=255)

                                                                                aGGCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACcagc         >  2:4348202/1‑93 (MQ=255)

                                                                                 ggCCACCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACcagca        <  2:804242/93‑1 (MQ=255)

                                                                                     aCCGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGAAGGCg    >  2:1657473/1‑93 (MQ=255)

                                                                                       cGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGaa  <  1:3796377/93‑1 (MQ=255)

                                                                                       cGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGaa  >  2:4262568/1‑93 (MQ=255)

                                                                                       cGCCTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGaa  <  1:780577/93‑1 (MQ=255)

CATAAGAGACAATTACCGGCGGCATCGGGCGTTTAGGGTCGATCCCCTGCTGCCACAATGAACGCGCCTGTGCTGCGTCAAGGCCACCGCTTTGCGTTAGCGCAATACGAGGGATGGCGTGACTGTAACCGCTCTGTTGTGTACTTTCGGCGAGCGAACCGGTATTAACCAGCAGGCGAA  >  NC_000913/3679216‑3679395

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