BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,799	T→C	V462A (GTA→GCA)	metL →	aspartokinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,799	0	T	C	99.7%	985.2 / ‑4.9	301	V462A (GTA→GCA)	metL	aspartokinase II
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (141/159); minor base A (1/0); total (142/159)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 306 total aligned reads displayed.

TGGTGGCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGT  >  NC_000913/4128709‑4128888

tggtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCAt                                                                                         <  1:1467633/93‑1 (MQ=255)

tggtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCAt                                                                                         >  1:4125768/1‑93 (MQ=255)

  gtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTg                                                                                       >  1:4736006/1‑93 (MQ=255)

  gtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTg                                                                                       <  2:1833325/93‑1 (MQ=255)

   tggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGt                                                                                      <  2:2532918/93‑1 (MQ=255)

   tggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGt                                                                                      >  1:3456850/1‑93 (MQ=255)

    ggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGtt                                                                                     <  1:2484926/93‑1 (MQ=255)

       aGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCgg                                                                                  <  1:1225895/93‑1 (MQ=255)

          aCTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTaa                                                                               >  1:1175269/1‑93 (MQ=255)

          aCTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTaa                                                                               >  2:2609264/1‑93 (MQ=255)

          aCTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGATaa                                                                               >  2:3680936/1‑93 (MQ=255)

           cTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAg                                                                              >  1:99467/1‑93 (MQ=255)

           cTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAg                                                                              <  2:478294/93‑1 (MQ=255)

            tGCGCGCCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAgg                                                                             <  2:1224536/93‑1 (MQ=255)

            tGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAgg                                                                             >  2:4522968/1‑93 (MQ=255)

            tGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAgg                                                                             >  2:160545/1‑93 (MQ=255)

            tGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAgg                                                                             <  2:4972421/93‑1 (MQ=255)

             gcgcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAggg                                                                            <  2:5306872/93‑1 (MQ=255)

               gcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCa                                                                          >  1:177220/1‑93 (MQ=255)

               gcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCa                                                                          <  1:1124407/93‑1 (MQ=255)

                cACCGGCCCGACCGAAAGCGTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCaa                                                                         <  2:6447728/93‑1 (MQ=255)

                 aCCGGCCCGACCGAAAGCCTGATTCAGGTGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAat                                                                        <  2:3745339/93‑1 (MQ=255)

                 aCCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAat                                                                        >  2:3219496/1‑93 (MQ=255)

                 aCCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAat                                                                        >  1:3099237/1‑93 (MQ=255)

                 aCCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAat                                                                        >  2:1987833/1‑93 (MQ=255)

                  ccGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAata                                                                       <  2:2446013/93‑1 (MQ=255)

                  ccGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAata                                                                       >  2:2151187/1‑93 (MQ=255)

                  ccGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAata                                                                       >  2:672692/1‑93 (MQ=255)

                   cGGCCCGCCCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAatat                                                                      <  1:2629868/93‑1 (MQ=255)

                   cGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAatat                                                                      <  2:2326887/93‑1 (MQ=255)

                    ggCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATc                                                                     >  2:4544889/1‑93 (MQ=255)

                    ggCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATc                                                                     <  1:6231580/93‑1 (MQ=255)

                    ggCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATc                                                                     <  1:5065166/93‑1 (MQ=255)

                    ggCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATc                                                                     >  1:2019196/1‑93 (MQ=255)

                    ggCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATc                                                                     >  2:5620630/1‑93 (MQ=255)

                      cccGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCgg                                                                   <  1:5464552/93‑1 (MQ=255)

                       ccgaccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGt                                                                  <  1:4986719/93‑1 (MQ=255)

                       ccgaccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGt                                                                  >  1:4002833/1‑93 (MQ=255)

                        cgaccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGtt                                                                 <  2:547464/93‑1 (MQ=255)

                        cgaccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGtt                                                                 <  1:3518542/93‑1 (MQ=255)

                         gaccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTc                                                                <  2:5547150/93‑1 (MQ=255)

                         gaccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTc                                                                >  2:6119595/1‑93 (MQ=255)

                             gaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGt                                                            <  2:344091/93‑1 (MQ=255)

                              aaaGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGtt                                                           <  1:218682/93‑1 (MQ=255)

                              aaaGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGtt                                                           >  1:4093489/1‑93 (MQ=255)

                              aaaGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGtt                                                           <  1:1728304/93‑1 (MQ=255)

                               aaGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTg                                                          <  1:2262985/93‑1 (MQ=255)

                                aGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTgg                                                         <  1:3222914/93‑1 (MQ=255)

                                aGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTgg                                                         >  1:2056997/1‑93 (MQ=255)

                                  ccTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCt                                                       <  1:1855599/93‑1 (MQ=255)

                                   cTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTg                                                      >  1:5310105/1‑93 (MQ=255)

                                   cTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTg                                                      <  1:2332505/93‑1 (MQ=255)

                                   cTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTg                                                      <  2:4896471/93‑1 (MQ=255)

                                   cTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTg                                                      >  1:2782042/1‑93 (MQ=255)

                                    tGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTgg                                                     >  2:2735863/1‑93 (MQ=255)

                                      aTTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGAATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGaa                                                   >  1:840833/1‑93 (MQ=255)

                                       ttCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAAc                                                  <  1:3701524/93‑1 (MQ=255)

                                          aGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGt                                               >  1:5951185/1‑93 (MQ=255)

                                          aGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGt                                               <  1:5602162/93‑1 (MQ=255)

                                            gggCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTc                                             <  1:6487078/93‑1 (MQ=255)

                                             ggCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCg                                            <  2:6567637/93‑1 (MQ=255)

                                                tGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGccc                                         <  1:4609244/93‑1 (MQ=255)

                                                tGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGccc                                         >  1:4938611/1‑93 (MQ=255)

                                                tGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGccc                                         >  2:6388601/1‑93 (MQ=255)

                                                  cATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGt                                       >  2:5222591/1‑93 (MQ=255)

                                                  cATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGt                                       <  2:3983446/93‑1 (MQ=255)

                                                   aTCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTg                                      >  2:6195114/1‑93 (MQ=255)

                                                     cAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgag                                    <  2:434367/93‑1 (MQ=255)

                                                      aGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagc                                   <  2:2987371/93‑1 (MQ=255)

                                                       gTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagca                                  >  1:2981424/1‑93 (MQ=255)

                                                       gTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagca                                  >  2:58867/1‑93 (MQ=255)

                                                        tCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcag                                 <  2:610179/93‑1 (MQ=255)

                                                           gTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagagc                              >  1:5788805/1‑93 (MQ=255)

                                                           gTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagagc                              <  2:2258880/93‑1 (MQ=255)

                                                           gTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagagc                              >  1:4396390/1‑93 (MQ=255)

                                                             cTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCAc                            >  1:679837/1‑93 (MQ=255)

                                                               tCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGTGACTGTTCGCCCGTGAGCCGAGGACGc                          >  2:4436786/1‑93 (MQ=255)

                                                               tCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGc                          <  2:4041300/93‑1 (MQ=255)

                                                               tCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGc                          >  2:5100374/1‑93 (MQ=255)

                                                                 cgcgCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCtt                        >  2:5767266/1‑93 (MQ=255)

                                                                 cgcgCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCtt                        <  2:109086/93‑1 (MQ=255)

                                                                  gcgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCttt                       <  2:1757065/93‑1 (MQ=255)

                                                                  gcgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCttt                       <  1:5568992/93‑1 (MQ=255)

                                                                    gcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCg                     <  1:2703319/93‑1 (MQ=255)

                                                                    gcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCg                     >  1:3925944/1‑93 (MQ=255)

                                                                      aGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGc                   <  2:3606204/93‑1 (MQ=255)

                                                                      aGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGc                   <  2:3139876/93‑1 (MQ=255)

                                                                         aaaaCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACg                >  2:1478567/1‑93 (MQ=255)

                                                                          aaaCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGt               >  2:3916094/1‑93 (MQ=255)

                                                                               cATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCgg          >  1:2426913/1‑93 (MQ=255)

                                                                               cATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCgg          >  2:757848/1‑93 (MQ=255)

                                                                                 tCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCt        <  2:4719056/93‑1 (MQ=255)

                                                                                 tCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCt        >  1:6649356/1‑93 (MQ=255)

                                                                                   ggCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCttt      >  2:2817226/1‑93 (MQ=255)

                                                                                     ccTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGa    <  2:4956511/93‑1 (MQ=255)

                                                                                     ccTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGa    <  2:3026879/93‑1 (MQ=255)

                                                                                      cTGGCATTGTTAGGTAAGGGCAATATCGGTTCCCGTTGTCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAg   <  2:4587109/93‑1 (MQ=255)

                                                                                       tGGCATTTTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGt  <  2:1301395/93‑1 (MQ=255)

                                                                                       tGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGt  <  2:1385662/93‑1 (MQ=255)

                                                                                       tGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGt  <  1:4540830/93‑1 (MQ=255)

TGGTGGCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGT  >  NC_000913/4128709‑4128888

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