Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG50.6% 99.1 / 142.5 77V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (2/36);  new base G (39/0);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 6.03e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-14

CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAC  >  NC_000913/599429‑599724
                                                                                                                                                   |                                                                                                                                                     
cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggg                                                                                                                                                    >  1:282324/1‑150 (MQ=255)
cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGggggggg                                                                                                                                                    >  1:728604/1‑150 (MQ=255)
 gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg                                                                                                                                                   >  2:909933/1‑149 (MQ=255)
   tCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGggggtgg                                                                                                                                                  >  2:914082/1‑149 (MQ=255)
     gccgccTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAAGTCCAGATTAGCCCGGAGACGCGCCGGGGCGTCGCCGAACTACACGGCCAAGGCGCCGTCGCGGGTGGGggggtgggg                                                                                                                                                >  1:60110/6‑147 (MQ=255)
       cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggggtc                                                                                                                                                >  1:226274/1‑143 (MQ=255)
        gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGAGcc                                                                                                                                               >  1:718643/1‑148 (MQ=255)
         acTTTAATCACATCGTTTTAGAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATcc                                                                                                                                               >  2:505961/1‑147 (MQ=255)
            ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGACCTGccg                                                                                                                                          >  1:667443/1‑147 (MQ=255)
             tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggggtctctcgaccc                                                                                                                                       >  2:1112616/1‑137 (MQ=255)
             tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATc                                                                                                                                        <  2:802052/150‑1 (MQ=255)
             tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTcccacgaacc                                                                                                                                       >  2:292000/1‑143 (MQ=255)
              aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATcc                                                                                                                                       >  1:472905/1‑150 (MQ=255)
                 cacaTCGATTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCGCGGTGGTGCGCCGGGGCCTTGCCGAACTAAACGGCGAAGGCGAAGTGGAGGGCGGGGGGGTGGTcccccgcac                                                                                                                                        >  1:184375/1‑139 (MQ=255)
                  acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGATcc                                                                                                                                       >  2:115060/1‑146 (MQ=255)
                    aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaaaa                                                                                                                                <  1:169559/151‑1 (MQ=255)
                     tCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCaaaa                                                                                                                                >  2:608076/1‑150 (MQ=255)
                           tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGCCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGGCTGCGGACCGGCaaaaa                                                                                                                               >  1:1061890/1‑145 (MQ=255)
                             aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgag                                                                                                                       <  2:547299/151‑1 (MQ=255)
                             aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgag                                                                                                                       <  2:1282524/151‑1 (MQ=255)
                             aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTTCGATCCGGCCAAAAAGccc                                                                                                                          >  2:222096/1‑148 (MQ=255)
                              aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCCGTCCCGCAAAAAcccc                                                                                                                           >  1:924922/1‑146 (MQ=255)
                                gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCCGGGAGCCGGGGAAAAAGcccca                                                                                                                        >  1:256060/1‑145 (MQ=255)
                                 ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGt                                                                                                                   <  1:1107445/151‑1 (MQ=255)
                                    gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGAt                                                                                                                <  1:509398/151‑1 (MQ=255)
                                     tcaaaaTGGCATGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGTGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGAt                                                                                                                <  2:222151/148‑1 (MQ=255)
                                     tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                               >  2:1236700/1‑151 (MQ=255)
                                       aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgc                                                                                                             <  1:276740/151‑1 (MQ=255)
                                         aaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCTGGGATTCGGCCAAAAAGccccgggagaggtcgcc                                                                                                            >  2:511594/1‑136 (MQ=255)
                                             gCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGt                                                                                                       <  1:909892/151‑1 (MQ=255)
                                               gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa                                                                                                     <  2:540449/151‑1 (MQ=255)
                                               gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa                                                                                                     <  2:821370/151‑1 (MQ=255)
                                                 gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGt                                                                                                       >  2:1263433/1‑147 (MQ=255)
                                                   ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGa                                                                                                 <  2:1156890/151‑1 (MQ=255)
                                                   ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCTGGGTCCGGCCAAAACGCCCGGGAAGGGGTCGCCGCCGTGGAggg                                                                                                 >  1:283471/1‑150 (MQ=255)
                                                      tttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACa                                                                                               <  1:91614/150‑1 (MQ=255)
                                                      tttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGCTCCGGCAAAAAACCCCGAGGAGGGATCGCCGCCGTGGAgg                                                                                                  >  1:218400/1‑147 (MQ=255)
                                                       ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaa                                                                                              <  1:71575/150‑1 (MQ=255)
                                                       ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaa                                                                                              <  2:33211/150‑1 (MQ=255)
                                                        tgtCTGCGCGATGTTGCAAAGGTCCAGATTGGCCCGGGGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGGGAAGTGCCGGGCGGGGTGGCGATCCTGAGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCAGTGAAGGACaaa                                                                                             <  1:792945/148‑1 (MQ=255)
                                                        tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGTTGATCCTGCGATCCGGAAAAATCGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                             <  2:1106259/150‑1 (MQ=255)
                                                         aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGTGCCCTCGATCCGGCAAAAAAGCCCGCGAAGGGATCGCCGCCGGGAAGGGCAAACt                                                                                           >  1:1080564/1‑151 (MQ=255)
                                                         aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGGGCTCCGGCAAAAAACCCCGAGAAGTGAGCGCCGCCGGGGAGGACaaa                                                                                             >  1:215309/1‑149 (MQ=255)
                                                            tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGGGAAGGGCaaa                                                                                             >  2:854521/1‑146 (MQ=255)
                                                             gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGTTATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                           <  2:1007571/147‑1 (MQ=255)
                                                             gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTg                                                                                          <  1:257004/148‑1 (MQ=255)
                                                               gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGGTTCCGGAAAAACGCCCGAGAAGTGATCGCCGCCGTGCAGGAGAAACTGGaaa                                                                                      >  1:126207/1‑150 (MQ=255)
                                                               gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                     <  1:463868/151‑1 (MQ=255)
                                                                   tGTTGCGAAGGTCCAGATGGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                                 <  1:1074966/150‑1 (MQ=255)
                                                                    gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGa                                                                               <  2:241561/151‑1 (MQ=255)
                                                                      tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTCCGAACCGGCGAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACAGGGAAAGCTGaa                                                                              >  2:661391/1‑150 (MQ=255)
                                                                        cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGAGGACGCCGCCGGGAAGGAGAAAAAGGAAAAGCTGgaa                                                                             >  1:776088/1‑149 (MQ=255)
                                                                        cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGGAGGACAAACTGGAAACGCTGaaa                                                                             >  2:1129186/1‑149 (MQ=255)
                                                                         gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCCGCAAAAACGCCCCAGAAATGGTCGCCGCCGTGAAGGACAAACTGGAAAAGCtcgaaa                                                                            >  1:1214669/1‑144 (MQ=255)
                                                                         gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGCAAAAAAGCCCGAGAAGTGAACGCCGCCGGGAAGGACAAAATGGAAAAGCTGgaa                                                                             >  1:1165062/1‑148 (MQ=255)
                                                                            ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                        <  1:534783/150‑1 (MQ=255)
                                                                               ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGCAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCCAAAACGCCCCGGGAGGGATTGCCGCCGGCGAGGGCCGACGAGAAACGCTGAGAACGAGTCt                                                                     >  1:437352/1‑150 (MQ=255)
                                                                               ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTTCTGCCGTCCCGGGAAAAAACCCGGGAAAGGATTGCCGCCGGGGAAGACAAAATGGGAAAAGTGaaaaagagagt                                                                     >  2:345378/1‑143 (MQ=255)
                                                                               ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCCAAAAAGCCCCCGGAGGGATCGCCGCCCGGGAGGACAAACTGGAAACGCTTGaaaata                                                                          >  1:625232/1‑144 (MQ=255)
                                                                                cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTg                                                                    <  2:283996/150‑1 (MQ=255)
                                                                                     tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGATCCGGCAAAAACACCCGGGCAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGcc                                                                  >  1:218371/1‑147 (MQ=255)
                                                                                      ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa                                                              <  2:226278/150‑1 (MQ=255)
                                                                                        cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                           <  1:855158/151‑1 (MQ=255)
                                                                                          cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGACCGCGCTCCGGCAAAAAAGCCCCAGAAGGGATCCGCGCCGTGAGGGGCAAACTGGAAAAACTGCAAAGTAGTTTTTCGGaa                                                              >  1:1054194/1‑146 (MQ=255)
                                                                                           ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                        <  1:209932/151‑1 (MQ=255)
                                                                                             agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGGAGGACAAACTGGGAAAGCTGAAAAGGAGTCTGCCGGAAGGCgggg                                                       >  2:658312/1‑150 (MQ=255)
                                                                                              gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag                                                     <  2:472934/151‑1 (MQ=255)
                                                                                                     cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGTGAACCGGCAAAAAAGCCCCAGAAGTGAGCGCAGCCGGGCAGGACAAAATAGGAACGGCGGAAAGGAGGCTGCCTGAAGGCGGGGAGATa                                                  >  1:551936/1‑147 (MQ=255)
                                                                                                         gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCCAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGGAAAGAAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                           >  1:571038/1‑150 (MQ=255)
                                                                                                           aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                        <  2:246017/151‑1 (MQ=255)
                                                                                                              gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGTGGGGGGGGGGTTCTCCCCTCCCGCAAAAAAAACCGAGAAGTGGTCGGCGCCGTGTAGGGCAAACCGGAAACGACGAAAAAAAGACTACCGGAAGGCGGGGAAGTAGGTacaa                                           >  2:1111255/1‑145 (MQ=255)
                                                                                                              gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                      <  2:604234/150‑1 (MQ=255)
                                                                                                              gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACg                                     <  1:881883/151‑1 (MQ=255)
                                                                                                               ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                      <  2:551967/149‑1 (MQ=255)
                                                                                                                   aCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAATTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                  <  1:200876/149‑1 (MQ=255)
                                                                                                                   aCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGc                                <  2:1130667/151‑1 (MQ=255)
                                                                                                                    cTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGGAGTGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGGAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGc                                >  2:1166474/1‑150 (MQ=255)
                                                                                                                       aaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcc                            <  2:893287/151‑1 (MQ=255)
                                                                                                                         cGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGCTCCGGCAAAAAAGCCCGAGGAGGGATCGCCGCCGGGGAGGACAAACTGGAAACGCTGGAAAGGAGTCTGCCCGGAGGGGGGGGGGAAGTTTAAACAAACGATCGcagcc                            >  2:1046995/1‑149 (MQ=255)
                                                                                                                            cgaaggcgaagTGGCGGGCGGGGGGGTGATCCTGCGAACCGGCAAAAACCCCCGAGAAGGGATCGCCGCCGTGGAGGACAAACTGGAAAAGCTGAAAAATAGTCTGCCGGAAGGCGGGGGGATAGGTACAACATACGATCGCAGCCAGCTc                       >  2:1159127/1‑151 (MQ=255)
                                                                                                                             gaaggcgaagTGGCGGGCGGGGGGGTGATCCTGCGATCCCGCCAAAACGCCCGAGAAGTGGTCCCCGCCGTGAAGGAAAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGGGGGGGGAGTTACAAAATACGATCGCAGCCAGCTc                       >  1:1158177/1‑150 (MQ=255)
                                                                                                                                          cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt         >  1:904905/1‑151 (MQ=255)
                                                                                                                                           gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc        <  2:107627/151‑1 (MQ=255)
                                                                                                                                               ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc     >  2:368764/1‑150 (MQ=255)
                                                                                                                                                 ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAc  >  1:1212131/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAC  >  NC_000913/599429‑599724

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: