Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	599,579	0	T	G	39.4%	106.3 / 11.5	66	V289G (GTG→GGG)	cusA	copper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (10/30);  new base G (26/0);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.57e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCA  >  NC_000913/599434‑599728                                                                                                                                                   |                                                                                                                                                          gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTcc                                                                                                                                                      >  2:1048832/1‑151 (MQ=255)  acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATcc                                                                                                                                                      >  2:671857/1‑150 (MQ=255)   cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATcc                                                                                                                                                      >  2:511655/1‑149 (MQ=255)     acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGGCCTGC‑g                                                                                                                                                 >  2:1124236/1‑151 (MQ=255)      cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGGTcc                                                                                                                                                      >  2:201653/1‑146 (MQ=255)          aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGGTCCTGC‑GGTc‑c                                                                                                                                            >  2:61754/1‑150 (MQ=255)           aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑Cgg                                                                                                                                          <  1:920820/151‑1 (MQ=255)            tCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGGGTCTGC‑GATc‑ccg                                                                                                                                          >  1:363612/1‑148 (MQ=255)               caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGGTCCTGC‑GATC‑Cggga                                                                                                                                        >  2:510390/1‑147 (MQ=255)                aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGATTCTGC‑GGTT‑CGGCaa                                                                                                                                       >  1:1166246/1‑149 (MQ=255)                aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGGGCTGG‑GCTC‑CGGCaaa                                                                                                                                      >  1:662185/1‑150 (MQ=255)                  cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGC‑GATC‑CCGCaaaaa                                                                                                                                    >  2:76679/1‑150 (MQ=255)                   gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGGGTCCGC‑GAGC‑CGGCaaaa                                                                                                                                     >  1:172301/1‑148 (MQ=255)                    ttttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGGTCCTGC‑GATC‑CGGCaaaaa                                                                                                                                    >  2:417672/1‑148 (MQ=255)                     tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGc                                                                                                                                 <  1:1062365/150‑1 (MQ=255)                     tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGC‑GATC‑CGGCaaaaa                                                                                                                                    >  2:315446/1‑147 (MQ=255)                     tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGAATTGCCGAACTAAACGGCGAAGGCGAAGGGGC‑GGGCGGGGGGGTGGTTCTGC‑GATG‑CGGCacaaaaacc                                                                                                                                >  1:718036/1‑147 (MQ=255)                       tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGGGTCCGC‑GAGC‑CCGCCAAAACAccc                                                                                                                               >  1:492148/1‑150 (MQ=255)                        aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCg                                                                                                                              <  2:957731/150‑1 (MQ=255)                        aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGACGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGcc                                                                                                                                <  2:1103229/148‑1 (MQ=255)                         aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGC‑GATC‑CCGCAAAAAAGccc                                                                                                                               >  2:1227620/1‑148 (MQ=255)                             ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGt                                                                                                                        <  2:297626/151‑1 (MQ=255)                              cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTg                                                                                                                       <  1:402419/151‑1 (MQ=255)                               aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCCGGGGGGGGATCCTGC‑GATT‑CCGCCAAAACGCCCGCGAAGTg                                                                                                                       >  2:1223366/1‑150 (MQ=255)                                gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGAt                                                                                                                     <  1:152019/151‑1 (MQ=255)                                    aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATcgcc                                                                                                                 <  2:303561/151‑1 (MQ=255)                                       tGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGC‑GATC‑CGGCAAAAAAGCCCGAGAAGGGATCgccgcc                                                                                                              >  1:547155/1‑151 (MQ=255)                                        ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGATCCCGC‑GAGC‑CCGCAAACACGCCCGAGAAGAGATTgccgcc                                                                                                              >  2:953321/1‑150 (MQ=255)                                             gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                         <  2:448698/150‑1 (MQ=255)                                             gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                         <  2:239045/150‑1 (MQ=255)                                                cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc                                                                                                     <  2:481803/151‑1 (MQ=255)                                                  tttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑a                                                                                                   <  1:854114/150‑1 (MQ=255)                                                   ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGAC‑a                                                                                                   >  2:1008830/1‑149 (MQ=255)                                                   ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGC‑GCTC‑CCGCAAAAAAGCCCGAGAAGGGATTGCCGCCGGGGAGGAC‑aa                                                                                                  >  2:747251/1‑150 (MQ=255)                                                        tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGGACCCTC‑GTTCTCGGCAAAAA‑ACCCGCGGGGGGGTTGCCGCCGTGAGGGAG‑AAAA‑Tg                                                                                             >  1:703220/1‑149 (MQ=255)                                                        tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGATCCTGC‑GGTC‑CCGGAAAAACGCCCGAGAAGGGATCGCCGCCCGGAAGGAC‑AAAA‑Tgg                                                                                            >  2:95063/1‑150 (MQ=255)                                                           gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGGTTC‑CGGCAAAAA‑GCCCGCGGAGGGGTCGCCGCCGGGGAGGAC‑ACAC‑CGGaaa                                                                                         >  1:664863/1‑150 (MQ=255)                                                            cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGaaa                                                                                         >  2:121625/1‑149 (MQ=255)                                                                   gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGA‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGaaaa                                                                                <  2:703250/151‑1 (MQ=255)                                                                    cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGTCCCG‑‑GCGC‑CCGCAAAAACAACCGAGAAGAGATGGGCGCCCTGGAGGACAAAAC‑TGGaaaaatggaaa                                                                                 >  1:449360/1‑142 (MQ=255)                                                                      aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAagta                                                                             <  2:493241/151‑1 (MQ=255)                                                                                    cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                               <  1:1102855/151‑1 (MQ=255)                                                                                      cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGt                                                             <  1:65300/151‑1 (MQ=255)                                                                                      cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGt                                                             <  2:1092348/151‑1 (MQ=255)                                                                                        gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGC‑GATC‑CGGCAAAAAAACCCGAGGAGTGATTGCCGGCGCGAAAGAC‑AAAC‑TGGGAAAGCTGAAAAAAAGGCTGCCGGGAggggggg                                                           >  2:204770/1‑146 (MQ=255)                                                                                             gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                      <  1:1087848/151‑1 (MQ=255)                                                                                               gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGt                                                    <  1:1058632/151‑1 (MQ=255)                                                                                                  ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGATCCTGC‑GAAC‑CGGCCAAAACGCCCGGGAAGGGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGAGGTAGTTac                                                 >  2:1218226/1‑151 (MQ=255)                                                                                                    ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                               <  2:378806/151‑1 (MQ=255)                                                                                                      cATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAGTTacaa                                               >  2:1197666/1‑149 (MQ=255)                                                                                                        ttGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGGGGGGATAGTTacaaca                                             >  2:381653/1‑149 (MQ=255)                                                                                                         tGCCGAACTAAACGGGGAAGGCTAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGTCAAAACCGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATa                                           <  1:763049/150‑1 (MQ=255)                                                                                                         tGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATa                                           <  2:127461/150‑1 (MQ=255)                                                                                                              aaCTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGC‑GATC‑CGGCGAAAACGCCCCAGAAGTGATCGCCGCCGTGAAGGAC‑AAAA‑TGGAAACGCTGGAAAGGAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCg                                     >  2:670601/1‑151 (MQ=255)                                                                                                               aCTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGc                                    <  1:60758/151‑1 (MQ=255)                                                                                                                  aaaCGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                                   <  1:417659/149‑1 (MQ=255)                                                                                                                  aaaCGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                                   <  2:1166310/149‑1 (MQ=255)                                                                                                                   aaCGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGATCCTGC‑GATC‑CGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGGC‑AAAA‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAACGcagcc                                >  2:583728/1‑151 (MQ=255)                                                                                                                    aCGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGGGTTCCCGC‑GAGC‑CGGGAAAAACACCCGAGAAGGGATCGCCCCCCGGAAGGAA‑AAACATGGAAA‑ACTGAAAAGTATTTTGCCGCAAGGGCGGGGGAGAGTTACAACAAACGAGCGcagcc                                >  1:1253883/1‑150 (MQ=255)                                                                                                                         gaaggcgaagTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCa                          <  1:953269/151‑1 (MQ=255)                                                                                                                                  gTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgcg                 <  1:370798/151‑1 (MQ=255)                                                                                                                                      c‑gggcgggGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt             <  1:379651/151‑1 (MQ=255)                                                                                                                                        gggcgggGTGGTGATCCTGC‑GATC‑CGGCAAACACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc            >  1:1183447/1‑151 (MQ=255)                                                                                                                                        gggcgggGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc            <  1:1118768/151‑1 (MQ=255)                                                                                                                                        gggcgggGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc            <  1:46645/151‑1 (MQ=255)                                                                                                                                            ggggTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc         >  2:45019/1‑150 (MQ=255)                                                                                                                                                tggtgATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCt    >  2:528533/1‑151 (MQ=255)                                                                                                                                                  gtgATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCa  <  1:428213/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                          GACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGC‑GATC‑CGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAC‑AAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCA  >  NC_000913/599434‑599728

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