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

ACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCA  >  NC_000913/599438‑599728                                                                                                                                               |                                                                                                                                                        acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCg                                                                                                                                                <  1:955680/151‑1 (MQ=255) acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCg                                                                                                                                                <  1:724377/151‑1 (MQ=255)   tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGa                                                                                                                                                       <  2:699468/142‑1 (MQ=255)   tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCg                                                                                                                                                >  1:1876/1‑149 (MQ=255)           caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGGGACCGCG‑ATcc                                                                                                                                           >  2:1188618/1‑145 (MQ=255)            aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGGCCTGCG‑GTCCCGCaaa                                                                                                                                     >  2:1051492/1‑150 (MQ=255)                 tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGCTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTTCTGCG‑GTCCCGCGAAAAAGcc                                                                                                                               >  2:1125281/1‑151 (MQ=255)                  ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTCCG‑GCCCGGCAAAAAc                                                                                                                                  >  2:21466/1‑147 (MQ=255)                       aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCG‑ATCCGGCAAAAACGCCCg                                                                                                                             >  1:397589/1‑147 (MQ=255)                         ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGt                                                                                                                       <  1:1091354/151‑1 (MQ=255)                         ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGt                                                                                                                       <  2:760721/151‑1 (MQ=255)                          cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGTGGCGGTCCTGCG‑GACTCGCCAAAAAGCCCGAGGAGTg                                                                                                                      >  1:494166/1‑151 (MQ=255)                           aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCG‑ATCCGGCAAAAAAGCCCGAGAAggg                                                                                                                      >  1:496351/1‑148 (MQ=255)                            gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGAt                                                                                                                    <  2:222616/151‑1 (MQ=255)                              gAAAATGGCGTGCCCGGTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                                   <  2:106039/150‑1 (MQ=255)                                aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTTCTGCG‑ATCCGGCAAAAACGCCCGGGAAGTGATcgc                                                                                                                 >  1:620731/1‑150 (MQ=255)                                  aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGAGCCCCCG‑GGCCCGCAAAAAAGACCGCGaaagggtcccc                                                                                                                >  2:1234014/1‑140 (MQ=255)                                   tGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTTCTCCG‑GTCCCGGAAAAACGCCCGGGGAGGGATTgccgcc                                                                                                             >  2:470396/1‑151 (MQ=255)                                   tGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGATCC‑GCG‑GTCCGGCAAAAACGCCCGAGAAGGGAGCgccgcc                                                                                                             >  2:865804/1‑151 (MQ=255)                                       gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCG‑GTGCCGCAAAAACACCCGAGAAGAGGTcgccg                                                                                                               >  2:688396/1‑145 (MQ=255)                                         gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                        <  2:416664/150‑1 (MQ=255)                                          cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                        <  1:423071/149‑1 (MQ=255)                                          cccGTTTATCGGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                                      <  2:772753/151‑1 (MQ=255)                                            cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTGCTTCGCAGCGGGCAAAAAC‑CCCGGGAAATGCCCCTCGCCGGGAAgg                                                                                                      >  2:895933/1‑149 (MQ=255)                                                tATCTGCGCGATGTCGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGGCCCTCG‑TGCCGGCAAAAAAGGCCCCGGAGGGATCgccgccccggagtgcaaa                                                                                                 >  2:81725/1‑138 (MQ=255)                                                 aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAc                                                                                                <  2:804859/150‑1 (MQ=255)                                                      cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGATCCTGCG‑ATCCCGCAAAAACGCCCGAGAAGTGATTGCCGCCGTGGAGGGCAAAA‑Tgg                                                                                            >  1:601198/1‑148 (MQ=255)                                                        cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCCGCG‑ATCCGGCCAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAAc                                                                                        >  1:1177137/1‑150 (MQ=255)                                                         gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCACTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCG‑ATCCGGCAAAAAAGCCCGAGAAGTGGTCGCCGCCGGGAAGGACCAAA‑TGGAAAAGc                                                                                      >  1:2483/1‑151 (MQ=255)                                                           tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCCGCG‑ATCCGGGCAAAACGCCCGAGAAGTGAGTGCCGCCGGGAAGGACAAAA‑TGGgaaacctg                                                                                    >  2:440457/1‑146 (MQ=255)                                                              tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCG‑ATCCGGCAAAAACGCCCGAGGAGGGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGaa                                                                                  >  1:919200/1‑150 (MQ=255)                                                                cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGGTCTGCG‑CACCGGCAAAAACCGCCGCGGAGTGATCGGCCACGGGAAGGGCAAAC‑TCGAAAAGCCGcaaa                                                                                >  1:687562/1‑150 (MQ=255)                                                                  aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCG‑GACCGGGAAAAAAGCCCGAGAAGGGATCGCCGCCGGGGAGGAGAAAC‑TGGaaa                                                                                         >  2:428353/1‑139 (MQ=255)                                                                   aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAagtag                                                                            <  2:1251972/151‑1 (MQ=255)                                                                     gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAagtag                                                                            >  1:160052/1‑149 (MQ=255)                                                                       ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGGCCCGCG‑GACCCGCCAAAACGCCCGAGAAAGGATCGCCGCCGTGGAGGGCCAACATGGAAAcccgaaaaataatc                                                                           >  2:385991/1‑136 (MQ=255)                                                                        cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCGGc                                                                       <  2:723300/151‑1 (MQ=255)                                                                            ttGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCgg                                                                    >  1:172987/1‑150 (MQ=255)                                                                             tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGGGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGa                                                                   <  2:566758/150‑1 (MQ=255)                                                                             tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCG‑ATCCCGCAAAAACGCCCGAGAAGGGATCGGCGCCGGGGAGGACAAAA‑TGGAAAAGCTGGAAAGGAGTCTGCCgg                                                                    >  1:551999/1‑149 (MQ=255)                                                                                      gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                         <  2:1038947/151‑1 (MQ=255)                                                                                        tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga                                                        <  1:948819/150‑1 (MQ=255)                                                                                              ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCCCCG‑ATCCCGCCAAAAACCCCCAGAAGGGGTCGCCCCCCTGGAGGGCAAAC‑TGGAAAAGCTGAAAAGTAGTCTGCCCGGAGGCGTGGAGATAGtt                                                   >  2:250471/1‑149 (MQ=255)                                                                                                 gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaac                                              <  2:208740/151‑1 (MQ=255)                                                                                                 gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaac                                              <  1:1209215/151‑1 (MQ=255)                                                                                                   aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                            <  1:331456/151‑1 (MQ=255)                                                                                                   aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                            <  2:142309/151‑1 (MQ=255)                                                                                                   aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGGGATCCTGCG‑ATCCGGCAAAAACGCCCGCGAAGTGATCGCCGCCCTGGAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGGGATAGTTacaaca                                             >  1:970545/1‑150 (MQ=255)                                                                                                      gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                         <  2:445364/151‑1 (MQ=255)                                                                                                      gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCG‑GTCCGGCAAAAAAGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAAA‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGGTAGTTACAACATa                                           >  1:129954/1‑149 (MQ=255)                                                                                                         gAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                       <  1:1028587/150‑1 (MQ=255)                                                                                                         gAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                      <  2:778444/151‑1 (MQ=255)                                                                                                          aaCTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCG‑ATCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGATAGTTACAACATACGATc                                      >  1:33991/1‑150 (MQ=255)                                                                                                               aaCGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATACTGCG‑ATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGGGGAGATAGTTACAACATACCATCGcagcc                                >  1:924395/1‑151 (MQ=255)                                                                                                                aCGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCG‑ATCCGGCAAAAAAGCCCGAGGAGTGATTGCCGCCGGGAAGGACAAAC‑TGGAAAAGCTGAAAAGGAGTCTGCCGGGAGGGGGGGGGGTAGTTACAAAAAAAGATCGcagcc                                >  2:1243017/1‑150 (MQ=255)                                                                                                                 cGGCGAAGGTGAAGTGGCGGGC‑GGGGGGGTGGTCCTGCG‑ATCCCGCAAAAAAGCCCGCGAAGGGATCGCCGCCGGGGAGGACAAAC‑CGGAAAAGCTGAAAAGTAGTCTGCCGGGAGGCGGGGGGAAAGTTACAACATACGATCGcagcc                                >  2:575685/1‑149 (MQ=255)                                                                                                                 cGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagccag                              <  1:381187/151‑1 (MQ=255)                                                                                                                      aaggcgaagTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAt                         <  2:405200/151‑1 (MQ=255)                                                                                                                       aggcgaagTGGCGGGC‑GGGGGGGGGATGCTTCG‑ATCCGGCAAAAAAGCCCGAGGAGTGAGCGCCGCCGTGTAGGACAAAC‑TGGAAACGCTGAAAAGAAGTCTGCCGGAAGGCGTGGGGAAAGTTACAACATAAGAACGCAGCCAGCCCAtt                        >  2:1084514/1‑151 (MQ=255)                                                                                                                           gaagTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc                    <  2:361518/151‑1 (MQ=255)                                                                                                                                   gggc‑gggGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc            >  2:265715/1‑151 (MQ=255)                                                                                                                                        ggggTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc         >  2:531447/1‑150 (MQ=255)                                                                                                                                              gtgATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTc   >  1:519240/1‑150 (MQ=255)                                                                                                                                              gtgATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCa  <  2:881042/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                        ACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCG‑ATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAC‑TGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCA  >  NC_000913/599438‑599728

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