Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	514,292	C→T	6.1%	S100S (TCC→TCT)	cueR →	copper‑responsive regulon transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	514,292	0	C	T	6.1%	207.9 / 2.6	65	S100S (TCC→TCT)	cueR	copper‑responsive regulon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base T (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTT  >  NC_000913/514146‑514438                                                                                                                                                   |                                                                                                                                                   cAGGCACGGCAGGTGGGCGTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCAt                                                                                                                                                  <  7:317622/149‑1 (MQ=255)   ggCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATgc                                                                                                                                                <  5:269456/149‑1 (MQ=255)   ggCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATgc                                                                                                                                                <  6:78342/149‑1 (MQ=255)                ggCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAgctgct                                                                                                                                   <  3:305644/149‑1 (MQ=255)                   tttAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTgg                                                                                                                                 >  2:46154/1‑148 (MQ=255)                   tttAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTgg                                                                                                                                 >  2:62897/1‑148 (MQ=255)                       aCCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACt                                                                                                                             <  1:229101/148‑1 (MQ=255)                        ccTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTgg                                                                                                                           >  7:304005/1‑149 (MQ=255)                          tGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCTATGCGCGACCAGCTGCTGGCACTGGc                                                                                                                          >  6:109389/1‑148 (MQ=255)                          tGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCTATGCGCGACCAGCTGCTGGCACTGGc                                                                                                                          >  8:18566/1‑148 (MQ=255)                             aaGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGATGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAt                                                                                                                      <  3:113475/149‑1 (MQ=255)                             aaGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAt                                                                                                                      <  5:301027/149‑1 (MQ=255)                              agagagCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtg                                                                                                                     >  1:238449/1‑149 (MQ=255)                               gagagCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGTCACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgc                                                                                                                    <  2:310657/149‑1 (MQ=255)                               gagagCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAAtgc                                                                                                                    >  2:264624/1‑149 (MQ=255)                                     ggcgAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGccc                                                                                                              <  2:82568/149‑1 (MQ=255)                                      gcgAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCt                                                                                                             >  8:308430/1‑149 (MQ=255)                                            tGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGAt                                                                                                       >  7:185313/1‑149 (MQ=255)                                             ggTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGAt                                                                                                       >  7:159141/1‑148 (MQ=255)                                               tGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGa                                                                                                     >  5:52899/1‑148 (MQ=255)                                                   tCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAgcg                                                                                                >  5:359679/1‑149 (MQ=255)                                                     tGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGcc                                                                                              >  5:310163/1‑149 (MQ=255)                                                       tttAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGcc                                                                                              <  4:227648/147‑1 (MQ=255)                                                         tAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGAc                                                                                           <  6:52899/148‑1 (MQ=255)                                                          aaCGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACt                                                                                          <  8:185313/148‑1 (MQ=255)                                                          aaCGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACt                                                                                          <  1:62897/148‑1 (MQ=255)                                                          aaCGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTg                                                                                         >  5:325880/1‑149 (MQ=255)                                                            cGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGcc                                                                                       >  7:85567/1‑149 (MQ=255)                                                                  gCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTa                                                                                 >  2:203269/1‑149 (MQ=255)                                                                   cagagtcacaGCGACGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTAt                                                                                <  6:121184/143‑1 (MQ=255)                                                                    aGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATc                                                                               <  6:156236/149‑1 (MQ=255)                                                                        gCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGaaa                                                                           <  7:308430/149‑1 (MQ=255)                                                                           caGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAAtc                                                                        <  6:359679/149‑1 (MQ=255)                                                                                   aCGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGCTCGAACGCCCCATTGCGGAGCTGCCATCCCTGCGCGCCCAGCTGCTGGCACGGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGAGTATCGAAAATCTCTCCGgc                                                                >  3:182299/1‑149 (MQ=255)                                                                                      tCACACGGCGCACGCTGGAGAAGGGGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGAACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGgctgc                                                             <  2:104413/149‑1 (MQ=255)                                                                                        aaaCGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGctgctg                                                           >  2:52835/1‑149 (MQ=255)                                                                                              cgcACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGtcatca                                                     >  2:222163/1‑149 (MQ=255)                                                                                                                gCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGAt                                                                                                       <  4:267496/81‑1 (MQ=255)                                                                                                                gCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGAt                                                                                                       >  3:267496/1‑81 (MQ=255)                                                                                                                 cGGAGATCGAACGACACATTGAGGAGCTGCAATCTATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAaga                                  <  5:109389/149‑1 (MQ=255)                                                                                                                 cGGAGATCGAACGACACATTGAGGAGCTGCAATCTATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAaga                                  <  7:18566/149‑1 (MQ=255)                                                                                                                 cGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAaga                                  <  1:203269/149‑1 (MQ=255)                                                                                                                 cGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAaga                                  <  1:46154/149‑1 (MQ=255)                                                                                                                    agaTCGAACGACACATTGAGGAGCGGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacga                               <  8:159418/149‑1 (MQ=255)                                                                                                                      aTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGtcatca                                                     <  2:161622/125‑1 (MQ=255)                                                                                                                      aTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGtcatca                                                     >  1:161622/1‑125 (MQ=255)                                                                                                                      aTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacga                               >  3:337976/1‑147 (MQ=255)                                                                                                                      aTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacga                               <  4:337976/147‑1 (MQ=255)                                                                                                                       tCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagacgaga                             <  8:85567/148‑1 (MQ=255)                                                                                                                         gAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACgg                          >  5:339708/1‑149 (MQ=255)                                                                                                                         gAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACgg                          >  5:31530/1‑149 (MQ=255)                                                                                                                          aaCGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGc                         >  4:56991/1‑149 (MQ=255)                                                                                                                          aaCGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGc                         <  3:56991/149‑1 (MQ=255)                                                                                                                           aCGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCt                        >  2:224817/1‑149 (MQ=255)                                                                                                                             gACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACTGCTCg                      >  3:289108/1‑149 (MQ=255)                                                                                                                                acaTTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGa                    >  5:167317/1‑148 (MQ=255)                                                                                                                                acaTTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGa                    >  8:137880/1‑148 (MQ=255)                                                                                                                                  aTTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagac                                 >  4:113799/1‑133 (MQ=255)                                                                                                                                  aTTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAagac                                 <  3:113799/133‑1 (MQ=255)                                                                                                                                      aggagCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAg             <  7:137880/149‑1 (MQ=255)                                                                                                                                        gagCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAggg           <  8:304005/149‑1 (MQ=255)                                                                                                                                          tcggCAAGCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAACTCTCTACGGCTGCTGTCATCAGCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGtt         <  4:182299/146‑1 (MQ=255)                                                                                                                                              cAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGa                                                                                    <  7:335556/70‑1 (MQ=255)                                                                                                                                              cAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGa                                                                                    >  8:335556/1‑70 (MQ=255)                                                                                                                                                 tCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCtt  >  8:345527/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CAGGCACGGCAGGTGGGCTTTAACCTGGAAGAGAGCGGCGAGCTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTT  >  NC_000913/514146‑514438

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑