BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	376,777	A→G	N272N (AAT→AAC)	yaiM ←	putative esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	376,777	0	A	G	99.7%	1232.7 / ‑5.2	385	N272N (AAT→AAC)	yaiM	putative esterase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (183/201); total (184/201)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

Only 100 of 392 total aligned reads displayed.

CAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTT  >  NC_000913/376687‑376866

cAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGtt                                                                                         <  1:5623896/93‑1 (MQ=255)

cAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGtt                                                                                         >  1:6608844/1‑93 (MQ=255)

 aGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTttgttg                                                                                        <  1:7737879/93‑1 (MQ=255)

     cTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATCTTCAGTTTGTTGGCGt                                                                                    <  2:7508588/93‑1 (MQ=255)

     cTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGt                                                                                    <  1:6130267/93‑1 (MQ=255)

     cTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGt                                                                                    <  1:7334209/93‑1 (MQ=255)

      tACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTg                                                                                   >  2:7920270/1‑93 (MQ=255)

      tACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTg                                                                                   >  1:4542641/1‑93 (MQ=255)

      tACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTg                                                                                   >  1:3985615/1‑93 (MQ=255)

      tACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTg                                                                                   >  1:2494001/1‑93 (MQ=255)

       aCGTCTTACGGGAGAAGTGTAGGGCGGATAAGGCGTTGACGCCGCATCCGGCAGTGGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTgg                                                                                  <  1:4342874/93‑1 (MQ=25)

         ggCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTa                                                                                >  2:6538952/1‑93 (MQ=255)

          gCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAg                                                                               >  1:3871163/1‑93 (MQ=255)

            ttACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGc                                                                             <  1:665403/93‑1 (MQ=255)

            ttACGGGCGAAGGGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGc                                                                             >  2:6369867/1‑93 (MQ=255)

                gggCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAata                                                                         <  1:6011085/93‑1 (MQ=255)

                gggCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAata                                                                         <  1:472583/93‑1 (MQ=255)

                gggCGAAGTGTAGCCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAata                                                                         <  2:6920139/93‑1 (MQ=255)

                 ggCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAatat                                                                        <  1:5320732/93‑1 (MQ=255)

                 ggCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAatat                                                                        >  2:6225974/1‑93 (MQ=255)

                  gCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATg                                                                       >  1:7514270/1‑93 (MQ=255)

                    gAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCt                                                                     <  1:2131700/93‑1 (MQ=255)

                     aaGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTc                                                                    >  2:700590/1‑93 (MQ=255)

                      aGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCg                                                                   >  2:307485/1‑93 (MQ=255)

                          tAGGCCGGATAAGGGGTTCACGCCGACTGCTGCAGTTGTTTTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCaa                                                               <  2:6444471/93‑1 (MQ=255)

                          tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCaa                                                               <  1:5789569/93‑1 (MQ=255)

                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATa                                                             <  2:2889061/93‑1 (MQ=255)

                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATa                                                             <  2:7280102/93‑1 (MQ=255)

                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATa                                                             <  2:278226/93‑1 (MQ=255)

                               cGGATAAGTCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATaaaa                                                          <  2:2433866/93‑1 (MQ=255)

                               cGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATaaaa                                                          <  1:7926353/93‑1 (MQ=255)

                               cGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATaaaa                                                          >  1:4855267/1‑93 (MQ=255)

                                  aTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTg                                                       <  1:1527170/93‑1 (MQ=255)

                                  aTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTg                                                       <  2:6797473/93‑1 (MQ=255)

                                   tAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTgg                                                      >  2:3977433/1‑93 (MQ=255)

                                   tAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTgg                                                      <  2:1203843/93‑1 (MQ=255)

                                     aGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGgag                                                    >  1:5674430/1‑93 (MQ=255)

                                     aGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGgag                                                    >  2:1038575/1‑93 (MQ=255)

                                          ttCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACaaa                                               >  2:6303313/1‑93 (MQ=255)

                                          ttCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACaaa                                               >  1:1717483/1‑93 (MQ=255)

                                          ttCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACaaa                                               <  1:3530143/93‑1 (MQ=255)

                                          ttCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACaaa                                               >  2:3322249/1‑93 (MQ=255)

                                          gtCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACaaa                                               >  2:4792600/2‑93 (MQ=255)

                                           tCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACaaaa                                              <  2:2119835/93‑1 (MQ=255)

                                             aCGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAaata                                            >  1:3258649/1‑93 (MQ=255)

                                               gccgcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAaataat                                          >  2:4299323/1‑93 (MQ=255)

                                               gccgcATCCGGCAGTCGTGCTATAATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAaataat                                          <  2:5147429/93‑1 (MQ=255)

                                                ccgcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAatagta                                         <  2:3527853/93‑1 (MQ=255)

                                                ccgcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAataata                                         <  1:2600234/93‑1 (MQ=255)

                                                ccgcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAataata                                         >  2:7364240/1‑93 (MQ=255)

                                                  gcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGc                                       <  2:3152146/93‑1 (MQ=255)

                                                  gcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGc                                       <  2:377236/93‑1 (MQ=255)

                                                  gcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGc                                       >  2:5382500/1‑93 (MQ=255)

                                                  gcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGc                                       <  2:9366/93‑1 (MQ=255)

                                                  gcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGc                                       >  1:7628052/1‑93 (MQ=255)

                                                  gcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGc                                       >  1:7063486/1‑93 (MQ=255)

                                                  gcATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGc                                       <  1:3720078/93‑1 (MQ=255)

                                                     tCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCtgt                                    <  1:6900690/93‑1 (MQ=255)

                                                       cGGCAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGa                                  >  2:7782120/1‑93 (MQ=255)

                                                          cAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCa                               <  1:1578271/93‑1 (MQ=255)

                                                          cAGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCa                               >  1:2321170/1‑93 (MQ=255)

                                                           aGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCAt                              <  2:554811/93‑1 (MQ=255)

                                                           aGTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCAt                              >  1:5850740/1‑93 (MQ=255)

                                                            gTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATa                             >  1:6471199/1‑93 (MQ=255)

                                                            gTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATa                             >  1:6330363/1‑93 (MQ=255)

                                                            gTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATa                             >  1:7259164/1‑93 (MQ=255)

                                                            gTCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAACAGCTGTGATCATa                             >  1:6283614/1‑93 (MQ=255)

                                                            gTCGTGCTATTATCAACGCATATTCAGTTTGATGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATa                             >  1:3111973/1‑93 (MQ=255)

                                                             tCGTGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATaa                            <  2:5966418/93‑1 (MQ=255)

                                                                tGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAccc                         <  2:7849307/93‑1 (MQ=255)

                                                                tGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAccc                         >  1:1459738/1‑93 (MQ=255)

                                                                tGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAccc                         >  2:5584377/1‑93 (MQ=255)

                                                                tGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAccc                         <  2:1845618/93‑1 (MQ=255)

                                                                tGCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAAccc                         >  1:5034818/1‑93 (MQ=255)

                                                                 gCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCt                        <  2:4352885/93‑1 (MQ=255)

                                                                 gCTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCt                        >  1:4947374/1‑93 (MQ=255)

                                                                  cTATTATCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTc                       <  1:5107717/93‑1 (MQ=255)

                                                                   tattatCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCg                      >  1:6404243/1‑93 (MQ=255)

                                                                   tattatCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCg                      >  1:542887/1‑93 (MQ=255)

                                                                    attatCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGt                     <  2:2194286/93‑1 (MQ=255)

                                                                     ttatCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTg                    >  2:6120696/1‑93 (MQ=255)

                                                                     ttatCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTg                    <  2:5668191/93‑1 (MQ=255)

                                                                      tatCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGa                   <  2:4610002/93‑1 (MQ=255)

                                                                        tCAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATa                 >  1:788810/1‑93 (MQ=255)

                                                                         cAACGCATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATaa                <  1:1317668/93‑1 (MQ=255)

                                                                             gCATATTCAGTTTGTTGGCGTGGGAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACgta            >  2:1314675/1‑91 (MQ=255)

                                                                              cATATTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGAt           <  1:6820159/93‑1 (MQ=255)

                                                                                tatTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTa         >  1:4162090/1‑93 (MQ=255)

                                                                                tatTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTa         >  1:3370860/1‑93 (MQ=255)

                                                                                tatTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTa         >  1:1163338/1‑93 (MQ=255)

                                                                                 atTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTaa        <  2:4844977/93‑1 (MQ=255)

                                                                                 atTCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTaa        <  1:2724493/93‑1 (MQ=255)

                                                                                  ttCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAAc       <  1:1017866/93‑1 (MQ=255)

                                                                                  ttCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAAc       <  2:414532/93‑1 (MQ=255)

                                                                                   tCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACg      <  2:5838539/93‑1 (MQ=255)

                                                                                   tCAGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACg      >  2:3404370/1‑93 (MQ=255)

                                                                                     aGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTc    >  1:1920724/1‑93 (MQ=255)

                                                                                     aGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTc    <  2:2627802/93‑1 (MQ=255)

                                                                                     aGTTTGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTc    >  1:7879881/1‑93 (MQ=255)

                                                                                       tttGTTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCtt  <  1:2789176/93‑1 (MQ=255)

CAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTT  >  NC_000913/376687‑376866

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`