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

CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599429‑599725
                                                                                                                                                    |                                                                                                                                                     
cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGACCTAACCGGCGAAGGCGAAGTGGCGGGCGGggtggt                                                                                                                                                    <  2:539185/151‑1 (MQ=255)
 gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggggg                                                                                                                                                   >  1:502131/1‑149 (MQ=255)
     gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGCGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTcc                                                                                                                                               >  1:358411/1‑151 (MQ=255)
      acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                              <  2:380507/151‑1 (MQ=255)
      acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt                                                                                                                                              <  1:362522/151‑1 (MQ=255)
       cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGCGggggggggtcc                                                                                                                                                >  2:482458/1‑144 (MQ=255)
          cTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGc                                                                                                                                            >  2:863174/1‑149 (MQ=255)
           tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGACCTGCg                                                                                                                                           >  2:668114/1‑149 (MQ=255)
            ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAt                                                                                                                                         <  1:669052/150‑1 (MQ=255)
              aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGCCTGCggt                                                                                                                                         >  1:581676/1‑146 (MQ=255)
                 cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGc                                                                                                                                    <  1:711792/150‑1 (MQ=255)
                 cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCa                                                                                                                                   <  2:405681/151‑1 (MQ=255)
                   caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCaa                                                                                                                                  >  2:349599/1‑150 (MQ=255)
                    aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTTCGGCaaa                                                                                                                                 >  1:685413/1‑150 (MQ=255)
                          ttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGcc                                                                                                                           <  2:918622/150‑1 (MQ=255)
                            aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCg                                                                                                                         <  1:597534/150‑1 (MQ=255)
                              aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCCGCAAAAAAGccc                                                                                                                          >  2:40553/1‑147 (MQ=255)
                               aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGaa                                                                                                                     <  2:215066/151‑1 (MQ=255)
                               aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGaa                                                                                                                     <  1:395070/151‑1 (MQ=255)
                                gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAg                                                                                                                    <  2:106602/151‑1 (MQ=255)
                                 ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGt                                                                                                                   <  1:81001/151‑1 (MQ=255)
                                   aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCCGCAAAAAAGCCCGGGaa                                                                                                                     >  2:659060/1‑147 (MQ=255)
                                    gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGGATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGa                                                                                                                 <  1:673704/150‑1 (MQ=255)
                                    gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCGCGATCCGGCAAAAACGCCCGAGaa                                                                                                                     >  1:84161/1‑146 (MQ=255)
                                      gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcg                                                                                                              <  2:670855/151‑1 (MQ=255)
                                       aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGCGGGGGGGGGTCCTGGGATCCGGCCAAAAAGCCCGGGAAggggatg                                                                                                              >  2:321494/1‑144 (MQ=255)
                                          aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                                                                                          <  2:680415/151‑1 (MQ=255)
                                            ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGccgccg                                                                                                        <  1:141951/151‑1 (MQ=255)
                                             gCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGGGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCCGCAAAAACGCCCGAGAAGTGAACgccgcc                                                                                                         >  1:642458/1‑149 (MQ=255)
                                              cGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTg                                                                                                      <  1:729255/151‑1 (MQ=255)
                                                  cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                                  <  2:752081/151‑1 (MQ=255)
                                                     gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGa                                                                                                 >  1:845381/1‑149 (MQ=255)
                                                      tttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAgg                                                                                                  >  1:514977/1‑147 (MQ=255)
                                                        tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                             >  1:354259/1‑150 (MQ=255)
                                                           cTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGGCAAACTgg                                                                                         >  1:727284/1‑151 (MQ=255)
                                                             gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa                                                                                       <  2:134366/151‑1 (MQ=255)
                                                                cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                                                                                     <  1:801890/150‑1 (MQ=255)
                                                                          aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGCGGGGAACCTGCGATTCGGCTAAAAAGCCCGAGAAGTGATGGGCGCCGGGATGGGCAAACTGGAAAcaca                                                                                  >  2:749392/1‑140 (MQ=255)
                                                                                cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTGCGGTCCGGGAAAAACGCCCGAGAAGCGGTCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGTAGTCt                                                                      >  1:327753/1‑149 (MQ=255)
                                                                                 aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAACCTGCGATCCGGCGAAAACGCCCGGGGAGTGATGGCCGCCGCGAAGGAAAAAAAGGAAACGCCGAAAAGGAGTATGc                                                                    >  2:162780/1‑150 (MQ=255)
                                                                                         ccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa                                                               <  1:186023/147‑1 (MQ=255)
                                                                                          cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGCGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCg                                                           >  2:78568/1‑150 (MQ=255)
                                                                                                   gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGGTCGCGCGGTCCGGCAAAAACGCCCGGGAAGAGATCGCCGGCGTGAAGGACCAAATGGGAACGGTGGAAAAGAGGCTTCCGGAAGGGGGGgaggaagt                                                 >  2:619719/1‑146 (MQ=255)
                                                                                                      ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTac                                              <  2:310648/151‑1 (MQ=255)
                                                                                                      gaaggaTTGCCGATCTAAACGGCGAAGGAGAAGTGGCGGGCGGAATGGTGATCCTGCGAACCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGAAGATAGt                                                 <  2:94560/143‑1 (MQ=255)
                                                                                                        ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                            <  2:292221/151‑1 (MQ=255)
                                                                                                          cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGTCCTGCGGTCCGGCAAAAACGCCCGCGAAGTGAATGCCGCCGTGGAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCCGAAAGCGTGGAGATTGTTaaa                                             >  2:729923/1‑146 (MQ=255)
                                                                                                            ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGAAAAAATGGAAACGCCGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaaca                                          >  2:636714/1‑149 (MQ=255)
                                                                                                             tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCCGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGGAGGGGGGGGGGAAGTTACAACATa                                        >  2:159121/1‑150 (MQ=255)
                                                                                                              gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGAGATAGTTACAAAATa                                        >  1:425351/1‑149 (MQ=255)
                                                                                                                cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGGAAAGCTGGAAAGTAGTCTGCCGGAAGGGGGGGGGATAGTTACAACAAACGa                                     >  2:439144/1‑150 (MQ=255)
                                                                                                                 gcagtcaaCGGCGAAGGCGAGGCGGGGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                       <  2:498178/141‑1 (MQ=255)
                                                                                                                  aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCg                                  <  2:192519/151‑1 (MQ=255)
                                                                                                                  aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGGGAAGGACAAACTTGAAACGCTGAAAAGTAGGCTTCCGGAAGGCGGGGGGATAGTTACAACAaacgaacg                                  >  1:326081/1‑151 (MQ=255)
                                                                                                                     tAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATACGATCGc                                 >  1:746273/1‑149 (MQ=255)
                                                                                                                      aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                                <  2:231897/149‑1 (MQ=255)
                                                                                                                               aggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                     <  2:716855/151‑1 (MQ=255)
                                                                                                                               aggcgaagTGGCGGCCGGGGTGGTGATCATGCGAACCGGCAAAAACGCCCATGAAGTGATCGCCGCCATGAAGGTCAAACTGGAAACGCTGAAAAGTAGTATGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCa                       <  2:708138/149‑1 (MQ=255)
                                                                                                                                    aagTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcg                >  2:70168/1‑151 (MQ=255)
                                                                                                                                       tGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCgcgc             >  1:485292/1‑151 (MQ=255)
                                                                                                                                          cgggcgggGTGGTGCTCCTGCGATCCGGAAAAAACGCCCGAGAAGTGACCGCCGCCGTGAAGGACAACCTGGAAACGCTGAAAAGTAGGCTGCCGGAAGGCGTGGAGATAGTTAAAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt          <  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)
                                                                                                                                                    |                                                                                                                                                     
CGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  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: