Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5790TG44.7% 127.0 / 14.0 77V289G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (7/35);  new base G (34/0);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 6.13e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCA  >  NC_000913/599432‑599728
                                                                                                                                                   |                                                                                                                                                      
tCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGggggtgg                                                                                                                                                      >  2:914082/1‑149 (MQ=255)
  gccgccTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAAGTCCAGATTAGCCCGGAGACGCGCCGGGGCGTCGCCGAACTACACGGCCAAGGCGCCGTCGCGGGTGGGggggtgggg                                                                                                                                                    >  1:60110/6‑147 (MQ=255)
     gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGAGcc                                                                                                                                                   >  1:718643/1‑148 (MQ=255)
      acTTTAATCACATCGTTTTAGAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATcc                                                                                                                                                   >  2:505961/1‑147 (MQ=255)
         ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGACCTGccg                                                                                                                                              >  1:667443/1‑147 (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:1282524/151‑1 (MQ=255)
                          aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgag                                                                                                                           <  2:547299/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:821370/151‑1 (MQ=255)
                                            gTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa                                                                                                         <  2:540449/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                                                                                                  <  2:33211/150‑1 (MQ=255)
                                                    ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaa                                                                                                  <  1:71575/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)
                                                                                                                                                  gtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCa  <  1:1097553/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTCA  >  NC_000913/599432‑599728

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