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

CGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCT  >  NC_000913/599433‑599726
                                                                                                                                               |                                                                                                                                                       
tgacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggggg                                                                                                                                                     >  2:489290/2‑146 (MQ=255)
cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                                                                                                                                                    <  1:64518/149‑1 (MQ=255)
  acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTcc                                                                                                                                                 >  1:893966/1‑150 (MQ=255)
  acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAAATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggg                                                                                                                                                       >  2:821426/1‑144 (MQ=255)
   cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggggtgccg                                                                                                                                               >  2:530508/1‑143 (MQ=255)
    gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATcc                                                                                                                                                 >  1:50059/1‑148 (MQ=255)
        ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCc‑gcggt                                                                                                                                          >  2:425669/1‑148 (MQ=255)
        ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAtgccg                                                                                                                                               >  2:689432/1‑142 (MQ=255)
         tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCT‑TCggatc                                                                                                                                        >  2:650410/1‑147 (MQ=255)
          aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATcc                                                                                                                                        <  1:451389/150‑1 (MQ=255)
          aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GCGATcc                                                                                                                                        >  1:1031141/1‑150 (MQ=255)
           aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGTTGCC‑GCGATGCgg                                                                                                                                      >  2:232129/1‑151 (MQ=255)
              acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCT‑GCGGGCCgg                                                                                                                                      >  2:143715/1‑148 (MQ=255)
                  cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCaaaaa                                                                                                                                <  2:186407/150‑1 (MQ=255)
                  cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCaaaaa                                                                                                                                <  1:339363/150‑1 (MQ=255)
                  cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCaaaaa                                                                                                                                <  2:959501/150‑1 (MQ=255)
                  cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCT‑GCGATCCGGCaaaaa                                                                                                                                >  1:268928/1‑150 (MQ=255)
                   gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCT‑GCGATCCGGCaaaaa                                                                                                                                >  1:404374/1‑149 (MQ=255)
                        aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACCAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCg                                                                                                                          <  2:1031623/150‑1 (MQ=255)
                         aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCC‑GCGGTCCCGCAAAAAAGccc                                                                                                                           >  1:962551/1‑148 (MQ=255)
                           aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGaa                                                                                                                      <  1:1066936/151‑1 (MQ=255)
                             ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGt                                                                                                                    <  1:833653/151‑1 (MQ=255)
                             ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCT‑GCGGTCCCGGAAAAACGccc                                                                                                                           >  2:377946/1‑144 (MQ=255)
                               aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGa                                                                                                                  <  1:1005463/151‑1 (MQ=255)
                               aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGCCT‑GCGGTCCGGCGAAAAACCCCCAAAAGTg                                                                                                                   >  2:883373/1‑150 (MQ=255)
                                  gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCT‑GCGGTCCGGCAAAAAAACCCCGGAAGGGAtt                                                                                                                >  2:647281/1‑149 (MQ=255)
                                    aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGCAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTCGGAATCCGCCAAAAAccccgggca                                                                                                                       >  2:137112/1‑137 (MQ=255)
                                    aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATcgcc                                                                                                             <  2:920316/151‑1 (MQ=255)
                                    aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCC‑GCGATCCGGCAAAAACGCCCGGGAAggg                                                                                                                   >  2:639247/1‑143 (MQ=255)
                                      aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCT‑TCGATCCGGCAAAAACGCCCGGGAAGAGATcgcc                                                                                                             >  1:972365/1‑149 (MQ=255)
                                        ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                          <  1:200183/150‑1 (MQ=255)
                                        ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGccgccg                                                                                                         <  1:542841/151‑1 (MQ=255)
                                         gCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑TCGATCCCGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                          >  1:741321/1‑149 (MQ=255)
                                           gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa                                                                                                      <  1:442666/151‑1 (MQ=255)
                                             gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGACCT‑GCGAACCGGCAAAAAAACCCGAGCAGTGATCGccgccggg                                                                                                       >  2:867895/1‑146 (MQ=255)
                                              cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                                   <  1:266433/151‑1 (MQ=255)
                                               ccGTTTATCTGCGCGATGTTGCGAAGGTGCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGa                                                                                                  <  2:942216/151‑1 (MQ=255)
                                                    tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                              <  1:31633/150‑1 (MQ=255)
                                                       cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTACT‑GTGAGCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGGGAAGGGGAAACTgg                                                                                          >  1:421583/1‑151 (MQ=255)
                                                         gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                            <  1:441704/147‑1 (MQ=255)
                                                         gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTgg                                                                                          <  1:676181/149‑1 (MQ=255)
                                                           gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                      <  1:506851/151‑1 (MQ=255)
                                                           gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCT‑GCGATCCGGCAAAAAAGCCCCAGAAGGGATCGCCGCCGGGAAGGACAAACTGGaaa                                                                                       >  1:685845/1‑150 (MQ=255)
                                                                      aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAag                                                                             >  1:530896/1‑149 (MQ=255)
                                                                          tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GCGATCCGGGAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGAAAAACTGGAAAAGCTGAAAagta                                                                           >  1:882480/1‑147 (MQ=255)
                                                                             aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCT‑TCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAAGACAAACTGGGAAAGCCGAAAAGTAGTCTGcc                                                                    >  1:408972/1‑151 (MQ=255)
                                                                              gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCg                                                                   <  2:718048/151‑1 (MQ=255)
                                                                               aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GCGATCCGGCAAAAAAGCCCGCGAAGGGATCGCCGCCGTGAAGGAAAAACTGGAAACGCTGAAAAGGAGTCTGCCgg                                                                  >  1:56068/1‑151 (MQ=255)
                                                                                  ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                               <  2:336304/151‑1 (MQ=255)
                                                                                     ccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GCGATCCGGCAAAAACGCCCGGGAAGGGATTGCCGCCGGGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCg                                                            >  1:128313/1‑151 (MQ=255)
                                                                                      cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GCGAGCCGGGAAAAAAGCCCGAGGAGTGATCGCCGCCGTGGAGGAACAACTGGAAAAGCCGAAAAGTAGTCTGCCGGGAgg                                                              >  2:941057/1‑148 (MQ=255)
                                                                                       ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                          <  1:488522/151‑1 (MQ=255)
                                                                                        gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCC‑GCGGTCCGGCAAAAACGCCCGCGAAGTGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTg                                                          >  1:56653/1‑150 (MQ=255)
                                                                                         agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCgggg                                                         >  1:41820/1‑150 (MQ=255)
                                                                                            tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGCGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga                                                      <  1:618032/150‑1 (MQ=255)
                                                                                                       aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                          <  1:672004/151‑1 (MQ=255)
                                                                                                       aTTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGGAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                          <  2:19337/151‑1 (MQ=255)
                                                                                                         tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATa                                         <  2:129321/150‑1 (MQ=255)
                                                                                                             gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGAt                                     <  1:319767/150‑1 (MQ=255)
                                                                                                             gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                    <  1:920830/151‑1 (MQ=255)
                                                                                                                  aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                                 <  1:117919/149‑1 (MQ=255)
                                                                                                                          aaggcgaagTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGTACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAt                       <  2:260474/151‑1 (MQ=255)
                                                                                                                           aggcgaagTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                      <  2:360644/151‑1 (MQ=255)
                                                                                                                               gaagTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc                  <  2:510084/151‑1 (MQ=255)
                                                                                                                                 agTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgc                <  1:82560/151‑1 (MQ=255)
                                                                                                                                   tGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCgcgc              >  1:639868/1‑151 (MQ=255)
                                                                                                                                        ggcgggGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCg         <  2:50058/151‑1 (MQ=255)
                                                                                                                                         gcgggGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCg         >  2:961829/1‑150 (MQ=255)
                                                                                                                                         gcgggGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGa        <  2:1046408/151‑1 (MQ=255)
                                                                                                                                             ggtggtGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAc    <  2:687703/151‑1 (MQ=255)
                                                                                                                                               tggtgATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCt  >  1:145878/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
CGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCT‑GCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCT  >  NC_000913/599433‑599726

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