Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG39.6% 86.8 / 10.1 56V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (3/29);  new base G (21/0);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 6.37e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAC  >  NC_000913/599426‑599721
                                                                                                                                                      |                                                                                                                                                 
aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt                                                                                                                                                   <  2:847074/151‑1 (MQ=255)
 aaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg                                                                                                                                                  >  1:446996/1‑151 (MQ=255)
  aCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtgg                                                                                                                                                 <  2:421054/151‑1 (MQ=255)
   cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtgg                                                                                                                                                 <  2:783809/150‑1 (MQ=255)
      tCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGAGGGGGGGGGAt                                                                                                                                             >  2:116533/1‑151 (MQ=255)
         acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATTTTGCGAAGGTCCAGATTGGCCCGGAGATGCGTCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                          <  1:252457/151‑1 (MQ=255)
                tAATCCCATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCAGTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAt                                                                                                                                     <  2:566861/149‑1 (MQ=255)
                tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATc                                                                                                                                    <  1:353719/150‑1 (MQ=255)
                  aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCgg                                                                                                                                 <  2:671669/151‑1 (MQ=255)
                  aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCgg                                                                                                                                 <  1:697963/151‑1 (MQ=255)
                         cGTTTTAAGAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAAc                                                                                                                          <  2:708226/151‑1 (MQ=255)
                          gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACg                                                                                                                         <  2:893690/151‑1 (MQ=255)
                            tttAAAAGCCAGCGAAAATGGCGAGCCCGTTTATCTGCGCGATGTTGCGAAGGCCCAGCTTGGCCCGGCGATGCGCCGGGGCATTGACGATATAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTAGTCCGATTCGGGaaaaaa                                                                                                                          >  2:136619/1‑147 (MQ=255)
                                   gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAg                                                                                                                <  1:920784/151‑1 (MQ=255)
                                     cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTg                                                                                                              <  1:951232/151‑1 (MQ=255)
                                         gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcg                                                                                                          <  2:6973/151‑1 (MQ=255)
                                          aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGGGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAg                                                                                                                <  2:659187/144‑1 (MQ=255)
                                           aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAACCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgcc                                                                                                        <  2:949021/151‑1 (MQ=255)
                                             aTGGCGTGCCCGTTTATCTGCGCGATGGTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGAACTGCGGGCCGGGAAAAAAGCCCGAGAAg                                                                                                                >  2:558354/1‑141 (MQ=255)
                                               ggCGTGCCCGTTTATCTTCGCGAGGTTGGGAAGGTCCAGATTGGCACGGAGAGGCGCCGGGGGATTGGCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGGTCCTGCGATCCTGCAAAAACGCCCGAGAAGTGATCGccgccg                                                                                                    <  2:1019274/151‑1 (MQ=255)
                                                 cGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGCAAAAAAGCCCGAGAAGGGAAAgccgcc                                                                                                     >  1:1021063/1‑148 (MQ=255)
                                                   tGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGGCCCGGCAAAAAAGCCCCAGAAGTGATCGCCGCCGTGa                                                                                                 >  1:103811/1‑150 (MQ=255)
                                                   tGCCCCTCTGTCGGCGCGGTGTTGCGAAGGTCCAGATTGGCGCCGCGATGCGCCGGGGTATTGCTGAACTAAGCGGCGAAGCCGAAGGGGCGGCCGGGGTGGTGGTCCTGCGATCCGGTAAATGCGCTCCAGAAGTGCTCGCCGCCGt                                                                                                   <  2:24397/148‑1 (MQ=255)
                                                      ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTCCGGTCCGGCCAAAACGCCCGCGAAGTGATCGCCGCCGTGAAgg                                                                                              >  1:780494/1‑150 (MQ=255)
                                                        gTTGATCTGCGCGATGTTGCGAAGGTGCAAATAGCGCCGGAGATTCGCCGGGGCATTGCCGGACTAAACGGTGAATGCGAAGTGGCGGGCGGGGCGGTGACCCTGCGATCCGGAAAAAAAGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                <  2:978465/146‑1 (MQ=255)
                                                            aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAGCCGGCAAAAAAGCCCGAGGAGGGATCGGCGCCGTGAAg                                                                                               >  2:695646/1‑143 (MQ=255)
                                                            aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACaaa                                                                                         >  1:880374/1‑149 (MQ=255)
                                                                gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCGGAACTAAACGGCGAAGGCGAGGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGATAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                       <  2:638912/147‑1 (MQ=255)
                                                                      tGTTGCAAAGGCCCAGATTGGCGCGGAGATGCGGCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGTGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCt                                                                              <  2:71560/150‑1 (MQ=255)
                                                                       gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                                                                             <  1:150705/150‑1 (MQ=255)
                                                                       gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGa                                                                            <  1:197518/151‑1 (MQ=255)
                                                                         tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTACCGCGATTCCGCAAAAACGCCCGAGAAGGGAACGCCGCCGGGGAGGACAAAATAGAAACGCTg                                                                             >  2:756217/1‑148 (MQ=255)
                                                                          gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCCAAAAAGCCCGCGGAGTGATCGCCGCCGTGAAGGACAAACTGGGAAAGCTGaaa                                                                          >  1:282409/1‑150 (MQ=255)
                                                                           cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGACACCCCCATCCCGCACAAACGCCCGGGCAGGGATCGCCGCCGTGAAAGGCAAACTGGGAACGCCGa                                                                            >  2:513500/1‑147 (MQ=255)
                                                                            gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                       <  2:745093/151‑1 (MQ=255)
                                                                            gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                       <  1:639683/151‑1 (MQ=255)
                                                                            gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGGAACGCTGaaa                                                                          >  1:642619/1‑148 (MQ=255)
                                                                               ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAAATGGAAACGCTGAAAagtag                                                                     >  1:943638/1‑150 (MQ=255)
                                                                                gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc                                                                   <  1:505833/151‑1 (MQ=255)
                                                                                  ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAACCCTCGGAACGGGCAAAACGCCCCGGGGGGGAACGCCGCCGCGCAAGGGAAACTGGAAACCCTGAAAaataagt                                                                   >  2:325130/1‑146 (MQ=255)
                                                                                      aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                             <  1:1027460/151‑1 (MQ=255)
                                                                                      aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGTGGTCGCCGCCGTGGAGGACAAAATGGGAAAGCTGGAAAGTAGTCTGCCgg                                                             >  1:981246/1‑151 (MQ=255)
                                                                                         ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                             >  1:985731/1‑148 (MQ=255)
                                                                                           cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGATCCCGCCAAAACGCCCGCGAAGGGATACCCGCCCGGAAGGCCAAACTGGGAAAGCCGAAAAGTAGTCCGCCGGGAg                                                          >  2:1028495/1‑149 (MQ=255)
                                                                                             cGGAGATGCGGCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCAGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGt                                                      <  2:233127/151‑1 (MQ=255)
                                                                                                      gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCGGCAAAAACGCCCGGGAAGGGATCGCCGCCGGGAAGGACAAACTGGGAACGCTGAAAAGTAGGCTGCCGGAAGGGGTGGAGATa                                               >  1:405216/1‑149 (MQ=255)
                                                                                                        cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTa                                           <  1:401701/151‑1 (MQ=255)
                                                                                                         ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGGGAAAGtt                                            >  2:521041/1‑149 (MQ=255)
                                                                                                                 gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                  <  1:521530/151‑1 (MQ=255)
                                                                                                                  ccGAACTAAACGGCGAAGGCGAAGTGGCGGGAGGGGAGGGGCTCCTGCGAACCCGCAAAAAAGCCCGAGAAGAGATCGGCGCCGGGAAGGACAAAATGGAAAAGCTGGAAAGTAGGCTGCCGGGAGGCGGGGAGAGAGTTGCAACATAc                                   >  2:328669/1‑149 (MQ=255)
                                                                                                                   cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGCGGGGGTCCCGCGATCCGGCAAAAAAGCCCGAGAAGGGGTCGCCGCCCTGACGGACAAACCGGAAACGCTGAAAAAAAGGCTGCCGGGAGGGGCGGGGGTAGTTAAAACATAc                                   >  2:1030487/1‑148 (MQ=255)
                                                                                                                         aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                            <  1:964897/149‑1 (MQ=255)
                                                                                                                           aCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcca                        <  1:434888/151‑1 (MQ=255)
                                                                                                                               cgaaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                    >  1:549875/1‑151 (MQ=255)
                                                                                                                                gaaggcgaagTGGCGGGCGGGGGCGGGTTACTTCGATCCGGCAAAAACGCCCGGGGAGGGGTCGCCGCCGCGGAGGACAAAATGGGAACGCTGCAAAGGAGTCTGTCGGGAGGCGGGGGGGTAGTTACAAAAAAAGGTCGcagggagccc                    >  2:362458/1‑143 (MQ=255)
                                                                                                                                  aggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                 <  1:176787/151‑1 (MQ=255)
                                                                                                                                  aggcgaagGGGCGGGCGGGGTTGTGAACCTGCGAACCGGCAAAAACGCCCGAGAAGTGATCGCCGACGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGATACAACATACGATCGCAGCAAGATc                    <  2:868569/148‑1 (MQ=255)
                                                                                                                                              gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc     >  1:110524/1‑151 (MQ=255)
                                                                                                                                               ggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc     >  2:528769/1‑150 (MQ=255)
                                                                                                                                                  ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc  >  1:233236/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAC  >  NC_000913/599426‑599721

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