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

CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599429‑599725
                                                                                                                                                    |                                                                                                                                                     
cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGACCTAACCGGCGAAGGCGAAGTGGC‑GGGCGGggtggt                                                                                                                                                    <  2:539185/151‑1 (MQ=255)
 gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑Ggggggggggggg                                                                                                                                                   >  1:502131/1‑149 (MQ=255)
     gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGCGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGGGGTcc                                                                                                                                               >  1:358411/1‑151 (MQ=255)
      acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCt                                                                                                                                              <  2:380507/151‑1 (MQ=255)
      acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCt                                                                                                                                              <  1:362522/151‑1 (MQ=255)
       cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGCGggggggggtcc                                                                                                                                                >  2:482458/1‑144 (MQ=255)
          cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCCGc                                                                                                                                            >  2:863174/1‑149 (MQ=255)
           tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGACCTGCg                                                                                                                                           >  2:668114/1‑149 (MQ=255)
            ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGAt                                                                                                                                         <  1:669052/150‑1 (MQ=255)
              aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGAGCCTGCggt                                                                                                                                         >  1:581676/1‑146 (MQ=255)
                 cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGc                                                                                                                                    <  1:711792/150‑1 (MQ=255)
                 cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCa                                                                                                                                   <  2:405681/151‑1 (MQ=255)
                   caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGATCCGGCaa                                                                                                                                  >  2:349599/1‑150 (MQ=255)
                    aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGGTTCGGCaaa                                                                                                                                 >  1:685413/1‑150 (MQ=255)
                          ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGcc                                                                                                                           <  2:918622/150‑1 (MQ=255)
                            aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCg                                                                                                                         <  1:597534/150‑1 (MQ=255)
                              aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGATCCCGCAAAAAAGccc                                                                                                                          >  2:40553/1‑147 (MQ=255)
                               aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGaa                                                                                                                     <  2:215066/151‑1 (MQ=255)
                               aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGaa                                                                                                                     <  1:395070/151‑1 (MQ=255)
                                gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAg                                                                                                                    <  2:106602/151‑1 (MQ=255)
                                 ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGt                                                                                                                   <  1:81001/151‑1 (MQ=255)
                                   aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGGTCCCGCAAAAAAGCCCGGGaa                                                                                                                     >  2:659060/1‑147 (MQ=255)
                                    gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGGATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGa                                                                                                                 <  1:673704/150‑1 (MQ=255)
                                    gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCCGCGATCCGGCAAAAACGCCCGAGaa                                                                                                                     >  1:84161/1‑146 (MQ=255)
                                      gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcg                                                                                                              <  2:670855/151‑1 (MQ=255)
                                       aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGCGGGGGGGGGTCCTGGGATCCGGCCAAAAAGCCCGGGAAggggatg                                                                                                              >  2:321494/1‑144 (MQ=255)
                                          aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                                                                                          <  2:680415/151‑1 (MQ=255)
                                            ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGccgccg                                                                                                        <  1:141951/151‑1 (MQ=255)
                                             gCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGGGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGATCCCGCAAAAACGCCCGAGAAGTGAACgccgcc                                                                                                         >  1:642458/1‑149 (MQ=255)
                                              cGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTg                                                                                                      <  1:729255/151‑1 (MQ=255)
                                                  cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                                  <  2:752081/151‑1 (MQ=255)
                                                     gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGa                                                                                                 >  1:845381/1‑149 (MQ=255)
                                                      tttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAgg                                                                                                  >  1:514977/1‑147 (MQ=255)
                                                        tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                             >  1:354259/1‑150 (MQ=255)
                                                           cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGGCAAACTgg                                                                                         >  1:727284/1‑151 (MQ=255)
                                                             gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa                                                                                       <  2:134366/151‑1 (MQ=255)
                                                                cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                     <  1:801890/150‑1 (MQ=255)
                                                                          aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGCGGGGAACCTGCGATTCGGCTAAAAAGCCCGAGAAGTGATGGGCGCCGGGATGGGCAAACTGGAAAcaca                                                                                  >  2:749392/1‑140 (MQ=255)
                                                                                cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATTCTGCGGTCCGGGAAAAACGCCCGAGAAGCGGTCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGTAGTCt                                                                      >  1:327753/1‑149 (MQ=255)
                                                                                 aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGAACCTGCGATCCGGCGAAAACGCCCGGGGAGTGATGGCCGCCGCGAAGGAAAAAAAGGAAACGCCGAAAAGGAGTATGc                                                                    >  2:162780/1‑150 (MQ=255)
                                                                                         ccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa                                                               <  1:186023/147‑1 (MQ=255)
                                                                                          cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGCGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCg                                                           >  2:78568/1‑150 (MQ=255)
                                                                                                   gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGGTCGCGCGGTCCGGCAAAAACGCCCGGGAAGAGATCGCCGGCGTGAAGGACCAAATGGGAACGGTGGAAAAGAGGCTTCCGGAAGGGGGGgaggaagt                                                 >  2:619719/1‑146 (MQ=255)
                                                                                                      ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTac                                              <  2:310648/151‑1 (MQ=255)
                                                                                                      gaaggaTTGCCGATCTAAACGGCGAAGGAGAAGTGGC‑GGGCGGAATGGTGATCCTGCGAACCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGAAGATAGt                                                 <  2:94560/143‑1 (MQ=255)
                                                                                                        ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                            <  2:292221/151‑1 (MQ=255)
                                                                                                          cATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGGGGGGGGGTGGTCCTGCGGTCCGGCAAAAACGCCCGCGAAGTGAATGCCGCCGTGGAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCCGAAAGCGTGGAGATTGTTaaa                                             >  2:729923/1‑146 (MQ=255)
                                                                                                            ttGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGAAAAAATGGAAACGCCGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaaca                                          >  2:636714/1‑149 (MQ=255)
                                                                                                             tGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCTGCGATCCCGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGGAGGGGGGGGGGAAGTTACAACATa                                        >  2:159121/1‑150 (MQ=255)
                                                                                                              gCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGAGATAGTTACAAAATa                                        >  1:425351/1‑149 (MQ=255)
                                                                                                                cGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGGGGTCCCGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGGAAAGCTGGAAAGTAGTCTGCCGGAAGGGGGGGGGATAGTTACAACAAACGa                                     >  2:439144/1‑150 (MQ=255)
                                                                                                                 gcagtcaaCGGCGAAGGCGAGGCGGG‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                       <  2:498178/141‑1 (MQ=255)
                                                                                                                  aaCTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCg                                  <  2:192519/151‑1 (MQ=255)
                                                                                                                  aaCTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGGGAAGGACAAACTTGAAACGCTGAAAAGTAGGCTTCCGGAAGGCGGGGGGATAGTTACAACAaacgaacg                                  >  1:326081/1‑151 (MQ=255)
                                                                                                                     tAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATACGATCGc                                 >  1:746273/1‑149 (MQ=255)
                                                                                                                      aaaCGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                                <  2:231897/149‑1 (MQ=255)
                                                                                                                               aggcgaagTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                     <  2:716855/151‑1 (MQ=255)
                                                                                                                               aggcgaagTGGC‑GGCCGGGGTGGTGATCATGCGAACCGGCAAAAACGCCCATGAAGTGATCGCCGCCATGAAGGTCAAACTGGAAACGCTGAAAAGTAGTATGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCa                       <  2:708138/149‑1 (MQ=255)
                                                                                                                                    aagTGGC‑GGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcg                >  2:70168/1‑151 (MQ=255)
                                                                                                                                       tGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCgcgc             >  1:485292/1‑151 (MQ=255)
                                                                                                                                          c‑gggcgggGTGGTGCTCCTGCGATCCGGAAAAAACGCCCGAGAAGTGACCGCCGCCGTGAAGGACAACCTGGAAACGCTGAAAAGTAGGCTGCCGGAAGGCGTGGAGATAGTTAAAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt          <  1:121485/151‑1 (MQ=255)
                                                                                                                                            gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc         <  1:501040/151‑1 (MQ=255)
                                                                                                                                            gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGACAt          >  1:290948/1‑150 (MQ=255)
                                                                                                                                               cgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc      >  2:901908/1‑151 (MQ=255)
                                                                                                                                                   gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATAGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCAGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTAATTGACCGCGCCATCGACAAcc  >  2:69702/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGC‑GGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599429‑599725

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