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.6%	851.7 / ‑6.6	265	V462A (GTA→GCA)	metL	aspartokinase II
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (141/123); minor base A (0/1); total (141/124)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 270 total aligned reads displayed.

TGGTGGCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGT  >  NC_000913/4128709‑4128888

 ggtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCAtt                                                                                        >  2:499396/1‑93 (MQ=255)

 ggtggCAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCAtt                                                                                        >  1:5051410/1‑93 (MQ=255)

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

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

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

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

      cAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCg                                                                                   >  2:5880510/1‑93 (MQ=255)

      cAGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCg                                                                                   <  1:3988331/93‑1 (MQ=255)

       aGTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCgg                                                                                  >  1:4425833/1‑93 (MQ=255)

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

        gTACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGt                                                                                 >  1:3738824/1‑93 (MQ=255)

         tACTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTa                                                                                >  2:633992/1‑93 (MQ=255)

          aCTGCGCACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTaa                                                                               <  1:4266898/93‑1 (MQ=255)

            tGCGCGCCGCCCACACCGAAAGCCTGATTGAGTGGCTGCATAAGTCCGTCTGCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAgg                                                                             <  1:2179400/93‑1 (MQ=255)

              cgcACCGCCCCGTCCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGTCATTGTTCGGTAAGGGc                                                                           <  2:1034528/93‑1 (MQ=255)

               gcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTCAGGGCa                                                                          >  2:3218443/1‑93 (MQ=255)

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

               gcACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTATTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCa                                                                          >  1:4987530/1‑93 (MQ=255)

                cACCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCaa                                                                         <  1:5307641/93‑1 (MQ=255)

                 aCCGGCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAat                                                                        <  2:3377537/93‑1 (MQ=255)

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

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

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

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

                    cgcccgcccGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATc                                                                     <  1:1640286/92‑1 (MQ=255)

                     gCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCg                                                                    <  1:4031126/93‑1 (MQ=255)

                     gCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCg                                                                    <  2:1610285/93‑1 (MQ=255)

                     gCCCGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCg                                                                    >  2:2312457/1‑93 (MQ=255)

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

                      cccGACCGAAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCgg                                                                   >  1:5633458/1‑93 (MQ=255)

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

                       ccgaccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGACTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGt                                                                  <  1:1155382/93‑1 (MQ=255)

                        cgaccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGtt                                                                 >  2:2137823/1‑93 (MQ=255)

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

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

                          accgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTcc                                                               >  1:3629639/1‑93 (MQ=255)

                          accgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTcc                                                               >  1:1696609/1‑93 (MQ=255)

                           ccgaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTccc                                                              >  1:467068/1‑93 (MQ=255)

                            cgaAAGCCTGATTCAGGGGCTGCACCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCg                                                             >  1:769034/1‑93 (MQ=255)

                             gaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCGGAAAAACGCATCGGCCTGGCATTGTTCGGTTAGGGCAATATTGGTTcccct                                                            >  1:71577/1‑91 (MQ=255)

                             gaAAGCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGt                                                            >  2:4287135/1‑93 (MQ=255)

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

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

                                 gCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGc                                                        >  2:2918954/1‑93 (MQ=255)

                                 gCCTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGc                                                        >  2:5335610/1‑93 (MQ=255)

                                  ccTGATTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCt                                                       >  2:5842089/1‑93 (MQ=255)

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

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

                                      aTTCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGaa                                                   <  2:1462740/93‑1 (MQ=255)

                                       ttCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAAc                                                  >  2:6218599/1‑93 (MQ=255)

                                       ttCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAAc                                                  >  2:385772/1‑93 (MQ=255)

                                       ttCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAAc                                                  >  2:1240086/1‑93 (MQ=255)

                                        tCAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACt                                                 >  2:458664/1‑93 (MQ=255)

                                         cAGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTg                                                >  2:2227258/1‑93 (MQ=255)

                                          aGGGGCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTGCCGTTGGCTGGAACTGt                                               >  2:6288825/1‑93 (MQ=255)

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

                                           ggggCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGtt                                              <  2:4190007/93‑1 (MQ=255)

                                           ggggCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGtt                                              <  1:2330746/93‑1 (MQ=255)

                                           ggggCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGtt                                              <  1:2843185/93‑1 (MQ=255)

                                           ggggCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGtt                                              <  2:3924908/93‑1 (MQ=255)

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

                                             ggCTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCg                                            >  1:2279507/1‑93 (MQ=255)

                                               cTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGcc                                          <  2:3163969/93‑1 (MQ=255)

                                               cTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGcc                                          <  2:6135909/93‑1 (MQ=255)

                                               cTGCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGcc                                          <  2:4050207/93‑1 (MQ=255)

                                                 gCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCg                                        <  2:1768673/93‑1 (MQ=255)

                                                 gCATCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGCTCGCCCg                                        >  1:2719854/1‑93 (MQ=255)

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

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

                                                    tCAGTCCGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTga                                     >  1:5436872/1‑93 (MQ=255)

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

                                                         ccGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcaga                                <  2:1973108/93‑1 (MQ=255)

                                                         ccGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcaga                                <  1:1709054/93‑1 (MQ=255)

                                                          cGTCTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTgagcagag                               >  1:721065/1‑93 (MQ=255)

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

                                                             cTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCAc                            <  1:4568730/93‑1 (MQ=255)

                                                             cTTCCGCGCAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCAc                            <  1:5207490/93‑1 (MQ=255)

                                                                  gcgcAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCttt                       >  1:2044733/1‑93 (MQ=255)

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

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

                                                                     cAGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCgg                    <  2:3470037/93‑1 (MQ=255)

                                                                      aGAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGc                   >  2:1551443/1‑93 (MQ=255)

                                                                       gAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCa                  >  2:3787588/1‑93 (MQ=255)

                                                                       gAAAAACGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCa                  >  2:4549772/1‑93 (MQ=255)

                                                                        aaaacCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCAc                 <  1:794704/93‑1 (MQ=255)

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

                                                                           aaCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTa              <  1:647523/93‑1 (MQ=255)

                                                                            aCGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTAc             >  1:1039863/1‑93 (MQ=255)

                                                                             cGCATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTAcc            <  1:2944276/93‑1 (MQ=255)

                                                                               cATCGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCTGCACGTACCgg          <  1:3164784/93‑1 (MQ=255)

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

                                                                                  cGGCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCtt       <  1:1537190/93‑1 (MQ=255)

                                                                                   ggCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTCCCGGCttt      >  1:2126631/1‑93 (MQ=255)

                                                                                    gCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTg     <  1:241662/93‑1 (MQ=255)

                                                                                    gCCTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTg     >  2:1722205/1‑93 (MQ=255)

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

                                                                                      cTGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAg   >  2:5049546/1‑93 (MQ=255)

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

                                                                                       tGGCATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGt  >  2:2700755/1‑93 (MQ=255)

                                                                                       tGGCATTGTTCGGGAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCCGTGAGCAGAGCACGCTTTCGGCACGTACCGGCTTTGAGt  >  2:3186512/1‑93 (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` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`