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

CGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTC  >  NC_000913/599433‑599727
                                                                                                                                                   |                                                                                                                                                    
cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATc                                                                                                                                                  <  2:757190/151‑1 (MQ=255)
 gacgacTTTAATCACATCGTTTTAAAAGCCCGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCCGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGggggtggcac                                                                                                                                                 >  1:748193/1‑147 (MQ=255)
       tttAACCACATCGTTTTAAAAGCCAGTGAAAATCGCGCGCTCGTTTATCTGCGCGATGTTGCGCAGGCCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGGAGTGGCGGGG‑GGGGGGGGGGTcc                                                                                                                                                 >  1:773096/1‑145 (MQ=255)
        ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGAGccctc                                                                                                                                              >  1:1180211/1‑144 (MQ=255)
           aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCAGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGACCTGccaccccg                                                                                                                                       >  2:104378/1‑144 (MQ=255)
            tCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTCCTGCGGTCCg                                                                                                                                        >  1:821997/1‑149 (MQ=255)
            tCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCGATCCgg                                                                                                                                       >  2:1182405/1‑150 (MQ=255)
              acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGGTCCgg                                                                                                                                       >  2:616747/1‑148 (MQ=255)
               caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCGATcc                                                                                                                                         >  1:357804/1‑145 (MQ=255)
                 tCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCGATCCGGCaaaa                                                                                                                                  >  2:439565/1‑150 (MQ=255)
                     tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGc                                                                                                                              <  1:214617/150‑1 (MQ=255)
                       tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGGCCTGCGATCCGGCAAAAACGccc                                                                                                                            >  1:1395846/1‑150 (MQ=255)
                        aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCga                                                                                                                          <  2:1210633/151‑1 (MQ=255)
                           agacAGTGAAAATGGCGTGACCGTTTATCTGCGCGATGTTGCGAAGGTACAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGAGAAGTGGCGGGC‑GGGGGGGGGGTTCTTGGATGCGGCaaaaaagagcgag                                                                                                                         >  1:264607/4‑141 (MQ=255)
                           aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGaa                                                                                                                       <  2:1208271/151‑1 (MQ=255)
                            gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCGATCCGGCCAAAACGCCCgcg                                                                                                                         >  1:227225/1‑146 (MQ=255)
                             ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGt                                                                                                                     <  2:1292820/151‑1 (MQ=255)
                               aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTCCTGCGATCCGGGCAAAACGGCCggggaggg                                                                                                                    >  1:345117/1‑143 (MQ=255)
                                gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGa                                                                                                                   <  2:161330/150‑1 (MQ=255)
                                 tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTCCTGCGGTCCGGCAAAAAAGCCCGGGAAggg                                                                                                                    >  1:382694/1‑146 (MQ=255)
                                 tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGAt                                                                                                                  <  1:1359369/150‑1 (MQ=255)
                                   aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgc                                                                                                               <  1:901177/151‑1 (MQ=255)
                                    aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgcc                                                                                                              <  1:251428/151‑1 (MQ=255)
                                             gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAg                                                                                                     <  2:332065/151‑1 (MQ=255)
                                               ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGTCCT‑CGATCCGGCAAAAAACCCCGGGGAGGGGTCGCCGCCGGGGAgg                                                                                                    >  1:230444/1‑150 (MQ=255)
                                                   ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGGCCTGCGATCCGGCCAAAACGCCCGAGGAGGGGTGGCCGCCGCGCAgggaa                                                                                                 >  1:334056/1‑146 (MQ=255)
                                                           gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGGGAAGGGATGGCCGCCGCGGAGGAGAAAACGGaa                                                                                         >  1:801548/1‑149 (MQ=255)
                                                              aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                                    <  2:352557/151‑1 (MQ=255)
                                                               tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGGGAACCCGCGGTCCGGCAAAAACGCCCGAGGAGTGATTGCCGCCGTGAAGGACAAACCGGaaa                                                                                        >  2:67791/1‑146 (MQ=255)
                                                                   gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAAAGCTGgaa                                                                                >  2:1330255/1‑150 (MQ=255)
                                                                      aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagta                                                                            <  2:1307066/151‑1 (MQ=255)
                                                                        ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtagt                                                                          <  2:1011573/151‑1 (MQ=255)
                                                                          tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGGGGAGGGATCGCCGCCCTGAAGGAGAAAATGGAAACGCTGaaaa                                                                               >  2:101332/1‑144 (MQ=255)
                                                                             aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGcc                                                                     <  1:430955/151‑1 (MQ=255)
                                                                             aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGcc                                                                     >  2:1139186/1‑151 (MQ=255)
                                                                              gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCg                                                                    <  1:671119/151‑1 (MQ=255)
                                                                                 tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa                                                                  <  2:203052/150‑1 (MQ=255)
                                                                                 tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGGAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGcc                                                                     >  2:860246/1‑147 (MQ=255)
                                                                                  ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                                <  2:156116/151‑1 (MQ=255)
                                                                                  ggCCCGGAGATGCGCCGGGGCATTGCCGAACGAAACGGGGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                                <  1:643819/151‑1 (MQ=255)
                                                                                     ccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGATCCTGCGATCCCGCAAAAAAGCCCGAGAAGGGATTGCCGCCGCGGAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCg                                                             >  1:1087866/1‑151 (MQ=255)
                                                                                      cGGAGATGCGCCGGGGCATTGCCGAACTAAAGGGGGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                              <  2:433455/149‑1 (MQ=255)
                                                                                      cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGt                                                            <  1:98683/151‑1 (MQ=255)
                                                                                           aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga                                                       <  2:1390100/151‑1 (MQ=255)
                                                                                                 cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGGGAACGCCGCCGGGAAGGACAAAATGGAAAAGCTGAAAAGTAGGCTGCCGGAAGGGGGGGGGATAGtt                                                  >  1:347687/1‑150 (MQ=255)
                                                                                                 cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAAATGGAAACGCTGAAAAGGAGTCTGCCGGGAGGCGGGGGGATAGtt                                                  >  2:1061041/1‑150 (MQ=255)
                                                                                                  ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTac                                                <  1:1233981/151‑1 (MQ=255)
                                                                                                    ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGCTCCTGCGATCCGGGAAAAACGCCCGCGGAGGGAGCGGCGCCGGGAAGGAGAAAACGGGAAAGCCGAAAAATAGCCCGGCGGGAGGCGGGGGGATAGTTaca                                               >  1:1277/1‑150 (MQ=255)
                                                                                                      cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGAGCCTGCGATCCGGCAAAAAAGCCCGCGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGGGGGGATAGTTtcaa                                              >  1:214350/1‑149 (MQ=255)
                                                                                                      cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGGAGGCGTGGAGATAGTTacaa                                              >  2:1326949/1‑149 (MQ=255)
                                                                                                       aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGG‑GGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGGTCGCCGCCGTGGAGGACAAACTGGGAAAGCTGAAAAGAAGTCTGCCGGAAGGGGTGGAGATAGTTacaaaa                                            >  1:80145/1‑148 (MQ=255)
                                                                                                         tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGA‑GGGGGGGGGGTCCTGCGCTCCCGCAAAAACGCCCCTGGAGGGATGGCCGCCGCGGAGGCCAAACTGGAAACGCCGTAAAGTAGTCTGCCCGGAGGAGTGGGGATAGGTtcaac                                             >  2:1386638/1‑147 (MQ=255)
                                                                                                             gAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGGGCCTGCGATCCGGGAAAAAAGGCCGAGGAGGGATGGCCGCCGTGAAGGAAAAACTGGGAAAGCTGAAAAGGAGTCTGCCGGGAGGGGGGGGGATAGGTACAACATaag                                        >  1:101333/1‑146 (MQ=255)
                                                                                                              aaCTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCg                                    <  2:580787/151‑1 (MQ=255)
                                                                                                              aaCTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGa                                       >  2:633054/1‑148 (MQ=255)
                                                                                                               aCTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCg                                    <  1:212759/150‑1 (MQ=255)
                                                                                                               aCTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGGGGGGATCCTGCGAACCGGCAAAAACGGCCGAGAAGCGATCGCCGCCGTGCAGGACCAAATGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGGGGAGATATTTacacaaatccattgc                                   >  2:384409/1‑139 (MQ=255)
                                                                                                                  aaaCGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                                  <  2:104193/149‑1 (MQ=255)
                                                                                                                    aCGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcca                              <  2:296696/151‑1 (MQ=255)
                                                                                                                      ggcgaaggcgaaGTGGCGGGC‑GGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagccagc                            >  2:451182/1‑151 (MQ=255)
                                                                                                                               gaagTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc                   <  2:212263/151‑1 (MQ=255)
                                                                                                                                     gcgggc‑ggGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGcc              >  2:38726/1‑150 (MQ=255)
                                                                                                                                     gcgggc‑ggGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCa             <  1:820041/151‑1 (MQ=255)
                                                                                                                                            ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc        <  1:1282966/150‑1 (MQ=255)
                                                                                                                                            ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa       <  1:884796/151‑1 (MQ=255)
                                                                                                                                                  gtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTc  >  2:1309514/1‑150 (MQ=255)
                                                                                                                                                  gtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTc  <  1:885243/150‑1 (MQ=255)
                                                                                                                                                   tgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTc  >  2:1065952/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGC‑GGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTC  >  NC_000913/599433‑599727

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