Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG50.0% 89.9 / 21.5 70V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (3/32);  new base G (35/0);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.50e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599427‑599725
                                                                                                                                                     |                                                                                                                                                     
aaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg                                                                                                                                                      <  1:570460/151‑1 (MQ=255)
      cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggggt                                                                                                                                                 >  2:716209/1‑146 (MQ=255)
        acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                              <  1:637514/151‑1 (MQ=255)
        acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                              <  1:1068055/151‑1 (MQ=255)
        acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                              <  2:201560/151‑1 (MQ=255)
        acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTcc                                                                                                                                               >  1:908084/1‑150 (MQ=255)
          gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTcc                                                                                                                                               >  2:359799/1‑148 (MQ=255)
            cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggggtac                                                                                                                                               >  2:411715/1‑140 (MQ=255)
               tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGaggcttgg                                                                                                                                           >  1:249150/1‑140 (MQ=255)
               tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggggtacaccgag                                                                                                                                         >  2:232548/1‑137 (MQ=255)
                aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCg                                                                                                                                      <  2:951951/151‑1 (MQ=255)
                aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCg                                                                                                                                      <  1:1033301/151‑1 (MQ=255)
                   cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGc                                                                                                                                    <  2:274135/150‑1 (MQ=255)
                    acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCGTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGGCCgg                                                                                                                                     >  2:1087200/1‑148 (MQ=255)
                       tCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGGaaaa                                                                                                                                >  1:599827/1‑150 (MQ=255)
                        cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCGGCaaaaa                                                                                                                               >  1:337302/1‑150 (MQ=255)
                         gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACg                                                                                                                             <  2:1053464/151‑1 (MQ=255)
                            ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATACTGCGATCCGGCCAAAAAGccc                                                                                                                          >  1:574735/1‑151 (MQ=255)
                              aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCg                                                                                                                         <  1:226645/150‑1 (MQ=255)
                                   ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAACCCCGAGaa                                                                                                                     >  2:868433/1‑149 (MQ=255)
                                   ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTACCGCGATCCGGCaaaaaaacccccgga                                                                                                                     >  2:983580/1‑139 (MQ=255)
                                    cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTg                                                                                                                  <  2:217429/151‑1 (MQ=255)
                                        gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                               <  1:595481/150‑1 (MQ=255)
                                          aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGAACCGGGCAAAAAGccccgggagggaacgcc                                                                                                            >  2:555965/1‑137 (MQ=255)
                                                   gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                    <  1:441681/150‑1 (MQ=255)
                                                      cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGATTCGGCAAAAACGCCCCAGAAGGGATCGCCCCCGGGaaag                                                                                                  >  1:707863/1‑147 (MQ=255)
                                                      cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAgg                                                                                                  >  2:868603/1‑149 (MQ=255)
                                                         ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCCCGGTCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGGGGAGGACaa                                                                                              >  1:618085/1‑150 (MQ=255)
                                                         ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTTCGGCAAAAACGCCCGCGAAGGGATCGCCGCCGGGAAGGACaa                                                                                              >  1:755924/1‑150 (MQ=255)
                                                               gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTg                                                                                          <  2:764328/148‑1 (MQ=255)
                                                               gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTg                                                                                          <  2:355794/148‑1 (MQ=255)
                                                               gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGATCCGGCCAAAACGCCCGAGAAGCGATCGCCGCCGTGAAGGACAAACTgg                                                                                         >  1:517681/1‑149 (MQ=255)
                                                                cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGAATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTgg                                                                                         <  1:236548/148‑1 (MQ=255)
                                                                 gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCCAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                      >  1:9366/1‑150 (MQ=255)
                                                                        tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGGGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaa                                                                               <  2:348790/150‑1 (MQ=255)
                                                                        tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATACCGCAAAAACGCCCGGAGAGGGGTCGGCGCCGGGGAGGGCCAAATGGGAAAGCTTaa                                                                               >  2:425959/1‑150 (MQ=255)
                                                                         gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaa                                                                               <  2:878251/149‑1 (MQ=255)
                                                                           gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                           <  1:436585/151‑1 (MQ=255)
                                                                             aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                         <  1:1070530/151‑1 (MQ=255)
                                                                              ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGAGCCTGCGATCCCGGAAAAAAGCCCGGGAAGGGATCGCCCCCGGGAAGGAAAAAATGGAAAAGCTGgaaa                                                                             >  2:659949/1‑146 (MQ=255)
                                                                               gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGGAGGACAAACTGGAAAAGCTGAAAagtagt                                                                        >  1:1003238/1‑150 (MQ=255)
                                                                               gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAAATGGAAAAGGTGAAAagta                                                                          >  2:32681/1‑148 (MQ=255)
                                                                                 ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                      <  2:387761/150‑1 (MQ=255)
                                                                                       tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa                                                                <  1:727740/150‑1 (MQ=255)
                                                                                       tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCgg                                                                 >  1:549093/1‑149 (MQ=255)
                                                                                        ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                                 >  1:519789/1‑148 (MQ=255)
                                                                                          cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                            <  2:607005/151‑1 (MQ=255)
                                                                                           ccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTCCGATCCCGCAAAAAAGACCGGGGAGAGATCGCCGCCGTGAAGGACAAACAGGAAACGCTGCAAAAGAGGTTGCCGGAAGGCg                                                           >  2:880725/1‑151 (MQ=255)
                                                                                           ccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCCGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGGAAAGTAGTCTGCCGGCAGGCg                                                           >  1:742753/1‑151 (MQ=255)
                                                                                            cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCCAGGAGGGATCGCCGCCGGGGAGGACAAACTGGGAAAGCTGAAAAGTAGTCTGCCGGAAgggg                                                           >  1:949699/1‑148 (MQ=255)
                                                                                              gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAAATGGAAAAGCTGAAAAGGAGTCTGCCGGGAggggggg                                                        >  1:192149/1‑146 (MQ=255)
                                                                                              gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCCGCAAAAAAGCCCGAGGAGAGATCGCCGCCGTGGAGGACAAACTGGAAAAGCTGAAAAGAAGTCTGCCGGGAGGCg                                                           >  2:623072/1‑148 (MQ=255)
                                                                                                gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                      <  1:359862/151‑1 (MQ=255)
                                                                                                 aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                       <  2:85849/149‑1 (MQ=255)
                                                                                                   gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                    <  2:753301/150‑1 (MQ=255)
                                                                                                   gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                    <  1:672827/150‑1 (MQ=255)
                                                                                                   gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGGAGGCGGGgag                                                      >  1:674007/1‑148 (MQ=255)
                                                                                                       cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGAACCTCCGATCCGGCGAAAACGCCCGAGGAGGGAACGCCGCCGCGGCGGAACAAATGGGAACCGTGTTCACGAGTCTGCCGGAAGGGGTGGGGGTAGt                                                 >  2:98246/1‑149 (MQ=255)
                                                                                                        ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGAGCCTGCGATCCGGCCAAAACGCCCGAGAAGTGATCGCCGCCGGGGAAGACAAACTGGAAAAGCTGGAAAGTAGGCTGGCCGAAGGGGGGGGGATAGtt                                                >  1:1019831/1‑149 (MQ=255)
                                                                                                             aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                         <  2:18818/151‑1 (MQ=255)
                                                                                                                 ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGGAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATACg                                      >  1:630430/1‑150 (MQ=255)
                                                                                                                  cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGACCTGCGATCCGGCCAAAACCCCCGCGAAGAGATCGCCGCCGCTGAGGACAAACTGGGAACGCCGGAAAGAAATCTGCCGGAAGGGGTGGGGATAGGTACCAAATACGa                                     >  2:603826/1‑150 (MQ=255)
                                                                                                                   gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                   <  2:218063/151‑1 (MQ=255)
                                                                                                                    aaCTAACCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCg                                  <  2:244565/151‑1 (MQ=255)
                                                                                                                    aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGCGGAGGGGTCGCCGCCGTGAAGGGCAAAATGGAAAAGCTGAAAAGGAGTCTTCCGGGAGGGGGGGGGATAGTTACAACATACGaac                                   >  2:618349/1‑148 (MQ=255)
                                                                                                                       tAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGAGCCGTCGATCCGGCAAAAACCACCGGGAAGTGAACGCCGCCGTGGAGGGCAAACTGGAAAAGCTGGAAAAAAGTCTGCCGGAGGGGGTGgaggtggattta                                             >  2:407763/1‑128 (MQ=255)
                                                                                                                            ggcgaaggcgaaGTGGCGGGGGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGAGATAGTTACAACATACGATCGcagcca                            >  1:782953/1‑149 (MQ=255)
                                                                                                                            ggcgaaggcgaaGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGGAGGAAAAAATGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGATAGTTTCAACATACGAACGcagcc                             >  2:536626/1‑148 (MQ=255)
                                                                                                                               gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTc                        <  2:629471/150‑1 (MQ=255)
                                                                                                                                    cgaagTGGCGGGCGGGGTGGGGATTCTGCGATCCGGCAAAAAAACCCGAGAAGGGATCGCCGCCGTGGAGGAAAAACCGGAAACGCCGAAAAGGAGTCTGCCGGAAGGCGGGGGGATAGGTACAAAATACGAAAGCCGCCAGCTCATTGAc                  >  2:652899/1‑151 (MQ=255)
                                                                                                                                        gTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgc               <  2:404795/150‑1 (MQ=255)
                                                                                                                                              ggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCg        <  2:544450/151‑1 (MQ=255)
                                                                                                                                                 ggggTGGTGAACCTGCGAGCCGGCAAAAACGCCCGAGAAGTGAGCGCCGCCGTGAAGGACAGGCTGGAAACGCTGAAAAGTGGTCTGCCGGAAGGCGTGGGGATAGTTACACCATTCGATCGCAGCCAGCTCATTGACCGCGCCa           <  2:694543/145‑1 (MQ=255)
                                                                                                                                                  gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa     >  2:748498/1‑150 (MQ=255)
                                                                                                                                                    gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc  <  2:337239/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599427‑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: