BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,994,771	A→G	H291R (CAT→CGT)	xerC →	XerC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,994,771	0	A	G	99.6%	845.0 / ‑4.8	240	H291R (CAT→CGT)	xerC	XerC
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (117/122); total (117/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

Only 100 of 246 total aligned reads displayed.

AGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTA  >  NC_000913/3994681‑3994860

aGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTc                                                                                         >  2:2845534/1‑93 (MQ=255)

  gAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCa                                                                                       <  1:6296290/93‑1 (MQ=255)

    gctgctGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCAcg                                                                                     >  1:4757603/1‑93 (MQ=255)

        ctgGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGcc                                                                                 <  2:4268718/93‑1 (MQ=255)

         tgGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCa                                                                                >  1:2619128/1‑93 (MQ=255)

          gggTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCaa                                                                               >  2:3720057/1‑93 (MQ=255)

            gTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAAc                                                                             >  1:1381454/1‑93 (MQ=255)

            gTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAAc                                                                             <  1:730968/93‑1 (MQ=255)

              cATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACgg                                                                           >  1:2836138/1‑93 (MQ=255)

               aTGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACggg                                                                          >  2:2712990/1‑93 (MQ=255)

                 gCCACCCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACggggg                                                                        <  1:2677656/93‑1 (MQ=255)

                  ccAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGa                                                                       >  2:4193885/1‑93 (MQ=255)

                    aaCCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGaaa                                                                     >  1:3342436/1‑93 (MQ=255)

                     aCCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAaat                                                                    <  1:2197852/93‑1 (MQ=255)

                      ccTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAaata                                                                   <  1:5427181/93‑1 (MQ=255)

                      ccTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAaata                                                                   >  2:5031950/1‑93 (MQ=255)

                         ctcCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATg                                                                >  1:6204098/1‑93 (MQ=255)

                         ctcCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATg                                                                >  2:5228564/1‑93 (MQ=255)

                         ctcCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATg                                                                <  2:3182424/93‑1 (MQ=255)

                            caccacGCAAATCTATACTCATCTTGATTTTCAACGCCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGt                                                             <  2:1864437/93‑1 (MQ=255)

                            caccacGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGt                                                             <  1:6080084/93‑1 (MQ=255)

                            caccacGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGt                                                             <  1:1001273/93‑1 (MQ=255)

                             accacGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGtt                                                            <  1:914534/93‑1 (MQ=255)

                             accacGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGtt                                                            <  1:3639142/93‑1 (MQ=255)

                              ccacGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGttt                                                           <  1:1278652/93‑1 (MQ=255)

                              ccacGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGttt                                                           <  1:2802555/93‑1 (MQ=255)

                              ccacGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGttt                                                           <  1:144471/93‑1 (MQ=255)

                                acGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGttttt                                                         >  2:4922533/1‑93 (MQ=255)

                                acGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGttttt                                                         >  2:4156045/1‑93 (MQ=255)

                                   cAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTAcc                                                      >  2:1136234/1‑93 (MQ=255)

                                   cAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTAcc                                                      <  2:785283/93‑1 (MQ=255)

                                    aaaTCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCg                                                     <  2:383013/93‑1 (MQ=255)

                                    aaaTCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCg                                                     >  2:886599/1‑93 (MQ=255)

                                     aaTCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCgg                                                    >  2:2225355/1‑93 (MQ=255)

                                      aTCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGc                                                   >  1:5007989/1‑93 (MQ=255)

                                       tCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGcc                                                  <  1:2965698/93‑1 (MQ=255)

                                       tCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGcc                                                  <  2:6292198/93‑1 (MQ=255)

                                        cTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCt                                                 <  2:1424979/93‑1 (MQ=255)

                                           taCTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTg                                              >  1:4284961/1‑93 (MQ=255)

                                            gcTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTgg                                             <  1:1840293/92‑1 (MQ=255)

                                            aCTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGACATAATGCGTTTTTACCGGGCTTTgg                                             >  2:2917002/1‑93 (MQ=255)

                                            aCTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTgg                                             <  2:5188877/93‑1 (MQ=255)

                                            aCTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTgg                                             <  1:339164/93‑1 (MQ=255)

                                            aCTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTgg                                             <  2:5807001/93‑1 (MQ=255)

                                              tCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTgggg                                           >  1:5595424/1‑93 (MQ=255)

                                              tCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTgggg                                           >  2:5484757/1‑93 (MQ=255)

                                               cATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGgc                                          <  1:4268494/93‑1 (MQ=255)

                                               cATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGgc                                          >  1:2382546/1‑93 (MQ=255)

                                               cATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGgc                                          <  1:5088729/93‑1 (MQ=255)

                                                 tCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGgcgc                                        <  1:167813/93‑1 (MQ=255)

                                                 tCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGgcgc                                        <  2:1641566/93‑1 (MQ=255)

                                                 tCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGgcgc                                        >  1:1977972/1‑93 (MQ=255)

                                                  cTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGACTTTGGGGCGCa                                       >  1:834430/1‑93 (MQ=255)

                                                   ttGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCAt                                      >  1:3167278/1‑93 (MQ=255)

                                                    tGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCAtc                                     <  1:3872410/93‑1 (MQ=255)

                                                     gATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCAtct                                    <  2:6110204/93‑1 (MQ=255)

                                                      aTTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCAtctc                                   <  2:3704072/93‑1 (MQ=255)

                                                           cAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGgcgc                              <  2:568596/93‑1 (MQ=255)

                                                           cAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGgcgc                              <  1:5791711/93‑1 (MQ=255)

                                                           cAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGgcgc                              <  2:5941875/93‑1 (MQ=255)

                                                           cAACACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGgcgc                              <  2:4777043/93‑1 (MQ=255)

                                                            aaCACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCt                             >  1:4233148/1‑93 (MQ=255)

                                                            aaCACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCt                             >  2:4857907/1‑93 (MQ=255)

                                                            aaCACCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCt                             >  2:1289082/1‑93 (MQ=255)

                                                             acacCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTc                            <  1:5948192/93‑1 (MQ=255)

                                                               acCTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCAc                          >  1:5270977/1‑93 (MQ=255)

                                                                ccTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCAcc                         <  1:3536456/93‑1 (MQ=255)

                                                                ccTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCAcc                         <  1:2468752/93‑1 (MQ=255)

                                                                 cTTGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCt                        >  1:3107743/1‑93 (MQ=255)

                                                                  ttGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCtt                       >  2:2972164/1‑93 (MQ=255)

                                                                  ttGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCtt                       >  2:1606999/1‑93 (MQ=255)

                                                                  ttGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCtt                       <  2:6323949/93‑1 (MQ=255)

                                                                   tGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCttt                      <  1:5784305/93‑1 (MQ=255)

                                                                   tGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCttt                      <  2:2309089/93‑1 (MQ=255)

                                                                   tGCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCttt                      >  2:4447223/1‑93 (MQ=255)

                                                                    gCCTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTg                     <  2:2830024/93‑1 (MQ=255)

                                                                     ccTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGa                    >  1:4185432/1‑93 (MQ=255)

                                                                      cTCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATACGTTTTTACCGGCCTTTTGGGCGCATCTCGGCGCTCCCCttttcc                   >  2:4617151/1‑90 (MQ=255)

                                                                       tCGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGAcc                  <  2:5399414/93‑1 (MQ=255)

                                                                        cGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCt                 >  2:2012525/1‑93 (MQ=255)

                                                                        cGGTGTACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCt                 <  1:98523/93‑1 (MQ=255)

                                                                           tgtACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGga              <  2:23363/93‑1 (MQ=255)

                                                                           tgtACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGga              <  2:5332071/93‑1 (MQ=255)

                                                                           tgtACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGga              <  1:5904216/93‑1 (MQ=255)

                                                                            gtACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGgat             <  1:2052116/93‑1 (MQ=255)

                                                                            gtACGATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGgat             <  1:1044576/93‑1 (MQ=255)

                                                                                gATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATa         <  2:2475882/93‑1 (MQ=255)

                                                                                gATGCGGCGCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATa         >  1:803858/1‑93 (MQ=255)

                                                                                   gcggcgCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGAGCTGGATGATAccc      >  2:3030245/1‑93 (MQ=255)

                                                                                    cgtcgcgTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCt     <  1:3938109/90‑1 (MQ=255)

                                                                                    cggcgCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCt     <  1:4047859/93‑1 (MQ=255)

                                                                                    cggcgCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCt     >  1:1518044/1‑93 (MQ=255)

                                                                                    cggcgCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCt     <  2:1206003/93‑1 (MQ=255)

                                                                                     ggcgCGTCCGCGCGCCAAACGGGGGACATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATGCCCtt    <  2:4628054/93‑1 (MQ=255)

                                                                                     ggcgCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCtt    >  2:3569019/1‑93 (MQ=255)

                                                                                     ggcgCGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCtt    <  2:3339047/93‑1 (MQ=255)

                                                                                      gcgcGTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCttt   >  1:4524298/1‑93 (MQ=255)

                                                                                      gcgcGTCCACGCGCCAAACGGGGGAAATAATGCGGTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGCTGATACCCttt   >  1:6093120/1‑93 (MQ=255)

                                                                                       cgcgTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTGCCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTa  <  1:1531580/93‑1 (MQ=255)

                                                                                       cgcgTCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTa  <  1:2990228/93‑1 (MQ=255)

AGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTA  >  NC_000913/3994681‑3994860

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