BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	609,564	T→C	A718A (GCA→GCG)	fepA ←	FepA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	609,564	0	T	C	82.9%	559.3 / 112.4	240	A718A (GCA→GCG)	fepA	FepA
Reads supporting (aligned to +/- strand): ref base T (23/18); new base C (102/97); total (125/115)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

Only 100 of 247 total aligned reads displayed.

AGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGACTGACATTCTT  >  NC_000913/609473‑609653

 gCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGc                                                                                         <  2:4299807/93‑1 (MQ=255)

 gCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGc                                                                                         >  1:207043/1‑93 (MQ=255)

 gCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGc                                                                                         >  2:6096394/1‑93 (MQ=255)

 gCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGc                                                                                         >  2:6040580/1‑93 (MQ=255)

   aTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCa                                                                                       >  2:3965419/1‑93 (MQ=255)

   aTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCa                                                                                       <  2:5504458/93‑1 (MQ=255)

   aTCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCa                                                                                       <  1:5760326/93‑1 (MQ=255)

    tCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCaa                                                                                      >  2:1170094/1‑93 (MQ=255)

    tCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCaa                                                                                      <  1:4733637/93‑1 (MQ=255)

    tCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCaa                                                                                      >  1:4361863/1‑93 (MQ=255)

     cAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCaaa                                                                                     <  1:5345106/93‑1 (MQ=255)

       gAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATc                                                                                   <  2:2713833/93‑1 (MQ=255)

       gAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATc                                                                                   <  2:1618925/93‑1 (MQ=255)

       gAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATc                                                                                   >  2:740419/1‑93 (MQ=255)

       gAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATc                                                                                   >  2:6356715/1‑93 (MQ=255)

        aaGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCg                                                                                  >  1:1563515/1‑93 (MQ=255)

          gTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGcc                                                                                >  1:2854846/1‑93 (MQ=255)

          gTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGcc                                                                                <  2:2141255/93‑1 (MQ=255)

           tGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGccc                                                                               >  2:6308755/1‑93 (MQ=255)

           tGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGccc                                                                               >  1:2363164/1‑93 (MQ=255)

            gggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCg                                                                              <  1:4141743/93‑1 (MQ=255)

             ggTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGt                                                                             <  2:5838421/93‑1 (MQ=255)

              gtgtTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTg                                                                            >  1:5400799/1‑93 (MQ=255)

                gtTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTggc                                                                          >  1:1715884/1‑92 (MQ=255)

                gtTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGt                                                                          <  1:3371690/93‑1 (MQ=255)

                 tttACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTc                                                                         >  1:2220774/1‑93 (MQ=255)

                  ttACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCt                                                                        <  1:24206/93‑1 (MQ=255)

                   tACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTg                                                                       >  2:2278938/1‑93 (MQ=255)

                   tACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTg                                                                       <  1:3005940/93‑1 (MQ=255)

                       cTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCa                                                                   >  2:4959991/1‑93 (MQ=255)

                       cTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCa                                                                   >  2:568795/1‑93 (MQ=255)

                          atatACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTa                                                                >  1:3402662/1‑93 (MQ=255)

                          atatACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTa                                                                >  1:6512000/1‑93 (MQ=255)

                          atatACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTa                                                                <  1:6052386/93‑1 (MQ=255)

                           tataCCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTAc                                                               <  2:1877552/93‑1 (MQ=255)

                             taCCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTAccc                                                             <  2:6190857/93‑1 (MQ=255)

                                cACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGcc                                                          >  1:294225/1‑92 (MQ=255)

                                cACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCa                                                          >  2:2629246/1‑93 (MQ=255)

                                cACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCa                                                          <  2:4039972/93‑1 (MQ=255)

                                 ccgtacgtCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCAc                                                         <  2:4417410/92‑1 (MQ=255)

                                    tacgtCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGcc                                                      >  2:2516941/1‑93 (MQ=255)

                                    tacgtCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGcc                                                      >  1:1094903/1‑93 (MQ=255)

                                     acgtCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCca                                                     >  1:512733/1‑93 (MQ=255)

                                     acgtCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCca                                                     >  2:2107147/1‑93 (MQ=255)

                                      cgtCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCcac                                                    >  2:4907843/1‑93 (MQ=255)

                                        tCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACaa                                                  >  1:940928/1‑93 (MQ=255)

                                        tCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGGCTGGGCATTACCCGCACGCCACaa                                                  >  1:5052572/1‑93 (MQ=255)

                                         cccGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACaaa                                                 >  2:3323289/1‑93 (MQ=255)

                                          ccGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAAc                                                >  2:1466566/1‑93 (MQ=255)

                                            ggCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGt                                              <  1:5869580/93‑1 (MQ=255)

                                            ggCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGt                                              <  2:3642709/93‑1 (MQ=255)

                                             gCTCGTTATAGATATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCCCCAGTGGTCTGGGCATTACCCGCACGCCACAAACGtt                                             <  1:605826/93‑1 (MQ=255)

                                               tCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTg                                           <  1:3085500/93‑1 (MQ=255)

                                               tCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTg                                           <  2:2462736/93‑1 (MQ=255)

                                                 gTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTc                                         >  2:1883023/1‑93 (MQ=255)

                                                 gTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTc                                         >  2:380029/1‑93 (MQ=255)

                                                 gTTATAGGGATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGACAAATCGCCCGTGGTTTGGGCATTACCCGCACGCCACAAACGTTTGTc                                         <  2:1776029/93‑1 (MQ=255)

                                                 gTTATAGGAAGACGCCCCGACACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTc                                         <  1:6330307/93‑1 (MQ=255)

                                                  ttATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCg                                        >  1:273019/1‑93 (MQ=255)

                                                  ttATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCg                                        >  1:2592083/1‑93 (MQ=255)

                                                  ttATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCg                                        >  2:6127387/1‑93 (MQ=255)

                                                   tataGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGa                                       <  1:1470346/93‑1 (MQ=255)

                                                     taGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAAc                                     >  1:5577410/1‑93 (MQ=255)

                                                      aGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACa                                    <  2:3205671/93‑1 (MQ=255)

                                                       ggTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAg                                   >  2:5138336/1‑93 (MQ=255)

                                                       ggTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAg                                   <  2:3773902/93‑1 (MQ=255)

                                                       ggTATACGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAg                                   <  2:2909205/93‑1 (MQ=255)

                                                         tataCGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGAt                                 <  1:4447404/93‑1 (MQ=255)

                                                          ataCGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGAtt                                >  1:4910617/1‑93 (MQ=255)

                                                          ataCGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGAtt                                <  1:5709409/93‑1 (MQ=255)

                                                           taCGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTg                               >  2:3037213/1‑93 (MQ=255)

                                                           taCGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTg                               <  1:750347/93‑1 (MQ=255)

                                                            aCGCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGt                              <  2:6299685/93‑1 (MQ=255)

                                                              gCCGCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTcc                            <  2:6561997/93‑1 (MQ=255)

                                                              gCCCCGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTCCCCGCACGCCACAAACGTTTGTCGAACAGATTGTcc                            <  1:3515313/93‑1 (MQ=255)

                                                                 ccGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCAcg                         <  2:3415742/93‑1 (MQ=255)

                                                                 ccGGCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCaag                         >  2:5460667/1‑91 (MQ=255)

                                                                    gCACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCAcgccg                      >  1:5936970/1‑93 (MQ=255)

                                                                    gCACCGGCGATATAGTTGGCCACTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCAcgccg                      >  1:5265490/1‑93 (MQ=255)

                                                                     cACCGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCAcgccgc                     <  1:1035651/93‑1 (MQ=255)

                                                                       ccGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGccgccg                   >  1:274263/1‑93 (MQ=255)

                                                                       ccGGCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGccgccg                   >  2:1452427/1‑93 (MQ=255)

                                                                        cGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCgg                  <  2:5122787/93‑1 (MQ=255)

                                                                         ggCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGt                 >  1:3708553/1‑93 (MQ=255)

                                                                          gCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTc                >  2:1407043/1‑93 (MQ=255)

                                                                          gCGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTc                <  1:4589013/93‑1 (MQ=255)

                                                                           cGATATAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCa               <  2:1079145/93‑1 (MQ=255)

                                                                             atatAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGa             <  1:4656733/93‑1 (MQ=255)

                                                                             atatAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGa             >  2:6637021/1‑93 (MQ=255)

                                                                             atatAGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGa             >  2:2747717/1‑93 (MQ=255)

                                                                              tataGTTGGCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGAc            <  1:3907320/93‑1 (MQ=255)

                                                                              tataGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTTTCGAACAGATTGTCCACGCCGCCGGTCAGAc            <  1:2178712/93‑1 (MQ=255)

                                                                              tataGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGAc            <  2:5248793/93‑1 (MQ=255)

                                                                              tataGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGAc            >  2:4005493/1‑93 (MQ=255)

                                                                              tataGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGAc            >  2:3865349/1‑93 (MQ=255)

                                                                              tataGTTGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGAc            <  1:4778323/93‑1 (MQ=255)

                                                                                   ttGGCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGACTGACa       <  2:3730685/93‑1 (MQ=255)

                                                                                     ggCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGACTGACAtt     <  1:12026/93‑1 (MQ=255)

                                                                                     ggCCCCCGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGACTGACAtt     >  2:5068494/1‑93 (MQ=255)

                                                                                        ccccGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGACTGACAttctt  >  2:4690842/1‑93 (MQ=255)

AGCATCAGAAGTGGGTGTTTACGCTCATATACCACGTACGTCCCGGCTCGTTATAGGTATACGCCCCGGCACCGGCGATATAGTTGGCCCCTGCCAAATCGCCCGTGGTCTGGGCATTACCCGCACGCCACAAACGTTTGTCGAACAGATTGTCCACGCCGCCGGTCAGACTGACATTCTT  >  NC_000913/609473‑609653

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