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%	100.5 / 32.7	80	V288G (GTG→GGG)	cusA	copper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (2/41);  new base G (36/0);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 1.77e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599426‑599722                                                                                                                                                         |                                                                                                                                                      aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGggt                                                                                                                                                       <  1:154701/151‑1 (MQ=255)       tCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgg‑g‑ggggggggg                                                                                                                                                   >  1:1227614/1‑147 (MQ=255)          acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGTAGATGGGCCGGGGGATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCt                                                                                                                                              <  1:853534/151‑1 (MQ=255)           cgacTTTAATCTCACCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCTATGTTGCGAAGGTGCAGATTGGCCCGGAGATGCGCGGGGGCATTGCCGAACTAAACGGCGAAGACGAAGTGGCGGG‑C‑GGGGTGGTGCTcc                                                                                                                                               <  1:127137/149‑1 (MQ=255)             acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑Ggggggggg                                                                                                                                                   >  1:464815/1‑141 (MQ=255)                ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑t                                                                                                                                       <  1:867945/150‑1 (MQ=255)                 tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑G‑GGGGGGGGGGTTCCGC‑GA‑Tcc                                                                                                                                     >  1:619403/1‑151 (MQ=255)                 tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑G‑GGGGGGGGGAGCCttggtctc                                                                                                                                       >  2:1097743/1‑144 (MQ=255)                  aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑G‑GGGGGGGGGGTCCTTC‑GA‑ccc                                                                                                                                     >  1:206228/1‑150 (MQ=255)                   aTCAAATCGTTTTAAAAGCCAGTGAAAATGACGTGGCCGTTTATCTGCGCGATGTTGAGAAGATCCAGATTGGCCCGGAGATGCGCCGGGGCAATGCCGAACTAAACAGCGAAGGCGAAGTGGCGGG‑C‑Gggggggaggatcagctaac                                                                                                                                        >  1:205699/1‑135 (MQ=255)                     cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGc                                                                                                                                 <  1:1006870/150‑1 (MQ=255)                       caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGAGGCGGG‑C‑GGGGGGGGGGGCGTGC‑GG‑TCC‑GGGaa                                                                                                                               >  1:599404/1‑150 (MQ=255)                        aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCaaaa                                                                                                                             <  2:1150886/151‑1 (MQ=255)                         tCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGATCCTCG‑GG‑TTC‑GGCaaaa                                                                                                                             >  2:405778/1‑150 (MQ=255)                          cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCaaaaa                                                                                                                            <  1:240574/150‑1 (MQ=255)                          cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGTGTTCTGG‑GA‑GCC‑GGCaaaaa                                                                                                                            >  2:917900/1‑150 (MQ=255)                                aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGGTCCCGC‑GA‑TCC‑CGCAAAAAAGccc                                                                                                                       >  1:236054/1‑149 (MQ=255)                                aaaaGCCAGTGAAAATGGCGTGACCGTTTATCTGCGCGATGTTGCGAACGACCAGATTGGCCCGGAGCTGCGCCGGGGCACTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGAC‑GGGGGGGAGTTCCGCC‑GC‑TCCGGGCAAAAAccc                                                                                                                         >  1:1292443/1‑147 (MQ=255)                                 aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑CGGGGGGGGGGACCCTCG‑GA‑CCC‑GCCAAAAACGcc                                                                                                                        >  1:324039/1‑148 (MQ=255)                                 aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCgag                                                                                                                    <  1:404804/151‑1 (MQ=255)                                  aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCgaga                                                                                                                   <  2:489702/151‑1 (MQ=255)                                   aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGaa                                                                                                                  <  2:1337366/151‑1 (MQ=255)                                    gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑G‑GGGGGGGGGGTCCTGCGGA‑TCC‑GCAAAAAACGccgag                                                                                                                     >  1:1117385/1‑145 (MQ=255)                                     ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGGTCTTGC‑GG‑TCC‑GGCgaaaaaacccccgga                                                                                                                  >  1:983856/1‑139 (MQ=255)                                        gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGAt                                                                                                             <  2:336001/151‑1 (MQ=255)                                          gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATc                                                                                                            <  2:4974/150‑1 (MQ=255)                                          gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATc                                                                                                            <  1:1292951/150‑1 (MQ=255)                                           aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑G‑GGGGGGGGGGGACCGC‑GC‑CCC‑GGCGCAAACGCCCCCGAAGGGGTcg                                                                                                           >  2:1007021/1‑150 (MQ=255)                                             aaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑G‑GGGGGGGGGGTCCCGC‑GATTCC‑GCCAAAAAACCCCgggaggggttgcc                                                                                                          >  2:574496/1‑141 (MQ=255)                                              aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGCGGGGAGCCCGG‑GA‑ACC‑GGCAAAAAACCCCgaggaggggcggcccc                                                                                                       >  2:772083/1‑141 (MQ=255)                                                ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑G‑GGGGGGGGGATGCTGA‑GA‑GCC‑GGCAAAAAAACCCGGGAAGGGAAcgcc                                                                                                         >  2:765785/1‑147 (MQ=255)                                                  cGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGTGATCCTGC‑GA‑TCC‑GGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTg                                                                                                   >  2:640915/1‑151 (MQ=255)                                                   gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGATCCTGC‑GG‑TCC‑GGCCAAAACGCCCGAGAAGGGATCGCCGCCGTg                                                                                                   >  2:217225/1‑150 (MQ=255)                                                     gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                 <  2:1037682/150‑1 (MQ=255)                                                     gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAg                                                                                                <  1:686195/151‑1 (MQ=255)                                                       ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGa                                                                                              <  2:965119/151‑1 (MQ=255)                                                        cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc                                                                                             <  1:693597/151‑1 (MQ=255)                                                            tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGGTCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGGGATCGCCGCCGTggaggaacaa                                                                                          >  2:434741/1‑146 (MQ=255)                                                             aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGGACCTGC‑GG‑TTC‑GGCGAAAAAGCCCGGGGAGAGATTGCCGCCGGGAAGGACaaa                                                                                          >  2:1373020/1‑149 (MQ=255)                                                             aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                          >  2:347897/1‑149 (MQ=255)                                                                   gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                  <  2:569184/151‑1 (MQ=255)                                                                   gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                  <  2:1358093/151‑1 (MQ=255)                                                                     gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGTGATTCTGC‑GA‑TCC‑GGGAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGGCAAACTGGaaa                                                                                   >  2:388944/1‑148 (MQ=255)                                                                     gATGTTGCGAAGCTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACg                                                                                 <  1:616138/150‑1 (MQ=255)                                                                      aTGTTGCGAAGGTCCAGGTTGGCCCGGAGAGGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCATTC‑TA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                               <  1:167851/151‑1 (MQ=255)                                                                              aaGGTCCAGATTGGCCCGGAGAGGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAag                                                                         <  1:1256990/149‑1 (MQ=255)                                                                               aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                      <  2:537570/151‑1 (MQ=255)                                                                               aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCCAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaaa                                                                          >  2:545695/1‑147 (MQ=255)                                                                                 gtacAGATTAGGCGGGCGAGGCTCCGGGGCATTGCCAAATTAACCGGAGAAGGCGAAGTGGCGGG‑T‑GGGTTGGTGTTCTTGA‑AA‑TCC‑AACAAAATCCCCCAAGAAGTGATCCCCGCCGTGAAGGACAAACTGGAAGCGCTGAAAagtag                                                                      <  1:969993/146‑1 (MQ=255)                                                                                  tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑G‑GGGGGGGGGGTACCGC‑GA‑GCC‑CGCCAAAACCCCCCAGAAAGGGTTGGCGCCGGGGAGGGCAAACCGGAAACGCCGAAAagtag                                                                      >  1:741990/1‑148 (MQ=255)                                                                                  tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGGTCCTGC‑GA‑TCC‑GGCAAAAACGCCCGCGAAGTGATCGCCGCCGTGAAGGGCAAACCGGAAACGCGGAAAagtagatt                                                                   >  2:608002/1‑148 (MQ=255)                                                                                   ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑G‑GGGGGGGTGGACCTGG‑GA‑TTC‑CGCAAAAAAGCCCGAGGAGGGAACGCCGCCGTGGAGGGCAAAATGGAAAAGCTGAAAAGGAGGCt                                                                   >  2:709301/1‑150 (MQ=255)                                                                                    cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGc                                                                 <  2:1350153/151‑1 (MQ=255)                                                                                       aTTGGCCCGGAGATGCGCCGTGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCg                                                               <  2:249/150‑1 (MQ=255)                                                                                       aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGATTCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGGAAAACTTGGAAAGCTGAAAAATAGTCTGCCgg                                                              >  2:557144/1‑151 (MQ=255)                                                                                          ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa                                                            <  1:107496/150‑1 (MQ=255)                                                                                           gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                          >  2:42604/1‑151 (MQ=255)                                                                                              cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGt                                                       <  2:861485/151‑1 (MQ=255)                                                                                                gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGAAACTGC‑GA‑ACC‑AACAAAAACCGCCCCGAAGTGATTGCCGCCGCGGAGGAAAAAATGGAACAGCTGGAAAGTAGTCCGCCGGGAAGCgggg                                                     >  2:774891/1‑151 (MQ=255)                                                                                                   aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGGGAAGGGATCGCCGCCGTGAAGGGCAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCgggggg                                                   >  2:170611/1‑148 (MQ=255)                                                                                                     gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                 <  1:1123461/150‑1 (MQ=255)                                                                                                       gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGGAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAg                                               >  2:1197805/1‑150 (MQ=255)                                                                                                         cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTa                                            <  2:1147517/151‑1 (MQ=255)                                                                                                          ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGTGATCCTGC‑GG‑TCC‑GGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAACGCCGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTa                                            >  2:1009303/1‑150 (MQ=255)                                                                                                            ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGGTACCGC‑GA‑TCC‑GGCAAAAAAGCCCGAGAAGGGATCCCCGCCGTGAAGGACAAACCGGAAACGCTGAAAAGTAGCCTGCCGGAAGGGGCGGGGATAGtt                                             >  2:697721/1‑147 (MQ=255)                                                                                                              cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGATCCTGC‑GA‑TCC‑GGCCAAAAAGCCCGAGAAGGGATCGGCGCCGGGAAGGACAAACTGGGAAAGCTGAAAAGTAGTCTGGCGGGAGGCGGGGGGATAGTTacaa                                         >  1:1146083/1‑149 (MQ=255)                                                                                                               aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaaca                                       <  2:518475/150‑1 (MQ=255)                                                                                                               aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                      <  1:976629/151‑1 (MQ=255)                                                                                                               aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCCCCGCGAAGGAAAAAAACGAAACGCTGAAAAGAAGACTGCCGGAAGGGGGGGGGGTAGTTaaaaaa                                       >  2:13846/1‑148 (MQ=255)                                                                                                               aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGGGATTCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGAGATCGCCGCCGAGAAGGACAAACTGGAAACGCGGCAAAGAAATCTGCCGGAAGGAGGGGAGATAGATacaaaa                                       >  2:541322/1‑148 (MQ=255)                                                                                                                ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                      <  1:480466/150‑1 (MQ=255)                                                                                                                  gCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                   <  2:40050/151‑1 (MQ=255)                                                                                                                   ccGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTTGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGa                                  <  1:607971/151‑1 (MQ=255)                                                                                                                   ccGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGATATAGTTACAACATACGa                                  <  1:209883/151‑1 (MQ=255)                                                                                                                     gAACTAAACGGCGAAGGCGAAGTGGCGGG‑G‑GGGGGGGGGATCCTGG‑CA‑TCC‑GGCAAAAAAGCCCGAGAAGGGATTGCCGCCGGGGAGGACAAACTTGAAAACGTGAAAAGGAGGCTGCCGGAAGGGGTGGAGATAGATACAAAAAACGa                                  >  2:972280/1‑149 (MQ=255)                                                                                                                     gAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                <  2:1146623/151‑1 (MQ=255)                                                                                                                      aaCTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGCAGGCGGGGAGATAGTTACAACATACGATc                                >  2:788502/1‑150 (MQ=255)                                                                                                                       aCTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGGGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATACGATCGc                              >  2:1281917/1‑151 (MQ=255)                                                                                                                             cGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagccag                        <  2:619434/151‑1 (MQ=255)                                                                                                                              ggaaaaggagaatTGGCGGG‑C‑GTGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAACCGCCCGAGAAGTGTTTTCAGCCGTGAAGGACAAACTGGAAACGCTGAACAGTAGTCTGCCGTAAGGCGTGTACATAGTTACCACATACAATCGcagccag                        <  1:1221131/146‑1 (MQ=255)                                                                                                                                 gaaggcgaagTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCa                    <  2:451379/151‑1 (MQ=255)                                                                                                                                       gaagTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc              <  2:611837/151‑1 (MQ=255)                                                                                                                                            ggcggg‑c‑gGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGcc         >  2:228316/1‑151 (MQ=255)                                                                                                                                                     ggggTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa  <  2:286706/151‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                      AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGG‑C‑GGGGTGGTGATCCTGC‑GA‑TCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599426‑599722

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