BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,055,398	A→G	S73P (TCC→CCC)	glnA ←	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,055,398	0	A	G	99.6%	1003.2 / ‑4.6	278	S73P (TCC→CCC)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (130/147); total (131/147)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

Only 100 of 299 total aligned reads displayed.

GGTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGA  >  NC_000913/4055308‑4055487

ggTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCggtgggtg                                                                                         >  2:3503383/1‑93 (MQ=255)

ggTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCggtgggtg                                                                                         >  2:5933290/1‑93 (MQ=255)

ggTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCggtgggtg                                                                                         <  2:6228805/93‑1 (MQ=255)

  tCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCg                                                                                       <  2:1749487/93‑1 (MQ=255)

    aTAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTTGGTGCGTc                                                                                     >  1:6174155/1‑93 (MQ=255)

    aTAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTc                                                                                     <  1:1190585/93‑1 (MQ=255)

     tAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCt                                                                                    >  2:2163970/1‑93 (MQ=255)

     tAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCt                                                                                    >  2:3732640/1‑93 (MQ=255)

      aGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTg                                                                                   <  2:684941/93‑1 (MQ=255)

      aGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTg                                                                                   <  1:4620496/93‑1 (MQ=255)

         cTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCa                                                                                >  1:5515422/1‑93 (MQ=255)

          ttGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCAt                                                                               >  2:4535177/1‑93 (MQ=255)

          ttGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCAt                                                                               >  2:1314965/1‑93 (MQ=255)

           tGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATc                                                                              >  2:3258413/1‑93 (MQ=255)

           tGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATc                                                                              <  2:4695853/93‑1 (MQ=255)

            gCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCa                                                                             <  1:5889578/93‑1 (MQ=255)

            gCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCa                                                                             >  2:4218665/1‑93 (MQ=255)

             cAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAg                                                                            <  1:2325723/93‑1 (MQ=255)

             cAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAg                                                                            <  2:6092736/93‑1 (MQ=255)

             cAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAg                                                                            <  1:580926/93‑1 (MQ=255)

              aGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGc                                                                           >  2:6534323/1‑93 (MQ=255)

                     cAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGt                                                                    >  2:5579920/1‑93 (MQ=255)

                     cAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGt                                                                    <  1:6043673/93‑1 (MQ=255)

                        gTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCTTCAGCACCATGTCgg                                                                 >  1:910740/1‑93 (MQ=255)

                        gTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCgg                                                                 >  1:1311315/1‑93 (MQ=255)

                         ttCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGa                                                                >  2:6486197/1‑93 (MQ=255)

                          tCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGAc                                                               <  1:3489954/93‑1 (MQ=255)

                          tCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGAc                                                               >  2:3986032/1‑93 (MQ=255)

                           cAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTGTGGCATCAGCACCATGTCGGAcc                                                              >  2:4845474/1‑92 (MQ=255)

                           cAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACt                                                              <  2:4041576/93‑1 (MQ=255)

                            aaGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTc                                                             <  2:5447340/93‑1 (MQ=255)

                            aaGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTc                                                             >  1:2434193/1‑93 (MQ=255)

                             aGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCg                                                            >  2:3136650/1‑93 (MQ=255)

                              ggATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGt                                                           <  1:1393661/93‑1 (MQ=255)

                                  gTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAAt                                                       >  2:3310814/1‑93 (MQ=255)

                                    cGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGc                                                     <  2:3413584/93‑1 (MQ=255)

                                    cGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGc                                                     <  2:909412/93‑1 (MQ=255)

                                     gCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGcc                                                    <  2:2055536/93‑1 (MQ=255)

                                     gCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGcc                                                    >  2:3028846/1‑93 (MQ=255)

                                         cGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTc                                                <  2:431176/93‑1 (MQ=255)

                                          ggATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTcc                                               >  1:5600405/1‑93 (MQ=255)

                                          ggATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTcc                                               <  1:4482869/93‑1 (MQ=255)

                                          ggATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTcc                                               <  2:1130985/93‑1 (MQ=255)

                                           gATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCa                                              <  2:3604007/93‑1 (MQ=255)

                                              aaTCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgcc                                           <  2:4384443/93‑1 (MQ=255)

                                               aTCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgccg                                          <  2:1733021/93‑1 (MQ=255)

                                                tCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgccgc                                         <  1:1242410/93‑1 (MQ=255)

                                                tCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgccgc                                         >  1:6519816/1‑93 (MQ=255)

                                                tCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgccgc                                         >  2:2593868/1‑93 (MQ=255)

                                                 cAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgccgcc                                        <  2:2310028/93‑1 (MQ=255)

                                                 cAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgccgcc                                        <  1:3332484/93‑1 (MQ=255)

                                                  aGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCa                                       <  2:1264481/93‑1 (MQ=255)

                                                  aGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCa                                       >  2:5628145/1‑93 (MQ=255)

                                                    ggTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAAt                                     <  1:3951004/93‑1 (MQ=255)

                                                    ggTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAAt                                     >  1:6582575/1‑93 (MQ=255)

                                                     gTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATc                                    >  1:2236551/1‑93 (MQ=255)

                                                      tGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATcc                                   <  1:6380934/93‑2 (MQ=255)

                                                      tGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCg                                   <  1:4260714/93‑1 (MQ=255)

                                                       ggAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgc                                  >  2:1645819/1‑92 (MQ=255)

                                                        gAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgag                                 <  1:1943401/93‑1 (MQ=255)

                                                        gAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgag                                 <  2:1056330/93‑1 (MQ=255)

                                                         aGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgagg                                <  1:1700432/93‑1 (MQ=255)

                                                          gTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgagga                               <  1:4355788/93‑1 (MQ=255)

                                                             ggCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGcc                            <  2:548603/93‑1 (MQ=255)

                                                             ggCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGcc                            <  2:3884016/93‑1 (MQ=255)

                                                                gaagaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTc                         >  1:5394972/1‑93 (MQ=255)

                                                                gaagaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTc                         <  2:1486055/93‑1 (MQ=255)

                                                                gaagaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTc                         <  2:6355666/93‑1 (MQ=255)

                                                                 aagaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCa                        >  1:5854712/1‑93 (MQ=255)

                                                                 aagaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCa                        >  2:4080458/1‑93 (MQ=255)

                                                                  agaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCaa                       >  1:4937436/1‑93 (MQ=255)

                                                                   gaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCaac                      >  2:2486859/1‑92 (MQ=255)

                                                                   gaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCaaa                      >  1:4138269/1‑93 (MQ=255)

                                                                     aCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACa                    >  1:6159438/1‑93 (MQ=255)

                                                                     aCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACa                    >  1:1765585/1‑93 (MQ=255)

                                                                      cGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACAt                   <  1:3723619/93‑1 (MQ=255)

                                                                       gggTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACAtt                  >  2:5721057/1‑93 (MQ=255)

                                                                         gTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACAtttt                <  2:1925635/93‑1 (MQ=255)

                                                                         gTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACAtttt                <  1:4764703/93‑1 (MQ=255)

                                                                          tCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTg               >  2:5380389/1‑93 (MQ=255)

                                                                            aaTCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGcc             <  2:3791976/93‑1 (MQ=255)

                                                                            aaTCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGcc             <  2:2738463/93‑1 (MQ=255)

                                                                            aaTCACTGAGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCTAGTAGCCGTCTAACATTTTgtc             >  1:6465667/1‑91 (MQ=255)

                                                                             aTCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCct            >  1:3198572/1‑93 (MQ=255)

                                                                             aTCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCct            >  1:4197848/1‑93 (MQ=255)

                                                                             aTCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCct            <  1:1515390/93‑1 (MQ=255)

                                                                              tCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCctt           <  2:1395889/93‑1 (MQ=255)

                                                                                aCTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCcttct         >  1:4228642/1‑93 (MQ=255)

                                                                                aCTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCcttct         <  2:2131151/93‑1 (MQ=255)

                                                                                aCTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCcttct         <  2:5827358/93‑1 (MQ=255)

                                                                                 cTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCcttctt        <  2:5173670/93‑1 (MQ=255)

                                                                                  tGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCttcttc       <  1:273879/93‑1 (MQ=255)

                                                                                  tGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCttcttc       <  1:1064651/93‑1 (MQ=255)

                                                                                   gCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCg      <  2:2371144/93‑1 (MQ=255)

                                                                                     ggtgggtgCGTCTGGCATCAGCACCATGTCGGCCTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCgaa    <  1:1007956/93‑1 (MQ=255)

                                                                                     ggtgggtgCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCgaa    <  2:5317173/93‑1 (MQ=255)

                                                                                     ggtgggtgCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCgaa    <  2:2948328/93‑1 (MQ=255)

                                                                                      gTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCgaag   <  1:785124/93‑1 (MQ=255)

                                                                                       tGGGTGCGTCTGGCATCAGCGCCATGTAGGTCTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCgaaga  <  2:4352925/93‑1 (MQ=255)

                                                                                       tGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCgaaga  <  1:6334937/93‑1 (MQ=255)

GGTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGA  >  NC_000913/4055308‑4055487

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