Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG41.8% 220.3 / 49.8 122V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (5/66);  new base G (51/0);  total (56/66)
Fisher's exact test for biased strand distribution p-value = 5.10e-29
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599426‑599722
                                                                                                                                                       |                                                                                                                                                  
aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGggt                                                                                                                                                    <  1:1534518/151‑1 (MQ=255)
   cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCTCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGgggggg                                                                                                                                                  >  4:1528961/1‑150 (MQ=255)
      tCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGgggggggg                                                                                                                                                >  4:1557186/1‑147 (MQ=255)
       cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGgggggggggt                                                                                                                                              >  2:527824/1‑146 (MQ=255)
       cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGgggggggg                                                                                                                                                >  1:835334/1‑146 (MQ=255)
       cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGa                                                                                                                                               <  1:4955/149‑1 (MQ=255)
        gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑Ggggggggggggggttc                                                                                                                                            >  2:326857/1‑145 (MQ=255)
        gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTCCCCGTTTATCTGCGCGATGTTGCGAACGTCCAGAATGGCCCGGAGATGCGCCGGGGCACTGCCGCACTAACCGGCGAAGCCGAAGTGGC‑GGGGGGggagggggtc                                                                                                                                             >  1:1427862/1‑145 (MQ=255)
          cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTg                                                                                                                                          <  1:1058550/151‑1 (MQ=255)
            acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGGTCCt                                                                                                                                           >  3:149077/1‑148 (MQ=255)
            acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCg                                                                                                                                        <  4:397209/151‑1 (MQ=255)
            acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCg                                                                                                                                        >  4:1069860/1‑151 (MQ=255)
               ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGGTCCTGCggt                                                                                                                                      >  1:1133476/1‑148 (MQ=255)
               ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTg                                                                                                                                          >  1:960679/1‑146 (MQ=255)
                tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATc                                                                                                                                     <  4:1222914/150‑1 (MQ=255)
                tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATc                                                                                                                                     <  1:691955/150‑1 (MQ=255)
                 aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGGTCCg                                                                                                                                   >  4:126460/1‑151 (MQ=255)
                    cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCa                                                                                                                                <  4:395663/151‑1 (MQ=255)
                    cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGagacgg                                                                                                                                  >  4:1081144/1‑144 (MQ=255)
                     acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTccatcggtac                                                                                                                                    >  3:495371/1‑138 (MQ=255)
                       aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGGTCCTTCGGTCCGGCaaa                                                                                                                              >  4:1116475/1‑150 (MQ=255)
                          gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAAc                                                                                                                           <  3:383898/150‑1 (MQ=255)
                          gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACg                                                                                                                          <  4:1787120/151‑1 (MQ=255)
                           ttttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGGTCCTGCGGTCCGGCaaaaa                                                                                                                            >  1:7299/1‑148 (MQ=255)
                            tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGGTCCCGCGGTCCCGCcaaaa                                                                                                                            >  1:21791/1‑147 (MQ=255)
                              tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGGTCCGGCAAAAAAGccc                                                                                                                       >  4:1030563/1‑150 (MQ=255)
                              tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGccc                                                                                                                       >  1:1627895/1‑150 (MQ=255)
                              tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGccc                                                                                                                       >  4:279299/1‑150 (MQ=255)
                               aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCga                                                                                                                     <  1:1739169/151‑1 (MQ=255)
                                aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgag                                                                                                                    <  1:581012/151‑1 (MQ=255)
                                 aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgaga                                                                                                                   <  3:139152/151‑1 (MQ=255)
                                  aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGGTCCGGCAAAAACGccc                                                                                                                       >  1:374442/1‑146 (MQ=255)
                                  aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGGCCTGCGATCCGGCAAAAAAGccc                                                                                                                       >  1:909269/1‑146 (MQ=255)
                                    ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCCGCGGTCCGGCAAAAACGCCCCGGAAgg                                                                                                                >  3:379274/1‑150 (MQ=255)
                                     cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAggg                                                                                                               >  2:514245/1‑149 (MQ=255)
                                      aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGGTCCGGCCAAAACGCCCCAGAAggg                                                                                                               >  1:462050/1‑148 (MQ=255)
                                           aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGTGGTGGTCCTGCGGGACGGAAAAAAAGCCAGAGAAgag                                                                                                               >  1:1663684/1‑143 (MQ=255)
                                           aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGGGAAGTGATcgcc                                                                                                         >  4:210763/1‑151 (MQ=255)
                                             aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                                                                                       <  3:527803/151‑1 (MQ=255)
                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                      <  1:1320936/150‑1 (MQ=255)
                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGAATCTGGGAACCGGCAAAAACGCCCGAGAAGTGGTCgccccc                                                                                                      >  4:346709/1‑150 (MQ=255)
                                                gCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGt                                                                                                    <  4:1377495/151‑1 (MQ=255)
                                                 tgTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGACTCGGCAAAAAAACCCGAGAAGTGATCGccgccggg                                                                                                   >  1:89466/2‑149 (MQ=255)
                                                  gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTg                                                                                                   >  2:206430/1‑150 (MQ=255)
                                                    gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAg                                                                                                <  4:705662/151‑1 (MQ=255)
                                                    gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAg                                                                                                <  4:870530/151‑1 (MQ=255)
                                                     cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                               <  2:768/151‑1 (MQ=255)
                                                      ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGAACCTGGGAACCAGCGAAAAAGCCCGAGAAGTGATCGccgccggg                                                                                                   >  4:1715759/1‑144 (MQ=255)
                                                        gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACa                                                                                            <  4:870310/151‑1 (MQ=255)
                                                                gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTgg                                                                                      <  4:1028677/149‑1 (MQ=255)
                                                                 cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTACCGCGATCCGGCAAAAACGCCCGAGGAGGGATCGCCGCCGTGGAGGACAAACTGGaa                                                                                    >  3:88663/1‑150 (MQ=255)
                                                                  gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                   <  1:1378825/150‑1 (MQ=255)
                                                                   cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                  >  4:1773324/1‑150 (MQ=255)
                                                                    gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAAATGGgaaag                                                                                 >  2:535799/1‑148 (MQ=255)
                                                                    gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTTCTGCGATCCGGCAAAAACGCCCCGGGCACGGGCGCCGCCGGGAAGGGAAAACTGGAAAc                                                                                  >  3:305611/1‑149 (MQ=255)
                                                                     aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                               <  4:1079542/151‑1 (MQ=255)
                                                                     aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                               <  4:1534463/151‑1 (MQ=255)
                                                                       gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGTAAACGCTGa                                                                             <  1:1836844/151‑1 (MQ=255)
                                                                       gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGa                                                                             <  4:1089445/151‑1 (MQ=255)
                                                                         tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTTCTGCGGTCCGGCAAAAAAGCCCGGGAAGGGGTCGGCGCCGGGGAGGACAAAACGGGAACACTgga                                                                            >  1:326899/1‑148 (MQ=255)
                                                                          gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa                                                                           <  4:565059/150‑1 (MQ=255)
                                                                             aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagta                                                                       <  1:534190/151‑1 (MQ=255)
                                                                             aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGaaaa                                                                          >  4:1321870/1‑148 (MQ=255)
                                                                             aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGCGGAGGGCAAACCGGAAAAGCCGgaaa                                                                          >  1:1912150/1‑148 (MQ=255)
                                                                               ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                      <  1:1604099/150‑1 (MQ=255)
                                                                                gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc                                                                    <  1:60374/151‑1 (MQ=255)
                                                                                      aTTGGCCCGGAGATGCGCCGGGGCATTGCCGCACTAAACGGCGAAGGCGCAGTGGCGGGGGGGGGGGGTGA‑CCCCCGAATCGGGCAAAACACCCGACGAAAGCTTGCCGCCCCGAAGGAGAAACTTGAAACGCTGCAAAAGAGTCCGcc                                                                >  1:438122/1‑149 (MQ=255)
                                                                                        tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGCGAAGTGATCGCCGCCGGGAAGGGCAAACTGGAAAAGCTGGAAAGGAGTCTGCCgg                                                              >  1:954769/1‑149 (MQ=255)
                                                                                         ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa                                                            <  4:1047146/150‑1 (MQ=255)
                                                                                          gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                          <  2:281980/151‑1 (MQ=255)
                                                                                           cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                          <  4:40726/150‑1 (MQ=255)
                                                                                           cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                         <  1:985646/151‑1 (MQ=255)
                                                                                           cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                         <  3:427985/151‑1 (MQ=255)
                                                                                              ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                      <  4:551581/151‑1 (MQ=255)
                                                                                                 gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                    <  1:1684896/150‑1 (MQ=255)
                                                                                                  aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGGGTTCCTGCGGTCCGGGCAAAACGCCCGGGAAGTGGTCGCCGCCGGGAAAGAGGAAACGGGAACGCTGAAAAGTAGACTCCCGGAAGGCGTggcg                                                   >  4:1034983/1‑148 (MQ=255)
                                                                                                    gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                 <  4:835693/150‑1 (MQ=255)
                                                                                                    gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                <  4:970783/151‑1 (MQ=255)
                                                                                                    gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                <  4:1003823/151‑1 (MQ=255)
                                                                                                      gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGAGCCGGCAAAAACGCCCGAGGAGGGAATGCCGCCGGGAAGGACAAAATGGAAAAGCCGAAAAGAAATATGCCGGGAGGCggggggggag                                               >  1:1239780/1‑144 (MQ=255)
                                                                                                      gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGGAGGACAAACTGGAAACGCTGAAAAGAAGTCTGCCGGAAgggggggggaaa                                                >  1:343652/1‑139 (MQ=255)
                                                                                                       ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGtt                                             <  4:1099513/151‑1 (MQ=255)
                                                                                                        cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTa                                            <  2:164487/151‑1 (MQ=255)
                                                                                                              aTTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaaca                                       <  1:1528360/150‑1 (MQ=255)
                                                                                                              aTTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                      <  4:462221/151‑1 (MQ=255)
                                                                                                                tGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATa                                     <  4:1812712/150‑1 (MQ=255)
                                                                                                                 gCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                    <  2:135328/150‑1 (MQ=255)
                                                                                                                        tAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGAAGGGGGGGGGATAGTTACAACATACGATCg                               >  1:629729/1‑148 (MQ=255)
                                                                                                                         aaaCGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGATCCCGCGAACCGGGAAAAACGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAACTGGGAAAGCCGGAAAGGAGTCTGCCGGGAGGGGGGGGGATAGTTAAAAAACACGAACGc                              >  1:347046/1‑148 (MQ=255)
                                                                                                                          aaCGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGATACGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGGAAAGTAGTCTGCCGGAAGGCGGGGGGGTAGTTACAACATACGATc                                >  1:1101368/1‑145 (MQ=255)
                                                                                                                            cGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcca                         <  1:1080680/150‑1 (MQ=255)
                                                                                                                              gcgaaggcgaagTGGC‑GGGCGGGGGGGGGATCCTGTCATCCCGCCAAAAAGCCCGCGAAATGGTGGCCGCCCCGGAGGGCAAACCGGAAACGCTGAAAAATAGTGCGCCCGGAGGCGTGAAGGTAGTTTCAACATACGATCGcacgc                          >  1:296663/1‑144 (MQ=255)
                                                                                                                               cgaaggcgaagTGGC‑GGGCGGGGTGGGGATCCTGGGATCCGGCAAAAACCCCCGAGAAGAGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGCAGTCTGCCGGAAGGCGGGGGGAAAGTTACAACATACGATCGcagccagccc                     >  4:1410756/1‑149 (MQ=255)
                                                                                                                                     cgaagTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAc               <  2:445941/151‑1 (MQ=255)
                                                                                                                                      gaagTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc              <  4:241191/151‑1 (MQ=255)
                                                                                                                                       aagTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc              <  3:514225/150‑1 (MQ=255)
                                                                                                                                            gc‑gggcggGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCa        <  1:1107086/151‑1 (MQ=255)
                                                                                                                                               gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc      <  4:764578/151‑1 (MQ=255)
                                                                                                                                                 gcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGa    <  2:229767/151‑1 (MQ=255)
                                                                                                                                                    gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa  >  1:1180103/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599426‑599722

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