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

AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599426‑599725                                                                                                                                                       |                                                                                                                                                       aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt                                                                                                                                                        <  1:477488/151‑1 (MQ=255)  aaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGAGTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg                                                                                                                                                       <  1:238186/151‑1 (MQ=255)          acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTcc                                                                                                                                                >  2:896208/1‑150 (MQ=255)          acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                               <  1:271659/151‑1 (MQ=255)                ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GAt                                                                                                                                         <  2:256286/150‑1 (MQ=255)                 tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATc                                                                                                                                        <  1:63694/150‑1 (MQ=255)                 tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGC‑GGTcc                                                                                                                                       >  1:430887/1‑151 (MQ=255)                          cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTTC‑GGCCCGGGaaaaa                                                                                                                               >  2:1120589/1‑150 (MQ=255)                                aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGCTGTTGCGAAGGTCCAGATTGGCCCGGTGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCCGCGGATCCGCCAAAAAAGcc                                                                                                                           >  2:606015/1‑149 (MQ=255)                                aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCg                                                                                                                         <  1:484078/150‑1 (MQ=255)                                  aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGC‑GGGCCCGCcaaaaaaccc                                                                                                                          >  2:517701/1‑142 (MQ=255)                                    gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGC‑GGTCCGGCaaaaaacccc                                                                                                                          >  2:1172848/1‑140 (MQ=255)                                       aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGa                                                                                                                 <  2:636379/151‑1 (MQ=255)                                       aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCGTGC‑GGACCCGCCAAAAAGCCCgggga                                                                                                                     >  1:645190/1‑143 (MQ=255)                                        gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGa                                                                                                                 <  1:901597/150‑1 (MQ=255)                                        gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGGTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATGCCGC‑GGTCCGGCCAAAAAACCCgcgggatggt                                                                                                                >  1:207856/1‑142 (MQ=255)                                           aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GATCCGGCCAAAACGCCCGAGAAGGGATcg                                                                                                              >  1:1094353/1‑150 (MQ=255)                                              aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                                                                                          <  2:115930/151‑1 (MQ=255)                                                ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                         <  1:988990/150‑1 (MQ=255)                                                ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                         <  1:1017261/150‑1 (MQ=255)                                                   gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa                                                                                                     <  1:230686/151‑1 (MQ=255)                                                     gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GAACCGGCCAAAAACCCCGAGAAGTGATCGCCCCCGTg                                                                                                      >  1:606694/1‑148 (MQ=255)                                                      cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGC‑CGTTCGGGAAAAAAGCCCCAGAAGGGATCGcccccgcg                                                                                                      >  1:886818/1‑145 (MQ=255)                                                         gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GGTCCGGCAAAAACGGCCGGGAAGGGATCGCCGCCGGGAAgggc                                                                                                >  1:572697/1‑148 (MQ=255)                                                           ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaa                                                                                              <  1:1068091/150‑1 (MQ=255)                                                             aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCGGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAc                                                                                            <  2:831023/150‑1 (MQ=255)                                                             aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                             >  1:914598/1‑149 (MQ=255)                                                                 gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGa                                                                                        <  1:1096339/150‑1 (MQ=255)                                                                       tGTTGCTAAGGTCCGGATTGGCCCGGAGATGCGCCGGTGCATTGCCGAACTAAACGGCGAAGTCAAAGTGGCGGGCGGGGTGGTGAGCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                     <  2:1098764/147‑1 (MQ=255)                                                                        gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                                                                                 <  1:409611/150‑1 (MQ=255)                                                                        gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                                                                                 <  2:1005884/150‑1 (MQ=255)                                                                        gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGa                                                                                <  1:805502/151‑1 (MQ=255)                                                                         ttGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGATCCTGC‑GATCCGGCCAAAAAACCCGAGAAGTGATCGCCGCCGCGAAGGACAAAATGGGAAAGccg                                                                                 >  2:166141/1‑147 (MQ=255)                                                                            cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAag                                                                            <  2:1051247/151‑1 (MQ=255)                                                                            cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAag                                                                            <  2:921143/151‑1 (MQ=255)                                                                              aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GATCCGGCAAAAAAGCCCGAGAAGAGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGaaaa                                                                             >  1:1139380/1‑148 (MQ=255)                                                                               aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                         <  1:57837/151‑1 (MQ=255)                                                                                  tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGCAGGACAAAATGGAAACGCTGAAAaggagt                                                                        >  1:392526/1‑149 (MQ=255)                                                                                    cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTTC‑GGTCCCGCAAAAACGCCCGGGAAGAGATCGCCGCCGCGAAGGGAAAAACGGAAACGCGGgaaa                                                                             >  2:906835/1‑142 (MQ=255)                                                                                    cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGC‑GATCCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGaaaaaaa                                                                          >  1:750861/1‑142 (MQ=255)                                                                                     aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGC‑GATCCGGCAAAAACGCCCGAGAAGGGGTCGCCGCCGGGGAGGACAAACTGGAAAAGCTGAAAAGGAGTCTGcc                                                                   >  1:990246/1‑151 (MQ=255)                                                                                          ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                              <  2:112830/151‑1 (MQ=255)                                                                                           gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGCCTGC‑GCTTCGGCGAAAAAACCCGCGGAGTGATCGGCGCCCCGAAGGGCAAACTGGAAACGCCGCAAAGAAGTCTGCCgg                                                                 >  2:704670/1‑147 (MQ=255)                                                                                              cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGCGGGGTTGCTGC‑GATCCGGCGAAAAAGGCCGGGAAGTGTTCGCCGCCGGGGAGGACAAAATGGAAAGGCTGGAAAGTAACCTGCCGGaa                                                               >  1:58041/1‑146 (MQ=255)                                                                                                 agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                       <  1:750740/151‑1 (MQ=255)                                                                                                  gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGC‑GATCCGGCCAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAAAGCCGAAAAAGAGTCTGCCGGGAggggggg                                                        >  1:941102/1‑144 (MQ=255)                                                                                                   aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GGTCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGGAACGCTGAAAAGTAGTCTGCCGGGAGGCGTgg                                                        >  1:957927/1‑148 (MQ=255)                                                                                                       gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                   <  2:1130668/149‑1 (MQ=255)                                                                                                         cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGC‑GATCCGGCAAAAACGCCCCAGAAGGGATCGCCCCCGTGGAGGACAAACTGGGAACGCTGAAAAATAGTCTGCCGGGAGGCgggg                                                        >  2:840122/1‑142 (MQ=255)                                                                                                          ggggTATTGCCGAACTAAACGGGGAAGGGGAAGTGGCGGGCGGGGTGGTGATCCGGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAACCGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                    <  2:912294/145‑1 (MQ=255)                                                                                                            ggCATTGCCGAACTAAACGACGAAGGCGAGGTGGCGGGCGGGGTGGTGATCCTGC‑GTTCCGGCAAAAACGCCCGAAAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCTGAAGGCGTGGAGATAGTTacaa                                            <  2:1081256/151‑1 (MQ=255)                                                                                                                ttGCCTAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GGTCCGGCCAAAAAGCCCGAGAAGGGATTGCCGCCGGGAAAGACAAACTGGAAACGCTGAAAAGTAGGCTGCCGGGAGGGGTGGGGATAGGTacaaca                                          >  1:65902/1‑149 (MQ=255)                                                                                                                ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCCC‑GATCCCGCAAAAAAACCCGAGAAGGGAATGCCGCCGTGGAGGACAAACTGGAAAAGCCGGAAAGTAGTCTGCCGGGAGGGGGGGGGATAGATacaaaa                                          >  2:151847/1‑147 (MQ=255)                                                                                                                         tAAACGGCGAAGGCGAAGTGGCGGGAGGGGGGGCGATCCTGC‑GACCCGGCAAAAACGCCCGGGAAGGGATCGCCGCCGGGAAGGACAAAATGGGAACGCTGGAAAGGAGTCTGCCGGGAGGGGGGGAGATAGTTACAACATACGATCGc                                 >  1:234478/1‑149 (MQ=255)                                                                                                                          aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCg                                  <  2:68334/147‑1 (MQ=255)                                                                                                                           aaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GATCCGGCCAAAACGCCCGAGAAGGGATCGCCGCCGCGAAGGGCAAACTGGAAACGCTGAAAAGTAGTCCGCCGGGAGGCGGGGGGATAGTTACAACATACGATCGcagcc                             >  1:962397/1‑151 (MQ=255)                                                                                                                               gcgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCt                         <  1:786054/151‑1 (MQ=255)                                                                                                                                 gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                        <  2:142570/150‑1 (MQ=255)                                                                                                                                       gaagTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAc                  <  2:666231/150‑1 (MQ=255)                                                                                                                                           tGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCgcgc             >  1:93029/1‑151 (MQ=255)                                                                                                                                               gggcgggGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc         <  2:946144/151‑1 (MQ=255)                                                                                                                                                ggcgggGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAtag        <  2:683543/151‑3 (MQ=255)                                                                                                                                                      gtggtgATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc  >  2:256983/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599426‑599725

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