Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG52.5% 72.6 / 105.2 60V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (2/26);  new base G (31/0);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 9.54e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-12

GCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599430‑599725
                                                                                                                                                  |                                                                                                                                                       
gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGgtggtg                                                                                                                                                     <  2:932077/151‑1 (MQ=255)
        acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCg                                                                                                                                             <  1:959565/151‑1 (MQ=255)
        acTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCt                                                                                                                                                >  1:785923/1‑148 (MQ=255)
          tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAt                                                                                                                                           >  2:448185/1‑151 (MQ=255)
            tAATCACATCGTTTTAAAAGCCAGTGAAACTGGCGTGCCCGTTTATCTGCGCGCTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTCAACGGCGAAGGCGAAGTGGCGGGCGgggggggggtc                                                                                                                                                  >  2:337009/1‑137 (MQ=255)
            tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGATcc                                                                                                                                         >  2:693949/1‑151 (MQ=255)
             aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATcc                                                                                                                                         >  1:479006/1‑150 (MQ=255)
              aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCgg                                                                                                                                       >  1:838547/1‑151 (MQ=255)
              aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGCCCTCGATcc                                                                                                                                         >  2:749843/1‑149 (MQ=255)
                cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGc                                                                                                                                      >  1:810989/1‑150 (MQ=255)
                  caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaaa                                                                                                                                   <  1:616748/151‑1 (MQ=255)
                     cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaaaaa                                                                                                                                 <  2:543761/150‑1 (MQ=255)
                     cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCaaaaa                                                                                                                                 >  1:321499/1‑150 (MQ=255)
                             aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCgcg                                                                                                                         >  1:928645/1‑148 (MQ=255)
                              aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGGTCCGGCAAAAACGCCCcag                                                                                                                         >  1:190879/1‑149 (MQ=255)
                              aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGaa                                                                                                                       <  1:15060/151‑1 (MQ=255)
                                   gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGCGACCCCGCGCTCCCGCAAAAAAGccccggaaatg                                                                                                                    >  1:109493/1‑141 (MQ=255)
                                    tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGGGGCGGGGGGGGGGGGGGTCCCGCGATCCGGCAAAAACGCCCGGGGAGt                                                                                                                     >  1:162948/1‑147 (MQ=255)
                                     gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                                 <  2:550152/150‑1 (MQ=255)
                                     gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                                 <  1:849479/150‑1 (MQ=255)
                                        aaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccg                                                                                                             <  2:45648/151‑1 (MQ=255)
                                         aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGTCTGCGGGCCGGGAAAAACGCCCGGGAAGGGATggccgc                                                                                                            >  2:208884/1‑151 (MQ=255)
                                         aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                                                                                            <  2:191636/151‑1 (MQ=255)
                                                   cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGGGAAGGGATCGCCGCCGTGAAgg                                                                                                    >  1:874177/1‑149 (MQ=255)
                                                      ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGGGCAGTGATCGCCGCCGGGAAGGACaa                                                                                               >  1:10193/1‑150 (MQ=255)
                                                       tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                              <  1:480277/150‑1 (MQ=255)
                                                            gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                            <  2:861597/147‑1 (MQ=255)
                                                            gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa                                                                                       <  2:92789/151‑1 (MQ=255)
                                                            gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa                                                                                       <  2:601573/151‑1 (MQ=255)
                                                            gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGa                                                                                        <  1:385265/150‑1 (MQ=255)
                                                              gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAACCCGCAAAAACACCCGAGAAGGGATCGCCGCCGGGAAGGACAAAATGGgaa                                                                                      >  2:827387/1‑150 (MQ=255)
                                                                     tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGATCCGGGAAAAAACCCCGAGAAGGGGGCACCGGCGGGGAAGACAAACTGGGAAAGCTg                                                                                 >  2:756247/1‑148 (MQ=255)
                                                                       cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGTACTGCGATCCGGCAAAAAAGCCCCAGAAGAGGTCGCCGCCGGGAAGGACAAACTCGGaaaacttgaa                                                                               >  2:461338/1‑142 (MQ=255)
                                                                          aGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag                                                                         <  1:154132/151‑1 (MQ=255)
                                                                           ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtagt                                                                        <  2:842087/151‑1 (MQ=255)
                                                                            gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc                                                                       <  1:300239/151‑1 (MQ=255)
                                                                             tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGGTCGCCGCCGGGGAGGGCAAACCGGAAAAGCTGAAAAGGAGTc                                                                       >  2:168657/1‑150 (MQ=255)
                                                                               cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTg                                                                     <  1:718930/150‑1 (MQ=255)
                                                                               cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGGAGGACAAACTGGAAAAGCTGAAAAGGAGTCTg                                                                     >  2:588921/1‑150 (MQ=255)
                                                                               cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCt                                                                      >  1:783141/1‑149 (MQ=255)
                                                                                aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGcc                                                                   <  1:91305/151‑1 (MQ=255)
                                                                                     ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGGAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                                 >  2:634752/1‑148 (MQ=255)
                                                                                       cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                            <  2:551297/151‑1 (MQ=255)
                                                                                               tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCCAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAAAGCTGAAAAGGAGGTTGCCGGGAGGGGTggg                                                       >  1:577750/1‑147 (MQ=255)
                                                                                               tGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGGGGAGGGATCGCCGCCGTGGAGGACAAACTGGAAACGCTGAAAAGTAGGCTGCCGGGAGGCgggg                                                        >  2:309648/1‑147 (MQ=255)
                                                                                                gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGAt                                                    <  2:52588/150‑1 (MQ=255)
                                                                                                    cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCCGCAAAAACGCCCGAGGAGGGATCGCCGCCGTGAAGGGAAAACTGGGAAAGCTGAAAAATAGTCTGCCGGGAGGGGGGGGGAAAGtt                                                >  1:357532/1‑150 (MQ=255)
                                                                                                       ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGGAGGGCAAAATGGAAAAGCTGAAAAGGAGTCTGCCGGAAGGCGGGGGGATAGTTaca                                             >  2:803381/1‑150 (MQ=255)
                                                                                                        gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGCGGAGGACAAAACGGAAACGCTGAAAAGGAGTCTGCCCGAAGGCGGGGAGATAGTTacaa                                            >  2:404052/1‑150 (MQ=255)
                                                                                                          aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAACCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAAATGGAAACGCTGAAAAGAAGTCTGCCGGGAGGGGTGGGGATAGTTtaaaca                                          >  2:351765/1‑150 (MQ=255)
                                                                                                           ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGGGCTCCCGCAAAAAAACCCGAGGAGTGAACGCCGCCGCGAAGGGCAAACTTGGAACCCTGCAAAGTAGTCTCCCGGACGGGGTGGAGATaagtaaaa                                            >  1:943856/1‑140 (MQ=255)
                                                                                                             gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCTGCGATCCGGCAAAAAAGACCGAGAAGCGATCGCCGCCGGGAAGGACAAAATGGAAACGCTGAAAAGTAGGCTGCCAGGAGGCGGGGAGAAAGTTACAACATa                                        >  1:180627/1‑149 (MQ=255)
                                                                                                                gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                   <  1:826483/151‑1 (MQ=255)
                                                                                                                     aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcag                               <  2:490625/150‑1 (MQ=255)
                                                                                                                     aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                                <  2:321497/149‑1 (MQ=255)
                                                                                                                      aaCGACGAAGGCGAAGTGGCGGGCGGGGCGGGGAATCTCCGATACGGCAAAAACGCCCGCGAAGAGATATCCCGCATGGAGGGCAAAAATGGGAACACGAAAAATAGTGAGCCGGAAGGGGTGGGAAGAGTTTCAACATACGAttccaccc                             >  2:537847/1‑144 (MQ=255)
                                                                                                                          gcgaaggcgaagTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcc                             >  1:179068/1‑147 (MQ=255)
                                                                                                                                       ggcgggcgGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGcc            >  1:953412/1‑151 (MQ=255)
                                                                                                                                         cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt          <  1:448190/151‑1 (MQ=255)
                                                                                                                                          gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt          <  1:168658/150‑1 (MQ=255)
                                                                                                                                              ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc      >  2:807292/1‑150 (MQ=255)
                                                                                                                                               gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaa    <  2:767311/151‑1 (MQ=255)
                                                                                                                                                 gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc  >  2:570552/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                       
GCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599430‑599725

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