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.5%	618.9 / ‑4.9	192	T555A (ACT→GCT)	thrA	aspartokinase I
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (103/88); total (104/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

Only 100 of 193 total aligned reads displayed.

CTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCG  >  NC_000913/1909‑2088

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

  ggCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTc                                                                                       >  2:4279387/1‑93 (MQ=255)

  ggCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTc                                                                                       >  2:1215881/1‑93 (MQ=255)

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

  ggCCCAAGCCAAAGAGCCGTGTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTc                                                                                       <  1:2857189/93‑1 (MQ=255)

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

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

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

    ggcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTcca                                                                                     <  2:658333/92‑1 (MQ=255)

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

     gcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccag                                                                                    >  2:3332520/1‑93 (MQ=255)

       aaGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagcc                                                                                  >  2:1034364/1‑93 (MQ=255)

       aaGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagcc                                                                                  <  1:2214232/93‑1 (MQ=255)

       aaGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagcc                                                                                  <  2:2927034/93‑1 (MQ=255)

       aaGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTccagcc                                                                                  <  1:2831336/93‑1 (MQ=255)

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

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

           cAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGc                                                                              >  1:3777704/1‑93 (MQ=255)

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

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

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

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

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

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

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

                  ccGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCg                                                                       >  2:3032425/1‑93 (MQ=255)

                  ccGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCg                                                                       <  1:2407478/93‑1 (MQ=255)

                   cGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCgg                                                                      <  2:1961828/93‑1 (MQ=255)

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

                   cGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCgg                                                                      <  1:1306996/93‑1 (MQ=255)

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

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

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

                        aaTCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCaa                                                                 >  2:3232224/1‑93 (MQ=255)

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

                                   cTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACtt                                                      <  2:4535513/93‑1 (MQ=255)

                                   cTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACtt                                                      >  2:3112761/1‑93 (MQ=255)

                                   cTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGAAtt                                                      >  2:375417/1‑93 (MQ=255)

                                    ttAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTc                                                     <  2:4551477/93‑1 (MQ=255)

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

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

                                    ttAATTCGCCTCGTGAAAGAATAGCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTc                                                     <  2:3395690/93‑1 (MQ=255)

                                      aaTTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCt                                                   <  2:414736/93‑1 (MQ=255)

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

                                       aTTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTg                                                  >  2:774853/1‑93 (MQ=255)

                                        ttCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTgc                                                 <  1:508078/93‑1 (MQ=255)

                                        ttCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTgc                                                 <  2:3081695/93‑1 (MQ=255)

                                         tCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTgcg                                                <  1:3665813/93‑1 (MQ=255)

                                           gCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGcgcg                                              >  1:3066066/1‑93 (MQ=255)

                                           gCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGcgcg                                              >  1:3644383/1‑93 (MQ=255)

                                           gCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGcgcg                                              >  1:2120711/1‑93 (MQ=255)

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

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

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

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

                                               cGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAgg                                          <  1:2392222/93‑1 (MQ=255)

                                               cGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAgg                                          <  2:1583542/93‑1 (MQ=255)

                                                gTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGt                                         <  1:1625628/93‑1 (MQ=255)

                                                 tGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGGTCAATATGCCGATTTCTTGCGCGAAGGtt                                        >  2:3589301/1‑93 (MQ=255)

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

                                                 tGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGtt                                        <  2:1775318/93‑1 (MQ=255)

                                                 tGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGtt                                        <  1:3707581/93‑1 (MQ=255)

                                                   aaaGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTc                                      <  2:3779991/93‑1 (MQ=255)

                                                    aaGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTcc                                     >  1:1187222/1‑93 (MQ=255)

                                                    aaGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTcc                                     >  1:25005/1‑93 (MQ=255)

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

                                                     aGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCa                                    >  1:9059/1‑93 (MQ=255)

                                                       aaTATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACg                                  <  1:4714318/93‑1 (MQ=255)

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

                                                         tatCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGtt                                <  2:3658776/93‑1 (MQ=255)

                                                         tatCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGtt                                <  2:4317458/93‑1 (MQ=255)

                                                          atcatcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCGATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTg                               >  1:4376377/1‑93 (MQ=255)

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

                                                             atcTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTcc                            >  1:339223/1‑92 (MQ=255)

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

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

                                                                 gctgAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGcc                        >  2:73915/1‑93 (MQ=255)

                                                                  ctgAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACCCCg                       >  1:322290/1‑93 (MQ=255)

                                                                   tgAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGa                      >  2:4676737/1‑93 (MQ=255)

                                                                   tgAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGa                      >  1:3664607/1‑93 (MQ=255)

                                                                    gAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGaa                     >  2:1967124/1‑93 (MQ=255)

                                                                     aaCCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGaaa                    >  1:3886836/1‑92 (MQ=255)

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

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

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

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

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

                                                                           gTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAg              <  2:173088/93‑1 (MQ=255)

                                                                            tCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGCTCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAgg             >  1:2402852/1‑93 (MQ=255)

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

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

                                                                               ttgttgACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACTCCGAACAAAAAGGCCa          >  1:3039052/1‑93 (MQ=255)

                                                                                ttttgACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCaa         <  2:4353780/91‑1 (MQ=255)

                                                                                  ttgACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAaca       <  2:1184609/93‑1 (MQ=255)

                                                                                  ttgACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAaca       <  1:147605/93‑1 (MQ=255)

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

                                                                                      cTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCtc   >  1:1116083/1‑93 (MQ=255)

                                                                                      cTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCtc   >  2:676536/1‑93 (MQ=255)

                                                                                       tGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCtcg  <  2:4175785/93‑1 (MQ=255)

                                                                                       tGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCtcg  <  2:2626973/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` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`