BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,999	A→G	T555A (ACT→GCT)	thrA →	aspartokinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,999	0	A	G	99.8%	1260.9 / NA	401	T555A (ACT→GCT)	thrA	aspartokinase I
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (193/207); minor base C (1/0); total (194/207)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 410 total aligned reads displayed.

CTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCG  >  NC_000913/1909‑2088

cTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCt                                                                                         <  1:103028/93‑1 (MQ=255)

 tGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCtt                                                                                        >  2:476554/1‑93 (MQ=255)

  ggCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTc                                                                                       <  2:712041/93‑1 (MQ=255)

   gcgcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTcc                                                                                      <  2:5885514/93‑1 (MQ=255)

    cgcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTcca                                                                                     <  1:1349833/93‑1 (MQ=255)

      cAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagc                                                                                   <  2:6602286/93‑1 (MQ=255)

      cAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagc                                                                                   <  2:3603337/93‑1 (MQ=255)

      cAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagc                                                                                   <  2:6950456/93‑1 (MQ=255)

        aGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagcca                                                                                 >  2:3662153/1‑93 (MQ=255)

        aGCCAAAGAGCCGTTTAATCTCGAGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagcca                                                                                 <  2:5656251/93‑1 (MQ=255)

         gCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagccag                                                                                <  2:675131/93‑1 (MQ=255)

         gCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagccag                                                                                <  1:4527396/93‑1 (MQ=255)

           cAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGc                                                                              <  2:2034906/93‑1 (MQ=255)

            aaaGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCa                                                                             <  2:6857986/93‑1 (MQ=255)

            aaaGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCa                                                                             >  1:6242887/1‑93 (MQ=255)

              agagCCGTTTAATCTCGGGCGCTTAATTCGCTTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGt                                                                           <  2:2277491/93‑1 (MQ=255)

              agagCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGt                                                                           <  2:5021670/93‑1 (MQ=255)

              agagCCGTTTAATCTCGGGCGCTTAATTAGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGt                                                                           >  1:1821742/1‑93 (MQ=255)

               gagCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTg                                                                          >  2:1820439/1‑93 (MQ=255)

                 gCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTTCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGc                                                                        >  1:1926349/1‑93 (MQ=255)

                 gCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGc                                                                        <  1:314795/93‑1 (MQ=255)

                 gCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGc                                                                        >  2:7299369/1‑93 (MQ=255)

                   cGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCgg                                                                      >  1:4948781/1‑93 (MQ=255)

                    gTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGa                                                                     >  1:7793690/1‑93 (MQ=255)

                    gTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGa                                                                     <  2:2368017/93‑1 (MQ=255)

                     tttAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCgggt                                                                    >  2:451720/1‑91 (MQ=255)

                     tttAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGAt                                                                    >  1:2420385/1‑93 (MQ=255)

                     tttAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGAt                                                                    >  2:4092573/1‑93 (MQ=255)

                      ttAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATc                                                                   <  1:5558255/93‑1 (MQ=255)

                       tAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCa                                                                  >  1:3673158/1‑93 (MQ=255)

                        aaTCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCaa                                                                 <  1:5139955/93‑1 (MQ=255)

                         aTCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAat                                                                >  1:6380347/1‑93 (MQ=255)

                         aTCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAat                                                                <  2:5989873/93‑1 (MQ=255)

                          tctcGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAata                                                               <  1:935465/93‑1 (MQ=255)

                          tctcGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAata                                                               <  1:755495/93‑1 (MQ=255)

                          tctcGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAata                                                               <  2:1500949/93‑1 (MQ=255)

                           ctcGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAatat                                                              <  2:6473920/93‑1 (MQ=255)

                           ctcGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAatat                                                              <  1:3114920/93‑1 (MQ=255)

                             cGGGCGCTTAATTCGCCTCGTGAATGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGc                                                            <  1:7026955/93‑1 (MQ=255)

                             cGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGc                                                            >  1:4862648/1‑93 (MQ=255)

                               ggCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCg                                                          >  2:352130/1‑93 (MQ=255)

                               ggCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCg                                                          >  2:2613158/1‑93 (MQ=255)

                               ggCGCTTAATACGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCg                                                          >  2:4318302/1‑93 (MQ=255)

                                gcgcTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGa                                                         <  2:7390392/93‑1 (MQ=255)

                                 cgcTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGAc                                                        <  2:6656778/93‑1 (MQ=255)

                                  gcTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTTATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACt                                                       >  2:7764271/1‑93 (MQ=255)

                                  gcTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACt                                                       >  2:2730402/1‑93 (MQ=255)

                                  gcTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACt                                                       >  1:5786303/1‑93 (MQ=255)

                                    ttAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTc                                                     >  2:2995829/1‑93 (MQ=255)

                                    gtAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTc                                                     <  1:5394094/92‑1 (MQ=255)

                                       aTTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTg                                                  <  2:30659/93‑1 (MQ=255)

                                          cGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTgcgc                                               <  2:2809519/93‑1 (MQ=255)

                                           gCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGcgcg                                              <  2:4891329/93‑1 (MQ=255)

                                            ccTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGa                                             <  1:2977810/93‑1 (MQ=255)

                                             cTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGaa                                            >  1:3919325/1‑93 (MQ=255)

                                             cTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGaa                                            <  2:1695347/93‑1 (MQ=255)

                                              tCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAg                                           <  1:4103257/93‑1 (MQ=255)

                                                 tGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGtt                                        >  1:6832338/1‑93 (MQ=255)

                                                 tGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGtt                                        >  1:1500685/1‑93 (MQ=255)

                                                  gAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGttt                                       >  1:2727721/1‑93 (MQ=255)

                                                     aGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCa                                    <  1:1723625/93‑1 (MQ=255)

                                                     aGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCa                                    <  1:6596566/93‑1 (MQ=255)

                                                     aGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCa                                    <  2:7705680/93‑1 (MQ=255)

                                                       aaTATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACg                                  >  2:3451678/1‑93 (MQ=255)

                                                        atatCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGt                                 >  2:4983406/1‑93 (MQ=255)

                                                        atatCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGt                                 <  1:7925531/93‑1 (MQ=255)

                                                        atatCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGt                                 >  1:2119641/1‑93 (MQ=255)

                                                          atcatcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTg                               <  1:6907009/93‑1 (MQ=255)

                                                           tcatcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGt                              >  1:7526028/1‑93 (MQ=255)

                                                           tcatcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGt                              <  1:2918323/93‑1 (MQ=255)

                                                             atcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCa                            <  1:1175391/93‑1 (MQ=255)

                                                             atcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCa                            <  2:7860406/93‑1 (MQ=255)

                                                             atcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCa                            <  1:3366473/93‑1 (MQ=255)

                                                             atcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCa                            >  1:6006425/1‑93 (MQ=255)

                                                             atcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCa                            >  1:6175294/1‑93 (MQ=255)

                                                              tcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCAc                           <  1:6470309/93‑1 (MQ=255)

                                                                 gctgAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGcc                        <  2:3925191/93‑1 (MQ=255)

                                                                   tgAACCCGGTCATTGTTGACTGCGCTTCGAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGa                      >  2:7569308/1‑93 (MQ=255)

                                                                    gAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGaa                     <  1:5114940/93‑1 (MQ=255)

                                                                     aaCCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAAc                    <  1:7235892/93‑1 (MQ=255)

                                                                      aCCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACa                   >  2:2455462/1‑93 (MQ=255)

                                                                        ccGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACaaa                 >  1:3498839/1‑93 (MQ=255)

                                                                         cGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACaaaa                <  1:4203018/93‑1 (MQ=255)

                                                                         cGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACaaaa                <  1:5235654/93‑1 (MQ=255)

                                                                         cGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACaaaa                >  1:533469/1‑93 (MQ=255)

                                                                         cGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACaaaa                >  2:3156311/1‑93 (MQ=255)

                                                                          ggTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACaaaaa               <  2:5137153/93‑1 (MQ=255)

                                                                          ggTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACaaaaa               <  1:5885327/93‑1 (MQ=255)

                                                                          ggTCATTGTTGACTGCGCTGACAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACaaaaa               <  2:4433561/93‑1 (MQ=255)

                                                                           gTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAg              >  2:6221221/1‑93 (MQ=255)

                                                                            tCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAgg             >  1:732207/1‑93 (MQ=255)

                                                                             cATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGc            <  2:1100460/93‑1 (MQ=255)

                                                                              aTTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGcc           >  1:4570634/1‑93 (MQ=255)

                                                                               ttgttgACTGCGCTTCCGGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGGAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCa          >  2:4625248/1‑93 (MQ=255)

                                                                               ttgttgACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCa          >  1:6801479/1‑93 (MQ=255)

                                                                               ttgttgACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCa          >  2:2200970/1‑93 (MQ=255)

                                                                                tgttgACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCaa         <  2:14023/93‑1 (MQ=255)

                                                                                    gACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACAcc     <  2:5210458/93‑1 (MQ=255)

                                                                                     aCTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCt    <  2:5841319/93‑1 (MQ=255)

                                                                                      cTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCtc   <  2:5446617/93‑1 (MQ=255)

CTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCG  >  NC_000913/1909‑2088

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