Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	599,576	0	T	G	43.1%	85.0 / 26.3	59	V288G (GTG→GGG)	cusA	copper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (2/31);  new base G (25/0);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 2.01e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCT  >  NC_000913/599429‑599726                                                                                                                                                    |                                                                                                                                                        cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGggggggg                                                                                                                                                      >  2:576693/1‑150 (MQ=255)  gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg                                                                                                                                                     >  2:920283/1‑149 (MQ=255)         gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTg                                                                                                                                              <  1:1000854/150‑1 (MQ=255)         gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GGTCCTg                                                                                                                                              >  1:681060/1‑150 (MQ=255)           cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGG‑GGTcc                                                                                                                                                >  1:51858/1‑146 (MQ=255)               aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCg                                                                                                                                       <  1:822578/151‑1 (MQ=255)                       cattttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCAAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCaaaacc                                                                                                                               <  2:934456/149‑3 (MQ=255)                         ttttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGTCGAAGGCGAAGTGGCGGGCGGGGGGGT‑GATCCTGCGATCCGGCCAAAAAGc                                                                                                                             >  1:106566/1‑151 (MQ=255)                          tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGc                                                                                                                             <  2:538786/150‑1 (MQ=255)                            tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGccc                                                                                                                           <  1:44518/150‑1 (MQ=255)                                aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCTGCGAGCCGGCAAAAACGCCCgag                                                                                                                        >  2:506815/1‑149 (MQ=255)                                  ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGG‑GGTCCTGCGAGCCGGGCAAAAAGCCCGAGaa                                                                                                                      >  2:441272/1‑149 (MQ=255)                                    aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCTGCGGTCCGGCAAAAACGCCCGGGAAggg                                                                                                                   >  1:520020/1‑148 (MQ=255)                                      tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGT‑GGTCCTGCGATCCCGGCAAAACGACCGCGAAggg                                                                                                                   >  2:262542/1‑146 (MQ=255)                                        aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GGTCCCGCGATCCGGCAAAAAAGCCCGGGAAGGGGTcg                                                                                                               >  1:887398/1‑150 (MQ=255)                                         aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATcgcc                                                                                                             >  1:979307/1‑151 (MQ=255)                                           aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GGTCCTGCGGACCGGCAAAAAAGCCCGCGAAGGGGTcgccgc                                                                                                           >  2:917393/1‑151 (MQ=255)                                                gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTg                                                                                                       <  2:532716/150‑1 (MQ=255)                                                  gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGG‑GGTCCTGCGATCCGGCAAAAACGCCCGGGAAGGGATCGCCGCCGTg                                                                                                       >  2:395410/1‑148 (MQ=255)                                                  gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGT‑GATACTGCGGTCCGGCAAAAACGCCCGCGAAGGGGTCGCCGCCGTg                                                                                                       >  2:986183/1‑148 (MQ=255)                                                  gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGGTCGccgccggg                                                                                                       >  1:873492/1‑146 (MQ=255)                                                  gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCCTCGAGGCGGCAAAAAAGCCCGAGAAGGGATAgcagcccgg                                                                                                       >  2:32323/1‑145 (MQ=255)                                                         tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGT‑GGTGCCGCGCTCCGGCAAAAAAGCCCGCGGAGTGAGCGCCGCCGTGAAGGGCaaa                                                                                              >  2:113402/1‑150 (MQ=255)                                                              gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                            <  2:198587/147‑1 (MQ=255)                                                              gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTgg                                                                                          <  1:758906/149‑1 (MQ=255)                                                               cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                       <  1:27032/151‑1 (MQ=255)                                                               cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GGTCCTGCGATCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGGGAAGGACaaaacgga                                                                                         >  2:591890/1‑144 (MQ=255)                                                                 cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAGCGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACg                                                                                     <  1:566093/151‑1 (MQ=255)                                                                   aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                                   <  2:78787/151‑1 (MQ=255)                                                                        gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaaa                                                                              <  2:176267/151‑1 (MQ=255)                                                                         cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAag                                                                             <  1:567136/151‑1 (MQ=255)                                                                          gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGT‑GGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGGCAAAATGGAAACGCTGAAAag                                                                             >  1:101721/1‑150 (MQ=255)                                                                          gAAGGCCCAGATTGGCCCGGGGATGCGCTGGGGAATTGCCGAACTAAACGTCGCAGGCGAAGTGGCGGGGGGGGTGGT‑GATCCTGCGATCCGGTAAAAACGCCTGAGACGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                            <  2:136217/151‑1 (MQ=255)                                                                           aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagta                                                                           <  2:407399/151‑1 (MQ=255)                                                                             ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtagt                                                                         <  1:465739/151‑1 (MQ=255)                                                                              gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                          >  1:600507/1‑149 (MQ=255)                                                                               tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGG‑GATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGGCAAACTGGGAACGCTGAAAaggaggc                                                                        >  2:919504/1‑148 (MQ=255)                                                                                ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTg                                                                      <  2:457779/151‑1 (MQ=255)                                                                                 cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGT‑GATCCCGCGATCCGGCAAAAAAGCCCCAGAAGTGCTCGCCCCCGTGAAGGGCAAACTGGAAAAGCTGCAAAAGAGTCTGc                                                                     >  2:331541/1‑151 (MQ=255)                                                                                       ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                               <  1:579661/151‑1 (MQ=255)                                                                                       ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGCGGT‑GACTCTGCGCGTCGGCAAAAACGCCCGAGGAATGGTCGCCGCCGGGAAGGACAAAATTGGAACGCTGAAAATTTGACTGcc                                                                    >  1:950286/1‑146 (MQ=255)                                                                                            ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCg                                                            <  1:734504/149‑1 (MQ=255)                                                                                            ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                          <  1:7339/151‑1 (MQ=255)                                                                                            ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTg                                                          >  1:303957/1‑151 (MQ=255)                                                                                                          gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGATAGTTacaa                                             >  1:97429/1‑150 (MQ=255)                                                                                                           cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaaca                                           <  2:436879/151‑1 (MQ=255)                                                                                                             ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGGAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaaca                                           >  1:957734/1‑149 (MQ=255)                                                                                                               gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                       <  2:681044/151‑1 (MQ=255)                                                                                                                ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGa                                      <  1:182903/151‑1 (MQ=255)                                                                                                                ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGa                                      <  1:430483/151‑1 (MQ=255)                                                                                                                ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGa                                      <  1:150240/151‑1 (MQ=255)                                                                                                                   aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCg                                   <  1:990293/151‑1 (MQ=255)                                                                                                                      tAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGATCC‑GCGATCCGGCAAAAACGCCCGAGAAGGGAGTGCCGCCGGGGAGGAAAAACTGGAAAAGCTGAAAAGTAGGCTGCCGGGAGGGGGGGAGATAGTTACAACATACGa                                      >  2:38496/1‑145 (MQ=255)                                                                                                                       aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcag                                <  2:387182/150‑1 (MQ=255)                                                                                                                       aaaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGT‑GGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGGTCGCCGCCGTGGAGGACAAACTGGGAACGCTGAAAAGTAGGCTGCCGGGAGGGGTGGGGATAGTTACAACAAAAGAACGca                                 >  2:797140/1‑149 (MQ=255)                                                                                                                                    gaagTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc                  <  2:419251/151‑1 (MQ=255)                                                                                                                                            gggcgggGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc          <  1:609735/151‑1 (MQ=255)                                                                                                                                            gggcgggGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc          >  2:724771/1‑151 (MQ=255)                                                                                                                                                ggggTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc       >  2:63564/1‑150 (MQ=255)                                                                                                                                                ggggTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa      <  1:951688/151‑1 (MQ=255)                                                                                                                                                 gggTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaa     <  1:171296/151‑1 (MQ=255)                                                                                                                                                    tggt‑gATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCt  >  1:276054/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                        CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCT  >  NC_000913/599429‑599726

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: ‑