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

AACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAC  >  NC_000913/599427‑599724                                                                                                                                                      |                                                                                                                                                       aaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg                                                                                                                                                       <  1:148918/151‑1 (MQ=255)   cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggg                                                                                                                                                      >  1:121409/1‑150 (MQ=255)      tCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg                                                                                                                                                    >  1:240317/1‑147 (MQ=255)        gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGggggtgggcc                                                                                                                                                >  2:771462/1‑147 (MQ=255)                 aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑g                                                                                                                                     <  1:221338/151‑1 (MQ=255)                    cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGc                                                                                                                                   <  1:125555/150‑1 (MQ=255)                      caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCaaa                                                                                                                                <  2:140273/151‑1 (MQ=255)                      caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CC‑GGCaa                                                                                                                                 >  1:309027/1‑150 (MQ=255)                       aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCaaaa                                                                                                                               <  1:234283/151‑1 (MQ=255)                         cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAAc                                                                                                                             <  1:74700/151‑1 (MQ=255)                          gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACg                                                                                                                            <  1:142975/151‑1 (MQ=255)                             ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGcc                                                                                                                          <  1:194373/150‑1 (MQ=255)                               aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGGT‑CC‑GGCCAAAACGccc                                                                                                                         >  1:564562/1‑149 (MQ=255)                               aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGAT‑CC‑GGCCAAAAAGccc                                                                                                                         >  1:333587/1‑149 (MQ=255)                               aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGAT‑CC‑CGCAAAAAAGccc                                                                                                                         >  2:564968/1‑149 (MQ=255)                                   gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGACCCGCGAT‑AC‑GGCacaaa                                                                                                                              >  2:585885/1‑140 (MQ=255)                                            aaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTACCGCGGT‑CCGGGCAAAAAGGccgcgaa                                                                                                                     >  1:313986/1‑136 (MQ=255)                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                        <  1:687732/150‑1 (MQ=255)                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CC‑GGCCAAAACGCCCGAGAAGGGATCgccgcc                                                                                                        >  1:279005/1‑150 (MQ=255)                                                    gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAg                                                                                                  <  1:302392/151‑1 (MQ=255)                                                      ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGT‑CC‑GGCAAAAACGCCCGCGAAGTGATCGCCGCCGTGAAgg                                                                                                 >  2:300423/1‑150 (MQ=255)                                                      ccGTTTATCTGCGCGATGTTGCGAAGGGCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGGGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGa                                                                                                <  1:126159/151‑1 (MQ=255)                                                          ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaa                                                                                             <  1:128639/150‑1 (MQ=255)                                                                    gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGCC‑CC‑GGGAAAAAAGCCCGCGAAGTGGTCGCCGCCGGGAAGGGCAAAGCGGGAACGc                                                                                  >  2:456969/1‑151 (MQ=255)                                                                      tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGT‑CC‑GGCAAAAAAGCCCGCGGAGGGATCGCCGCCGGGGAGGACAAAATGGgaa                                                                                     >  2:106422/1‑146 (MQ=255)                                                                       gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAACCGCTGa                                                                               <  2:275462/151‑1 (MQ=255)                                                                       gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                                                                                <  1:513153/150‑1 (MQ=255)                                                                       gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CC‑GGCAAAAACGCCCGGGAAGGGATCGCCGCCGGGAAGGACAAAATGGAAACGCTg                                                                                >  1:674819/1‑150 (MQ=255)                                                                         tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CC‑GGCAAAAAAGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGaa                                                                              >  1:583091/1‑150 (MQ=255)                                                                          gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGG‑CC‑CGCAAAAAAGCCCGAGAAGTGATCGCCGCCGCGAAGGACAAACTGGGAACGCt                                                                                 >  1:326322/1‑146 (MQ=255)                                                                            gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAag                                                                           <  2:511539/150‑1 (MQ=255)                                                                              aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGGAGAGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaaa                                                                            >  1:673504/1‑147 (MQ=255)                                                                               ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGCT‑CC‑GGCAAAAACGGCCGCGAAGGGATCGGCGCCGCGGAGGACAAACTGGGAAAGCTGcaaaagag                                                                        >  2:426009/1‑146 (MQ=255)                                                                               ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGaaa                                                                             >  2:454361/1‑145 (MQ=255)                                                                                  ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTg                                                                    <  1:346506/151‑1 (MQ=255)                                                                                   cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGc                                                                   <  1:745910/151‑1 (MQ=255)                                                                                    aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGcc                                                                  <  1:728921/151‑1 (MQ=255)                                                                                    aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGCCCCTCGGATCCC‑GGCAAAAA‑GACCCAGGAGTGGATGGCGCCGGTACGGGACAACTGGAAACGCTGGAAAGTAGTCGGcc                                                                  >  2:88435/1‑151 (MQ=255)                                                                                    aGATTGGCCAGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGCCCGCCAT‑CA‑GGCAAAAACGCACGGGAAGGGAGTGCCGCCGGGAAGGGCAAACCGGGAAACCGGAAAagtcggct                                                                     >  2:88770/1‑143 (MQ=255)                                                                                      aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCCGCGAT‑AC‑GGCAAAAAAGCCCGGGAAGTGATTGCGGCCGCGGAGGACAAAATGGGAACACTGGAAAATAGTCTGcc                                                                  >  2:202354/1‑149 (MQ=255)                                                                                        tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa                                                               <  2:140275/150‑1 (MQ=255)                                                                                          gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                            <  1:390268/151‑1 (MQ=255)                                                                                             cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGTCCTGCGAC‑CC‑GGCAAAAACGCCCGAGGAGTGATAGCCGCCGGGAAGGGCAAACTGGAAAAGCCGAAAAGGAGTCTGCCGGAAGGCg                                                          >  2:758647/1‑150 (MQ=255)                                                                                              ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                        <  2:222895/151‑1 (MQ=255)                                                                                                           ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                           <  1:115165/151‑1 (MQ=255)                                                                                                                 gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGCCTTCGAT‑CC‑CGGAAAAAAGCCCGAGAAGGGGTCGCCGCCGTGAAGGACAAACTGGAAACGCCGGAAAGGAGGCTGCCCGGAGGCGTGGGGGTAGTTACAACCTACg                                     >  2:267662/1‑151 (MQ=255)                                                                                                                 gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAAAACGCGAT‑CC‑GGCAAAAACGCCCGCGAAGTGGTCGCCGCCGGGAAGGAGCAACCGGGAAAGCTGCAAAGAAGGCTGCCCGAAGGGGGGGAGATAGTTacaacaca                                       >  2:535554/1‑147 (MQ=255)                                                                                                                   cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                   <  2:464307/151‑1 (MQ=255)                                                                                                                    gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGATAGTTACAACATACGa                                    >  1:402078/1‑149 (MQ=255)                                                                                                                         aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                               <  1:487074/149‑1 (MQ=255)                                                                                                                             ggcgaaggcgaaGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAAAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagccagc                         <  1:65605/151‑1 (MQ=255)                                                                                                                                   ggcgaagTGGGGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAttt                   <  2:165338/151‑2 (MQ=255)                                                                                                                                     cgaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAc                 <  2:453362/151‑1 (MQ=255)                                                                                                                                      gaagTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc                <  2:158269/151‑1 (MQ=255)                                                                                                                                            gcgggcggGGCGGCGATCCTGCCAT‑CC‑GGCAAAAACGCCCGAGAAGCGATCGCCGCCGTGAAGGACAAACTGGAAACGCCGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGcc           >  2:509351/1‑150 (MQ=255)                                                                                                                                              gggcgggGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc        <  2:375524/151‑1 (MQ=255)                                                                                                                                              gggcgggGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc        >  2:774950/1‑151 (MQ=255)                                                                                                                                                  ggggTGCAGATCATGCGAT‑AC‑GGCAAAAGCGCCCGAGACGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGCAGTCTGCCGGTAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGCCCGCGCCATCGAc     >  2:737100/1‑150 (MQ=255)                                                                                                                                                    ggtggtGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAc  >  1:578398/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       AACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAT‑CC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAC  >  NC_000913/599427‑599724

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