Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	599,576	0	T	G	50.7%	66.6 / 111.1	68	V288G (GTG→GGG)	cusA	copper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (1/32);  new base G (34/0);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 2.46e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-11

CGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599433‑599722                                                                                                                                                |                                                                                                                                                      cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggg                                                                                                                                                     >  2:582685/1‑146 (MQ=255) cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggggtc                                                                                                                                                >  2:609734/1‑146 (MQ=255)  gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAACCGGCGAAGGCGAAGTGGCGGGCGgggggggg                                                                                                                                                   >  2:454287/1‑145 (MQ=255)   acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                              <  2:1042320/151‑1 (MQ=255)    cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑g                                                                                                                                            <  2:609493/151‑1 (MQ=255)    cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGCAGgggggggggaa                                                                                                                                                >  2:161708/1‑143 (MQ=255)     gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑Gc                                                                                                                                           <  2:610758/151‑1 (MQ=255)          tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATc                                                                                                                                      <  2:1055386/150‑1 (MQ=255)           aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATcc                                                                                                                                     <  1:1030188/150‑1 (MQ=255)           aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCg                                                                                                                                    <  1:293353/151‑1 (MQ=255)            aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCC‑GGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCT‑GC‑GATCCgg                                                                                                                                   >  2:87430/1‑150 (MQ=255)             tCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTccctccgtcc                                                                                                                                       >  1:491477/1‑140 (MQ=255)              cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCT‑GC‑gggcc                                                                                                                                     >  1:747500/1‑143 (MQ=255)              cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGc                                                                                                                                  <  2:841669/150‑1 (MQ=255)                caTCGTTTTAAAAGCCAGTGAAAATGGCTTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCaaa                                                                                                                               <  1:589986/151‑1 (MQ=255)                caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCaaa                                                                                                                               <  2:880067/151‑1 (MQ=255)                      tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCT‑TC‑GATCCCGGAAAAACGcc                                                                                                                         >  1:555270/1‑151 (MQ=255)                      tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCT‑GC‑GATCCGGCAAAAACGcc                                                                                                                         >  2:688223/1‑151 (MQ=255)                      tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GC‑GATCCGGCAAAAACGcc                                                                                                                         >  2:783434/1‑151 (MQ=255)                         aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCg                                                                                                                       <  2:696964/150‑1 (MQ=255)                         aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCT‑GC‑GGTCCCGCAAAAACGccc                                                                                                                        >  2:126794/1‑149 (MQ=255)                          aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGACGAAGGCGAAGTGGCGGGCGGGGGGGGGGCCCT‑GCGGCCCCGGCaaaaagcc                                                                                                                          >  1:241509/1‑144 (MQ=255)                                aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCC‑TC‑GAGCCGGCAAAAAAGcc                                                                                                                         >  1:464832/1‑141 (MQ=255)                                 gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGAGTCT‑GC‑GATCCGGCACAAAAGCCCGGGAA‑GTg                                                                                                               >  1:62971/1‑149 (MQ=255)                                  tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCC‑GC‑GATCCGGCAAAAAAGCCCGAGAA‑GTg                                                                                                               >  1:1059176/1‑148 (MQ=255)                                  tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCT‑GC‑GATCCGGCAAAAAAGCCCGAGAA‑ggg                                                                                                               >  2:362072/1‑146 (MQ=255)                                  tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GC‑GAGCCGGCAAAAACGGCCGCGAA‑ggg                                                                                                               >  1:1089284/1‑146 (MQ=255)                                               cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAgg                                                                                               <  1:710511/151‑1 (MQ=255)                                                    ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACaa                                                                                           <  2:893217/150‑1 (MQ=255)                                                        cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GC‑GATCCGGCGAAAACGCCCGGGAA‑GGGATCGCCGCCGGGAAGGACAAACTgg                                                                                      >  2:359264/1‑151 (MQ=255)                                                        cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GC‑GATCCCGCCAAAACGCCCGCGAA‑GGGATCGCCGCCGGGAAGGACAAACCgg                                                                                      >  1:1068782/1‑151 (MQ=255)                                                          gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGGGAAGGCGAAGTGGCGGGCGGGTTGGTGATCCT‑GC‑GATCCGGCAAAAGCGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAGCt                                                                                        <  1:829449/147‑1 (MQ=255)                                                          gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACt                                                                                        <  2:1040405/147‑1 (MQ=255)                                                           cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCT‑GC‑GCGCCCGCAAAAACGCCCCCGAA‑GCGATCGCCGCCCGGGAGGACCAAATGGaa                                                                                    >  1:663856/1‑150 (MQ=255)                                                           cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GC‑GAACCCGCAAAAACGCCCGAGGA‑GAGAATGCCGCCCTGAAGGGCAAACTGGaa                                                                                    >  1:203053/1‑150 (MQ=255)                                                             cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGGG‑GAGCGGGCAAAAA‑GCCCGAGAAGGGGATCGCCCCGGGAAAGGAGAACCTGaaa                                                                                    >  1:190506/1‑149 (MQ=255)                                                              gATGTTGGGAAGGTGCAGATTGGCCCGGAGATGCGCGGGGCGATTGCCGAACTAAACGGCGAAGGCGAAGTGGGGTGCGGGGTGGTGATCCT‑TC‑GATCCGGCAAAATCCCCCGAGAT‑GTGATCGCCGCCGTGAAGGACAAACTGCATACg                                                                                 <  1:885213/150‑1 (MQ=255)                                                               aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATACC‑CC‑GAACCCGAAAAAACAACCGGGAA‑AGGGTCGCCGCCGTGGAGGGCAAACAGGaaa                                                                                   >  1:825999/1‑147 (MQ=255)                                                                 gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCT‑TC‑CATCCGGCAAAAACGCCCGAGAA‑GGGATCGCCGCCGTGAAGGGCAACCTGGAAAAGCTg                                                                              >  1:41513/1‑150 (MQ=255)                                                                  ttGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGGTCCT‑GC‑GATCCGGCAAAAACGCCCGGGAA‑GGGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTg                                                                              >  2:12654/1‑149 (MQ=255)                                                                  ttGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCT‑GC‑GATCCGGCAAAAAAGCCCGCGAA‑GTGAGCGCCGCCGTGAAGGACAAACTGGAAAAGCTg                                                                              >  1:144671/1‑149 (MQ=255)                                                                          gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCAT‑GC‑GAACCGGCAAAAAACCCCGAGAA‑GGGGTCGCCGCCGTGAAGGGGAAAATGGAAAAGCTGaaaa                                                                          >  1:1075897/1‑145 (MQ=255)                                                                            ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                   <  2:518971/150‑1 (MQ=255)                                                                                aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCT‑GC‑GGTCCGGCCAAAAAGCCCGAGAA‑GGGGTCGCCGCCGGGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGcccg                                                              >  2:81640/1‑149 (MQ=255)                                                                                  tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa                                                             <  1:239399/150‑1 (MQ=255)                                                                                   ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                           <  2:57705/151‑1 (MQ=255)                                                                                        ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                      <  1:994950/151‑1 (MQ=255)                                                                                          agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCCTGAAGGACAAACTGGAAAAGCTGAAAAAGAGTCTGCCGGAAGGCGTgg                                                     >  2:89163/1‑150 (MQ=255)                                                                                              gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                 <  2:302922/150‑1 (MQ=255)                                                                                              gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                <  1:52870/151‑1 (MQ=255)                                                                                              gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                <  2:788476/151‑1 (MQ=255)                                                                                               cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAg                                               <  1:69471/151‑1 (MQ=255)                                                                                                gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGACCT‑GC‑GATCCGGCAAAAAAGCCCGGGGA‑GGGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCgggggggta                                                >  1:90779/1‑144 (MQ=255)                                                                                                 ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCT‑TC‑GGTCCCGCAAAAACGCCCGAGAA‑GGGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGGGATAGtt                                             >  2:729583/1‑151 (MQ=255)                                                                                                  cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTa                                            <  1:41810/151‑1 (MQ=255)                                                                                                  cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCT‑GC‑GATCCGGGAAAAACGCCCGAGAA‑AAGGTCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAggt                                             >  2:395998/1‑148 (MQ=255)                                                                                                  cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCT‑GC‑GATCCGGCAAAAACGCCCGGGGA‑GGGATAGCCGCCGGGGAGGAAAAACTGGAAAAGCTGAAAAGGAGTCTGCCGGAAGGGGGGGAGATAggt                                             >  1:756468/1‑148 (MQ=255)                                                                                                  cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GC‑GATCCGGCAAAAACGCCCGGGAA‑GGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAAAAGTCTGCCGGAAGGCGTGGGGATAGtt                                             >  2:725526/1‑150 (MQ=255)                                                                                                    gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTaca                                          <  2:588379/151‑1 (MQ=255)                                                                                                     ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                         <  1:1039623/151‑1 (MQ=255)                                                                                                         ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GC‑GAACCGGGAAAAACGCCCGAGAA‑GTGATCGCCGCCGGGGAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGGGTGGAGATAGTTacaa                                         >  2:490152/1‑147 (MQ=255)                                                                                                         ttGACGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCT‑GC‑GAACCCGCAAAAACGCCCGAGAA‑GGGATCGCCCCCGGGAAGGACAAAATGGAAAAGCTGAAAAGTAGGCTGCCGGGAGGGGTGGAGATAGGTacaaca                                       >  1:645768/1‑149 (MQ=255)                                                                                                           gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                   <  2:707344/151‑1 (MQ=255)                                                                                                           gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGAGTCT‑CC‑GATCCCGCAAAAAAGCCCGCGGA‑AGGGTCGGCGCCGCGGTTGACCAACTGGAAACGCCGCAAAGGAGCGACCCGGCAGGCGGGGGGATGGTTacaac                                        >  1:963271/1‑146 (MQ=255)                                                                                                            ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCG‑GC‑GGTCCCGGCAAAACGCCCGGGAA‑GTGAACGCCGCCCCCAGGGAGCAACCGGAAAAACCGGAAAGAAGTCTGCCGGCAAGCCTCGGGATATTTTCAACCTAc                                    >  1:507819/1‑149 (MQ=255)                                                                                                              gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GC‑GATCCGGCAAAAAAGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATACGa                                  >  2:719678/1‑149 (MQ=255)                                                                                                                          gaaggcgaagTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                     <  1:582418/150‑1 (MQ=255)                                                                                                                                   gTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgc            <  1:892243/150‑1 (MQ=255)                                                                                                                                      gcgggcggGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCa        <  2:863764/151‑1 (MQ=255)                                                                                                                                       cgggcgggGAGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGCGGAGATAGTTACAACATACGATTGCAGCCAGCTCATTGACCGCGCCAt       >  1:156100/1‑151 (MQ=255)                                                                                                                                           cgggGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc   <  2:999769/151‑1 (MQ=255)                                                                                                                                            ggggTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa  <  2:364086/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      CGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GC‑GATCCGGCAAAAACGCCCGAGAA‑GTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599433‑599722

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