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	98.3%	563.1 / ‑1.5	180	R119R (CGT→CGC)	ycdQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (72/104); total (74/106)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 184 total aligned reads displayed.

GTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCACGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGTTTCCTCAAC  >  NC_000913/1086624‑1086802

gTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCgg                                                                                        <  1:1648807/93‑1 (MQ=255)

gTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCgg                                                                                        >  2:4139602/1‑93 (MQ=255)

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

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

  tttAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGt                                                                                      >  1:4431116/1‑93 (MQ=255)

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

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

      aGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTa                                                                                  >  2:1620084/1‑93 (MQ=255)

      aGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTa                                                                                  <  2:345597/93‑1 (MQ=255)

       gcgcAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCACGGGTTTTAt                                                                                 >  2:1568837/1‑93 (MQ=255)

         gcAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCt                                                                               <  1:3869905/93‑1 (MQ=255)

         gcAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCt                                                                               <  1:3347759/93‑1 (MQ=255)

           aaTGTCTTTCCCCTGGTTTTGCGCCAGATGAATGGCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGt                                                                             <  2:3854863/93‑1 (MQ=255)

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

               gCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGaa                                                                         >  2:2224196/1‑93 (MQ=255)

               gCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGaa                                                                         <  1:2789414/93‑1 (MQ=255)

               gCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGaa                                                                         <  2:2427126/93‑1 (MQ=255)

               gCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCACGGGTTTTATCTGTTGaa                                                                         >  2:745570/1‑93 (MQ=255)

                cTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAAc                                                                        >  2:1624436/1‑93 (MQ=255)

                  ttCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCg                                                                      <  1:2550704/93‑1 (MQ=255)

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

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

                   tCCCCTGGTTTTGCGCCATATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATGTGTTGAACCGt                                                                     <  2:405397/93‑1 (MQ=255)

                   tCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGt                                                                     >  1:466267/1‑93 (MQ=255)

                   tCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGt                                                                     <  1:396562/93‑1 (MQ=255)

                      ccTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCAt                                                                  >  1:392199/1‑93 (MQ=255)

                      ccTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCAt                                                                  >  1:3421359/1‑93 (MQ=255)

                      ccTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCAt                                                                  >  1:2984600/1‑93 (MQ=255)

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

                       cTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGGCAtt                                                                 >  1:932405/1‑93 (MQ=255)

                        tGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCAttt                                                                <  1:751859/93‑1 (MQ=255)

                        tGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCAttt                                                                <  1:1012483/93‑1 (MQ=255)

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

                          gTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTAc                                                              >  1:3593496/1‑93 (MQ=255)

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

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

                          gTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTAc                                                              >  1:2760263/1‑93 (MQ=255)

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

                            tttGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACgg                                                            >  2:2493932/1‑93 (MQ=255)

                             ttGCGCCAGAGGAATGACCCGCAGATGGGGAATTGGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGc                                                           <  2:2584695/93‑1 (MQ=255)

                              tGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCa                                                          >  1:331523/1‑93 (MQ=255)

                                cgcCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAAt                                                        >  1:2024393/1‑93 (MQ=255)

                                 gcCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATa                                                       >  2:3725905/1‑93 (MQ=255)

                                 gcCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATa                                                       >  2:240239/1‑93 (MQ=255)

                                   cAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAAc                                                     <  1:483615/93‑1 (MQ=255)

                                    aGATGAAGGACCCGAAAATGGGTATTTTGTGGAGCCAGGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACt                                                    <  2:3303207/93‑1 (MQ=255)

                                     gATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACtt                                                   <  2:1718486/93‑1 (MQ=255)

                                     gATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACtt                                                   >  2:1042254/1‑93 (MQ=255)

                                      atgaatgaCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTc                                                  <  1:4118739/93‑1 (MQ=255)

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

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

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

                                           ttaCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGt                                             <  2:3263605/91‑1 (MQ=255)

                                           tgaCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGt                                             >  1:4354126/1‑93 (MQ=255)

                                           tgaCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGt                                             <  1:2287000/93‑1 (MQ=255)

                                             aCCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTtc                                           <  2:3452658/93‑1 (MQ=255)

                                             aCCCGCAAATGGGGAATTTGTGCAGACATGCGATCCCGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATACCTTCCATGTtc                                           >  2:3329691/1‑93 (MQ=255)

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

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

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

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

                                                  cAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATa                                      >  2:1955127/1‑93 (MQ=255)

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

                                                    aaTGGTGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAAc                                    <  1:4230819/93‑1 (MQ=255)

                                                    aaTGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAAc                                    >  2:3632876/1‑93 (MQ=255)

                                                     aTGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACg                                   <  2:1178413/93‑1 (MQ=255)

                                                     aTGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACg                                   >  2:1828438/1‑93 (MQ=255)

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

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

                                                        gggAATTTGTGCAGCCATGCGATCCAGGATGGCGAGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCtg                                <  2:952848/93‑1 (MQ=255)

                                                         ggAATTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCtgt                               <  2:516316/93‑1 (MQ=255)

                                                           aaTTTGTGCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGc                             <  1:789831/93‑1 (MQ=255)

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

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

                                                                ttgCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAg                        <  2:3562617/92‑1 (MQ=255)

                                                                gtgCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTACAg                        >  1:1234729/1‑93 (MQ=255)

                                                                 tgCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAgc                       <  2:2507872/93‑1 (MQ=255)

                                                                 tgCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAgc                       >  1:2494424/1‑93 (MQ=255)

                                                                  gCAGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAgcg                      <  2:1930862/93‑1 (MQ=255)

                                                                    aGCCATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAgcggc                    >  1:2693772/1‑93 (MQ=255)

                                                                       cATGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCgtg                 >  1:1710706/1‑93 (MQ=255)

                                                                        aTGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCgtgt                <  1:4089157/93‑1 (MQ=255)

                                                                         tGCGATCCAGGATGGCTCGGTGTTTATCTGTTGAACCGTCAGTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTa               <  1:253875/93‑1 (MQ=255)

                                                                         tGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTa               <  2:777013/93‑1 (MQ=255)

                                                                         tGCGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTa               <  1:1422089/93‑1 (MQ=255)

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

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

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

                                                                           cGATCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATg             <  1:606663/93‑1 (MQ=255)

                                                                             aTCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTAGGGCAATAACTTCAATGGTCTCATAACGCTGTGCTAAAGCGGCGTGGATGGt           >  1:200082/1‑93 (MQ=255)

                                                                             aTCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGt           >  2:3970083/1‑93 (MQ=255)

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

                                                                              tCCAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACGTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGtt          >  1:3230980/1‑93 (MQ=255)

                                                                              gccAGGATGGCGCGGGTTTTATCTGTTGAAGCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGtt          <  2:3106352/92‑1 (MQ=255)

                                                                               ccAGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGttt         <  1:3095122/93‑1 (MQ=255)

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

                                                                                 aGGATGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGTTTcc       <  2:4200301/93‑1 (MQ=255)

                                                                                    aTGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGTTTCCTCa    >  2:260559/1‑93 (MQ=255)

                                                                                     tGGCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGTTTCCTCaa   <  2:1798089/93‑1 (MQ=255)

                                                                                      ggCGCGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGTTTCCTCAAc  <  1:496026/93‑1 (MQ=255)

GTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGATGGCACGGGTTTTATCTGTTGAACCGTCATTTACGGCAATAACTTCAATGTTCTCATAACGCTGTGCTAAAGCGGCGTGTATGGTTTCCTCAAC  >  NC_000913/1086624‑1086802

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