BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,086,713	A→G	R119R (CGT→CGC)	ycdQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,086,713	0	A	G	99.7%	1081.3 / ‑5.5	337	R119R (CGT→CGC)	ycdQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (152/183); total (153/184)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 341 total aligned reads displayed.

GGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCACGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGTTTCCTCAA  >  NC_000913/1086623‑1086801

  ttttAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCggg                                                                                      >  1:8437031/1‑93 (MQ=255)

   tttAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGt                                                                                     <  1:4521395/93‑1 (MQ=255)

     tAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGttt                                                                                   <  2:8138326/93‑1 (MQ=255)

      aaGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGtttt                                                                                  <  2:8552563/93‑1 (MQ=255)

      aaGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGtttt                                                                                  <  2:6989781/93‑1 (MQ=255)

      aaGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGAATGGCGCGGGtttt                                                                                  <  1:3160379/93‑1 (MQ=255)

        gcgcAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTAt                                                                                >  2:9808118/1‑93 (MQ=255)

        acgcAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTAt                                                                                <  1:7977069/92‑1 (MQ=255)

         cgcAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATc                                                                               <  2:8996378/93‑1 (MQ=255)

          gcAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCt                                                                              >  2:6109608/1‑93 (MQ=255)

            aaTGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGt                                                                            <  2:818581/93‑1 (MQ=255)

             aTGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCGGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGtt                                                                           >  1:1478871/1‑93 (MQ=255)

             aTGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGtt                                                                           >  2:6227302/1‑93 (MQ=255)

              tGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTg                                                                          >  1:1129317/1‑93 (MQ=255)

              tGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTg                                                                          <  2:1882587/93‑1 (MQ=255)

              tGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTg                                                                          <  1:1733929/93‑1 (MQ=255)

              tGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGGTTTATCTGTTg                                                                          >  2:7248524/1‑93 (MQ=255)

               ggCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGa                                                                         >  1:6561474/1‑93 (MQ=255)

                 cTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAAc                                                                       <  2:2904783/93‑1 (MQ=255)

                  tttCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAAcc                                                                      >  2:5159435/1‑93 (MQ=255)

                   ttCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCg                                                                     >  2:2096055/1‑93 (MQ=255)

                   ttCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCg                                                                     >  1:8655561/1‑93 (MQ=255)

                     ccccTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTc                                                                   >  1:5861248/1‑93 (MQ=255)

                     ccccTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTc                                                                   <  2:9619992/93‑1 (MQ=255)

                     ccccTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGCAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTc                                                                   <  1:4670574/93‑1 (MQ=255)

                      cccTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCa                                                                  >  1:4437691/1‑93 (MQ=255)

                       ccTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCAt                                                                 <  1:8038585/93‑1 (MQ=255)

                       ccTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCAt                                                                 <  1:3361276/93‑1 (MQ=255)

                       ccTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCAt                                                                 <  1:2351336/93‑1 (MQ=255)

                        cgggTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCAtt                                                                <  1:4251490/91‑1 (MQ=255)

                        cTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCAtt                                                                >  2:2327093/1‑93 (MQ=255)

                        cTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCAtt                                                                >  1:4874073/1‑93 (MQ=255)

                         tGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCAttt                                                               >  1:1584373/1‑93 (MQ=255)

                          tgTTGTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATACAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTa                                                              <  2:1307578/92‑1 (MQ=255)

                          ggTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTa                                                              <  1:9730190/93‑1 (MQ=255)

                          ggTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTa                                                              <  2:4427138/93‑1 (MQ=255)

                          ggTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTa                                                              <  2:250513/93‑1 (MQ=255)

                          ggTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTa                                                              <  2:2425979/93‑1 (MQ=255)

                           gTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTAc                                                             <  1:3782268/93‑1 (MQ=255)

                            ttttGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACg                                                            <  2:8726044/93‑1 (MQ=255)

                            ttttGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACg                                                            <  2:4953089/93‑1 (MQ=255)

                            ttttGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACg                                                            >  1:8994867/1‑93 (MQ=255)

                            ttttGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACg                                                            >  2:3821212/1‑93 (MQ=255)

                            ttttGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACg                                                            <  1:8529264/93‑1 (MQ=255)

                             tttGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCGAGGATGGCACGGGTTTTATCTGTTGAACCGTCATTTACgg                                                           <  2:9207170/93‑1 (MQ=255)

                              ttGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGc                                                          <  1:6928593/93‑1 (MQ=255)

                               tGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCa                                                         <  1:1217641/93‑1 (MQ=255)

                                gcgcCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCaa                                                        <  1:782335/93‑1 (MQ=255)

                                gcgcCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCaa                                                        <  2:5472315/93‑1 (MQ=255)

                                  gcCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATa                                                      <  2:2228244/93‑1 (MQ=255)

                                   ccAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATaa                                                     <  2:504006/93‑1 (MQ=255)

                                     aGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACt                                                   >  2:7430970/1‑93 (MQ=255)

                                        tgaatgaCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCa                                                >  1:4357246/1‑93 (MQ=255)

                                        tgaatgaCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCa                                                <  1:3997729/93‑1 (MQ=255)

                                        tgaatgaCCCGCAAATGGGGAATTTGTGCAGACATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCa                                                >  1:5368489/1‑93 (MQ=255)

                                         gaatgaTCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCaa                                               <  2:7334751/93‑1 (MQ=255)

                                         gaatgaCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCaa                                               >  2:6026456/1‑93 (MQ=255)

                                         gaatgaCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCaa                                               <  1:9133797/93‑1 (MQ=255)

                                         gaatgaCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCaa                                               <  2:933263/93‑1 (MQ=255)

                                          aatgaCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAAt                                              >  2:456274/1‑93 (MQ=255)

                                           atgaCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATg                                             <  1:2974387/93‑1 (MQ=255)

                                             gaCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGtt                                           >  1:6857978/1‑93 (MQ=255)

                                              aCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCTCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTtc                                          >  1:2591154/1‑93 (MQ=255)

                                               ctcGCAAATGGGTAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTtct                                         <  2:5353967/91‑1 (MQ=255)

                                                ccGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTtctc                                        >  2:6803714/1‑93 (MQ=255)

                                                 cGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCa                                       <  2:7677087/93‑1 (MQ=255)

                                                 cGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCa                                       <  1:3812339/93‑1 (MQ=255)

                                                  gCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCAt                                      <  1:3669763/93‑1 (MQ=255)

                                                  gCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCAt                                      <  1:9320154/93‑1 (MQ=255)

                                                  gCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCAt                                      <  1:6152201/93‑1 (MQ=255)

                                                  gCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCAt                                      >  1:8376064/1‑93 (MQ=255)

                                                  gCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCAt                                      <  1:6314646/93‑1 (MQ=255)

                                                    aaaTGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATaa                                    <  1:1921175/93‑1 (MQ=255)

                                                    aaaTGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATaa                                    >  2:3908734/1‑93 (MQ=255)

                                                     agtGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAAc                                   <  1:8250410/91‑1 (MQ=255)

                                                     aaTGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAAc                                   <  1:4771965/93‑1 (MQ=255)

                                                       tGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGc                                 <  1:5489257/93‑1 (MQ=255)

                                                        ggggAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCt                                <  1:10000292/93‑1 (MQ=255)

                                                        ggggAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCt                                >  1:4223663/1‑93 (MQ=255)

                                                        ggggAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCt                                <  1:5226954/93‑1 (MQ=255)

                                                        ggggAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCt                                <  2:8462588/93‑1 (MQ=255)

                                                         gggAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCtg                               <  1:7179923/93‑1 (MQ=255)

                                                              tttGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTa                          >  2:1672589/1‑93 (MQ=255)

                                                               ttGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTaa                         <  1:2439278/93‑1 (MQ=255)

                                                                tgtgCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTaaa                        <  1:5627505/93‑1 (MQ=255)

                                                                 gtgCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAg                       <  2:5854633/93‑1 (MQ=255)

                                                                    cAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAgcgg                    <  1:9651960/93‑1 (MQ=255)

                                                                          tGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTa              >  1:2191237/1‑93 (MQ=255)

                                                                          tGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTa              >  2:8839685/1‑93 (MQ=255)

                                                                           gCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTAt             <  2:5672365/93‑1 (MQ=255)

                                                                           gCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTAt             <  2:7895074/93‑1 (MQ=255)

                                                                            cGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATg            >  2:469055/1‑93 (MQ=255)

                                                                             gATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCATTGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATgg           <  2:3561111/93‑1 (MQ=255)

                                                                             gATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATgg           <  2:7500111/93‑1 (MQ=255)

                                                                              aTCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGt          <  1:9914874/93‑1 (MQ=255)

                                                                               tCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGtt         >  1:6010621/1‑93 (MQ=255)

                                                                               tCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGtt         <  2:6527076/93‑1 (MQ=255)

                                                                                  aGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGTTTcc      >  2:8331036/1‑93 (MQ=255)

                                                                                    gATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGCGTATGGTGTCCTc    >  1:7283105/1‑93 (MQ=255)

                                                                                      tGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGTTTCCTCaa  >  2:3336622/1‑93 (MQ=255)

GGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCACGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGTTTCCTCAA  >  NC_000913/1086623‑1086801

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: `‑`