BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,067,441	(A)_6→5	coding (138/174 nt)	b4516 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,067,436	0	A	.	99.0%	930.3 / 0.9	206	coding (133/174 nt)	b4516	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base . (91/113); total (93/113)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 226 total aligned reads displayed.

ACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAA  >  NC_000913/1067344‑1067524

 ccGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTaa                                                                                         >  2:692038/1‑93 (MQ=255)

   gcgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTaaaa                                                                                       >  1:3532479/1‑93 (MQ=255)

      aaaGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTaaaaagg                                                                                    <  2:3000804/93‑3 (MQ=255)

       aaGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAggg                                                                                   >  2:4874039/1‑93 (MQ=255)

        agagCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTaaaaaggtg                                                                                  >  1:4426160/1‑89 (MQ=255)

          agCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAggtggt                                                                               <  1:2997108/93‑1 (MQ=255)

            cATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTaa                                                                             <  1:7518283/93‑1 (MQ=255)

            cATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTaa                                                                             <  1:3426548/93‑1 (MQ=255)

             aTCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTTGTAag                                                                            <  2:5376053/93‑1 (MQ=255)

             aTCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAag                                                                            >  2:974778/1‑93 (MQ=255)

              tCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAaga                                                                           <  1:5688394/93‑1 (MQ=255)

              tCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAaga                                                                           >  2:7524814/1‑93 (MQ=255)

                aGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGc                                                                         <  2:4289254/93‑1 (MQ=255)

                 gAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCa                                                                        >  2:1012516/1‑93 (MQ=255)

                 gAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGGGGTAAGAGCa                                                                        >  1:2251972/1‑93 (MQ=255)

                  aaGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAg                                                                       >  2:169362/1‑93 (MQ=255)

                   aGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGt                                                                      >  1:5932624/1‑93 (MQ=255)

                   aGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGt                                                                      <  2:5280096/93‑1 (MQ=255)

                   aGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGt                                                                      >  2:4001397/1‑93 (MQ=255)

                     cAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCa                                                                    <  1:5468033/93‑1 (MQ=255)

                     cAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCa                                                                    <  2:4461786/93‑1 (MQ=255)

                      aGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCAc                                                                   <  2:7858140/93‑1 (MQ=255)

                      aGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCAc                                                                   <  1:6097916/93‑1 (MQ=255)

                        gTAAAAAAGGTGGACAGCAGAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACgg                                                                 <  1:5267177/93‑1 (MQ=255)

                        gTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACgg                                                                 >  2:5677071/1‑93 (MQ=255)

                          aaaaaaGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCa                                                               >  2:3421912/1‑93 (MQ=255)

                           aaaaaGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCaa                                                              <  1:381909/93‑1 (MQ=255)

                           aaaaaGGTGGACAGAACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCaa                                                              >  2:2784909/1‑93 (MQ=255)

                            aaaaGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCaaa                                                             >  1:7804037/1‑93 (MQ=255)

                            aaaaGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCaaa                                                             <  1:7638169/93‑1 (MQ=255)

                             aaaGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCaaaa                                                            <  1:113687/93‑1 (MQ=255)

                             aaaGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCaaaa                                                            <  2:6155129/93‑1 (MQ=255)

                               aGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGc                                                          >  1:2341125/1‑93 (MQ=255)

                                ggTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCg                                                         <  2:3662281/93‑1 (MQ=255)

                                  tGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGAc                                                       >  2:7231572/1‑93 (MQ=255)

                                    gACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACaa                                                     <  2:5145294/93‑1 (MQ=255)

                                     acagcacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAAc                                                    <  2:7117944/93‑1 (MQ=255)

                                      cagcacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACt                                                   <  2:7661883/93‑1 (MQ=255)

                                      cagcacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACt                                                   <  1:1000712/93‑1 (MQ=255)

                                         cacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGc                                                >  2:3449632/1‑93 (MQ=255)

                                         cacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGc                                                <  1:4831341/93‑1 (MQ=255)

                                           cagcGGGGGTATTTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCg                                              >  2:3279184/1‑93 (MQ=255)

                                           cagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCg                                              <  1:306965/93‑1 (MQ=255)

                                            agcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCgg                                             <  1:6205483/93‑1 (MQ=255)

                                            agcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCgg                                             >  2:4108018/1‑93 (MQ=255)

                                             gcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCggg                                            <  1:451666/93‑1 (MQ=255)

                                             gcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCggg                                            <  1:4493581/93‑1 (MQ=255)

                                              cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGc                                           <  1:7146818/93‑1 (MQ=255)

                                               gggggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCt                                          <  1:2591823/93‑1 (MQ=255)

                                               gggggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCt                                          >  1:5223363/1‑93 (MQ=255)

                                               gggggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCt                                          <  2:7397370/93‑1 (MQ=255)

                                                ggggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTa                                         >  2:5846226/1‑93 (MQ=255)

                                                 gggTAATTTCAAGAATGACCCGCAGCGCGCATCTTAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTaa                                        <  1:2053088/93‑1 (MQ=255)

                                                  ggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAAt                                       <  1:5736862/93‑1 (MQ=255)

                                                  ggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAAt                                       <  1:4946208/93‑1 (MQ=255)

                                                  ggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAAt                                       >  1:4486865/1‑93 (MQ=255)

                                                  ggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAAt                                       >  1:3949184/1‑93 (MQ=255)

                                                   gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATc                                      <  1:7866242/93‑1 (MQ=255)

                                                    tAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCa                                     >  1:2625885/1‑93 (MQ=255)

                                                     aaTTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCaa                                    <  1:562869/93‑1 (MQ=255)

                                                        ttCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGa                                 >  2:6853488/1‑93 (MQ=255)

                                                        ttCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGa                                 <  1:1247749/93‑1 (MQ=255)

                                                         tCAAGAATGGCCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGAc                                <  1:6001613/93‑1 (MQ=255)

                                                           aaGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGa                              >  1:6039039/1‑93 (MQ=255)

                                                           aaGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGa                              >  2:5106728/1‑93 (MQ=255)

                                                           aaGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGa                              <  1:7177756/93‑1 (MQ=255)

                                                            aGAATGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCATGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGaa                             <  2:2937528/93‑1 (MQ=255)

                                                              aaTGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATg                           <  1:3206924/93‑1 (MQ=255)

                                                              aaTGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATg                           <  2:1899977/93‑1 (MQ=255)

                                                               aTGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGc                          <  1:2818418/93‑1 (MQ=255)

                                                                tGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCa                         >  2:4565012/1‑93 (MQ=255)

                                                                tGACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCa                         >  2:1312276/1‑93 (MQ=255)

                                                                 gACCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAt                        >  2:7015402/1‑93 (MQ=255)

                                                                  aCCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAtt                       <  2:5762541/93‑1 (MQ=255)

                                                                  aCCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAtt                       <  2:2108662/93‑1 (MQ=255)

                                                                  aCCCGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAtt                       <  2:4552307/93‑1 (MQ=255)

                                                                   cccGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTTGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttt                      <  1:1638741/93‑1 (MQ=255)

                                                                   cccGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttt                      >  1:6827401/1‑93 (MQ=255)

                                                                   cccGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttt                      >  2:5232793/1‑93 (MQ=255)

                                                                    ccGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAtttt                     <  2:2278700/93‑1 (MQ=255)

                                                                     cGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttttt                    <  1:4974263/93‑1 (MQ=255)

                                                                     cGCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttttt                    >  1:3567105/1‑93 (MQ=255)

                                                                      gCAGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTg                   >  1:1670188/1‑93 (MQ=255)

                                                                      gCAGAGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTg                   >  2:4756548/1‑93 (MQ=255)

                                                                        aGCGCGCATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTc                 >  1:7143943/1‑93 (MQ=255)

                                                                           gcgcATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGt              >  1:553331/1‑93 (MQ=255)

                                                                           gcgcATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGt              <  1:1034948/93‑1 (MQ=255)

                                                                            cgcATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTa             <  1:6569069/93‑1 (MQ=255)

                                                                            cgcATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTa             <  1:5117389/93‑1 (MQ=255)

                                                                             gcATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAg            <  1:6299421/93‑1 (MQ=255)

                                                                              cATCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGc           >  1:4253824/1‑93 (MQ=255)

                                                                               aTCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCt          <  2:641223/93‑1 (MQ=255)

                                                                               aTCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCt          >  2:5992256/1‑93 (MQ=255)

                                                                               aTCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCt          >  2:236080/1‑93 (MQ=255)

                                                                                tCTGAAGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTc         <  1:3955497/93‑1 (MQ=255)

                                                                                  tGAAGCAGGT‑AAAAAGGTGGGAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGt       >  2:2463785/1‑93 (MQ=255)

                                                                                    aaGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCa     >  1:1912999/1‑93 (MQ=255)

                                                                                    aaGCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCa     >  2:2544470/1‑93 (MQ=255)

                                                                                      gCAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCaaa   >  2:753937/1‑93 (MQ=255)

                                                                                       cAGGT‑AAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCaaaa  <  1:6554004/93‑1 (MQ=255)

ACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAA  >  NC_000913/1067344‑1067524

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