Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	599,576	0	T	G	45.6%	79.3 / 22.3	57	V288G (GTG→GGG)	cusA	copper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (7/24);  new base G (26/0);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 3.60e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCT  >  NC_000913/599434‑599726                                                                                                                                               |                                                                                                                                                        gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTcc                                                                                                                                                 >  1:788088/1‑151 (MQ=255) gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATcc                                                                                                                                                 >  1:594701/1‑151 (MQ=255)    gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGc                                                                                                                                              <  2:457047/151‑1 (MQ=255)      cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGggggtggggccgc                                                                                                                                              >  2:800096/1‑142 (MQ=255)                aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGCaaa                                                                                                                                   >  1:49472/1‑150 (MQ=255)                  cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaaaaa                                                                                                                                 <  2:604412/150‑1 (MQ=255)                    ttttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGCACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGGCTGGGATCCGGCgaaaaag                                                                                                                               >  2:942935/1‑148 (MQ=255)                      ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTTCGGCAAAAAAGccc                                                                                                                            >  1:732960/1‑151 (MQ=255)                       tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGccc                                                                                                                            >  2:105691/1‑150 (MQ=255)                            gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTGCTGCGATCCGGCCAAAAAGccc                                                                                                                            >  2:575466/1‑145 (MQ=255)                             ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCaaaaaaaccc                                                                                                                            >  1:52633/1‑139 (MQ=255)                             ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGt                                                                                                                     <  2:282434/151‑1 (MQ=255)                              cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTg                                                                                                                    <  2:608210/151‑1 (MQ=255)                              cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTg                                                                                                                    <  2:219556/151‑1 (MQ=255)                                  gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                                 <  2:174677/150‑1 (MQ=255)                                  gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                                 <  2:406769/150‑1 (MQ=255)                                  gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATcg                                                                                                                >  1:206343/1‑151 (MQ=255)                                     aaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccg                                                                                                             <  1:921202/151‑1 (MQ=255)                                     aaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccg                                                                                                             <  1:99959/151‑1 (MQ=255)                                      aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGGCCCGCAAAAACGCCCCAGGAGGGATcggc                                                                                                              >  2:348375/1‑147 (MQ=255)                                            tGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa                                                                                                       <  2:871722/150‑1 (MQ=255)                                              cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                                    <  1:453636/151‑1 (MQ=255)                                               ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCGGGAGAGGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAg                                                                                                     <  2:296438/149‑1 (MQ=255)                                                 gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTTCGGTCCGGCAAAAAAACCCCCGGAGTGAACGGCGCCGAGCAgg                                                                                                    >  1:226185/1‑148 (MQ=255)                                                    tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGTTCCTGCGATCCGGCAAAAAACCCCGGGAAGGGGTCGCCGCCGTGAAGGACaaa                                                                                               >  2:137082/1‑150 (MQ=255)                                                    tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGACGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGTGATCCGGCACAAAAGGCCGGGGAGGGATCGGCGCCGTTAAgggaaa                                                                                                >  1:644072/1‑145 (MQ=255)                                                          cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑gg                                                                                          <  2:791772/148‑1 (MQ=255)                                                            cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGGAGGACAAACT‑GGaaa                                                                                       >  1:790137/1‑149 (MQ=255)                                                               tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCCAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGGCAAACT‑GGAAACGCTg                                                                                  >  1:187564/1‑151 (MQ=255)                                                               tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGGTCGGCGCCGGGGAGGACAAACT‑GGAAAccgt                                                                                   >  1:104909/1‑147 (MQ=255)                                                                gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACg                                                                                     <  1:498452/147‑1 (MQ=255)                                                                gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTg                                                                                  <  1:409452/150‑1 (MQ=255)                                                                     gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCAAAAACGCCCGCGGAGGGATTGCCGCCGTGGAGGAGAAAATGGGAAACCCGGaaa                                                                               >  2:689629/1‑149 (MQ=255)                                                                      aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCCAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGGCCAACT‑GGAAAAGCTGAAAagg                                                                            >  2:614463/1‑149 (MQ=255)                                                                         gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTc                                                                        <  1:307707/151‑1 (MQ=255)                                                                            cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACT‑GGAAAAGCTGAAAAGTAGTCTg                                                                      >  1:903295/1‑150 (MQ=255)                                                                                  ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGaa                                                                <  1:803122/150‑1 (MQ=255)                                                                                         agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTgg                                                         >  1:213040/1‑150 (MQ=255)                                                                                          gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                       <  2:282156/151‑1 (MQ=255)                                                                                            tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga                                                      <  1:970444/150‑1 (MQ=255)                                                                                                 cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCCAGAAGGGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGAAGTCTGCCGGAAGGCGTGGGGATAGtt                                                 >  1:182454/1‑150 (MQ=255)                                                                                                  ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTa                                                <  1:800101/150‑1 (MQ=255)                                                                                                      cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGGGATAGTTacaa                                             >  1:924575/1‑149 (MQ=255)                                                                                                       aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaaca                                           <  1:886622/150‑1 (MQ=255)                                                                                                        ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAAGACAAACT‑GGAAAAGCTGAAAAGTAGGCTGCCGGGAGGCGGGGGGGTAGGTaaaaaa                                           >  2:254551/1‑147 (MQ=255)                                                                                                         tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACT‑GGGAAAGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAGTTacaacaaa                                         >  1:280166/1‑148 (MQ=255)                                                                                                               aCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGc                                  <  1:256288/151‑1 (MQ=255)                                                                                                                   aaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcc                              <  1:335078/151‑1 (MQ=255)                                                                                                                    aCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCGAAAAAGCCCGGGGAGGGATCGCCGCCGGGGAGGGAAAAAC‑GGGAACGCTGGAAAAGAGGTTGCCGGGAGGGGGGGGGGGAAGTACAACATaacgccgcagcc                              >  2:368769/1‑139 (MQ=255)                                                                                                                         gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCa                        <  1:326959/151‑1 (MQ=255)                                                                                                                          aaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGGAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGATAGTTACAACATACGATCGCAGCCAGCTCa                        >  1:404751/1‑150 (MQ=255)                                                                                                                                  gTGGCGGGGGGGGGGGGGGTCCTGCGATCCGGCCAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACT‑GGAAACGCTGGAAAGGAGTCTGCCGGGAGGCGGGGGGATAGTTACAAAATACGATCGCAGCCAGCTCATTTAcc                  >  2:862502/1‑148 (MQ=255)                                                                                                                                      cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt           >  1:735052/1‑151 (MQ=255)                                                                                                                                      cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGAGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt           >  2:859753/1‑151 (MQ=255)                                                                                                                                         gcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCg         >  1:645961/1‑150 (MQ=255)                                                                                                                                          cgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc       >  1:784547/1‑151 (MQ=255)                                                                                                                                           ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc       >  2:956266/1‑150 (MQ=255)                                                                                                                                            gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa      >  1:454638/1‑150 (MQ=255)                                                                                                                                              gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc   <  2:292007/151‑1 (MQ=255)                                                                                                                                               tggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCt  >  1:111371/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                        GACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACT‑GGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCT  >  NC_000913/599434‑599726

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