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

CGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTC  >  NC_000913/599433‑599727
                                                                                                                                                   |                                                                                                                                                    
cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATc                                                                                                                                                  <  2:757190/151‑1 (MQ=255)
 gacgacTTTAATCACATCGTTTTAAAAGCCCGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCCGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGggggtggcac                                                                                                                                                 >  1:748193/1‑147 (MQ=255)
       tttAACCACATCGTTTTAAAAGCCAGTGAAAATCGCGCGCTCGTTTATCTGCGCGATGTTGCGCAGGCCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGGAGTGGCGGGGGGGGGGGGGGTcc                                                                                                                                                 >  1:773096/1‑145 (MQ=255)
        ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGAGccctc                                                                                                                                              >  1:1180211/1‑144 (MQ=255)
           aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCAGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGACCTGccaccccg                                                                                                                                       >  2:104378/1‑144 (MQ=255)
            tCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCg                                                                                                                                        >  1:821997/1‑149 (MQ=255)
            tCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCgg                                                                                                                                       >  2:1182405/1‑150 (MQ=255)
              acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCgg                                                                                                                                       >  2:616747/1‑148 (MQ=255)
               caTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATcc                                                                                                                                         >  1:357804/1‑145 (MQ=255)
                 tCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCaaaa                                                                                                                                  >  2:439565/1‑150 (MQ=255)
                     tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGc                                                                                                                              <  1:214617/150‑1 (MQ=255)
                       tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGATCCGGCAAAAACGccc                                                                                                                            >  1:1395846/1‑150 (MQ=255)
                        aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCga                                                                                                                          <  2:1210633/151‑1 (MQ=255)
                           agacAGTGAAAATGGCGTGACCGTTTATCTGCGCGATGTTGCGAAGGTACAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGAGAAGTGGCGGGCGGGGGGGGGGTTCTTGGATGCGGCaaaaaagagcgag                                                                                                                         >  1:264607/4‑141 (MQ=255)
                           aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGaa                                                                                                                       <  2:1208271/151‑1 (MQ=255)
                            gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCCAAAACGCCCgcg                                                                                                                         >  1:227225/1‑146 (MQ=255)
                             ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGt                                                                                                                     <  2:1292820/151‑1 (MQ=255)
                               aGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGGGCAAAACGGCCggggaggg                                                                                                                    >  1:345117/1‑143 (MQ=255)
                                gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGa                                                                                                                   <  2:161330/150‑1 (MQ=255)
                                 tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCAAAAAAGCCCGGGAAggg                                                                                                                    >  1:382694/1‑146 (MQ=255)
                                 tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGAt                                                                                                                  <  1:1359369/150‑1 (MQ=255)
                                   aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgc                                                                                                               <  1:901177/151‑1 (MQ=255)
                                    aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgcc                                                                                                              <  1:251428/151‑1 (MQ=255)
                                             gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAg                                                                                                     <  2:332065/151‑1 (MQ=255)
                                               ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGTCCTCGATCCGGCAAAAAACCCCGGGGAGGGGTCGCCGCCGGGGAgg                                                                                                    >  1:230444/1‑150 (MQ=255)
                                                   ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGATCCGGCCAAAACGCCCGAGGAGGGGTGGCCGCCGCGCAgggaa                                                                                                 >  1:334056/1‑146 (MQ=255)
                                                           gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGGGAAGGGATGGCCGCCGCGGAGGAGAAAACGGaa                                                                                         >  1:801548/1‑149 (MQ=255)
                                                              aTGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCt                                                                                    <  2:352557/151‑1 (MQ=255)
                                                               tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGAACCCGCGGTCCGGCAAAAACGCCCGAGGAGTGATTGCCGCCGTGAAGGACAAACCGGaaa                                                                                        >  2:67791/1‑146 (MQ=255)
                                                                   gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAAAGCTGgaa                                                                                >  2:1330255/1‑150 (MQ=255)
                                                                      aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagta                                                                            <  2:1307066/151‑1 (MQ=255)
                                                                        ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtagt                                                                          <  2:1011573/151‑1 (MQ=255)
                                                                          tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGGGGAGGGATCGCCGCCCTGAAGGAGAAAATGGAAACGCTGaaaa                                                                               >  2:101332/1‑144 (MQ=255)
                                                                             aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGcc                                                                     <  1:430955/151‑1 (MQ=255)
                                                                             aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGcc                                                                     >  2:1139186/1‑151 (MQ=255)
                                                                              gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCg                                                                    <  1:671119/151‑1 (MQ=255)
                                                                                 tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa                                                                  <  2:203052/150‑1 (MQ=255)
                                                                                 tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGGAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGcc                                                                     >  2:860246/1‑147 (MQ=255)
                                                                                  ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                                <  2:156116/151‑1 (MQ=255)
                                                                                  ggCCCGGAGATGCGCCGGGGCATTGCCGAACGAAACGGGGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                                <  1:643819/151‑1 (MQ=255)
                                                                                     ccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGATCCCGCAAAAAAGCCCGAGAAGGGATTGCCGCCGCGGAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCg                                                             >  1:1087866/1‑151 (MQ=255)
                                                                                      cGGAGATGCGCCGGGGCATTGCCGAACTAAAGGGGGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                              <  2:433455/149‑1 (MQ=255)
                                                                                      cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGt                                                            <  1:98683/151‑1 (MQ=255)
                                                                                           aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga                                                       <  2:1390100/151‑1 (MQ=255)
                                                                                                 cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCGAGAAGGGAACGCCGCCGGGAAGGACAAAATGGAAAAGCTGAAAAGTAGGCTGCCGGAAGGGGGGGGGATAGtt                                                  >  1:347687/1‑150 (MQ=255)
                                                                                                 cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGACAAAATGGAAACGCTGAAAAGGAGTCTGCCGGGAGGCGGGGGGATAGtt                                                  >  2:1061041/1‑150 (MQ=255)
                                                                                                  ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTac                                                <  1:1233981/151‑1 (MQ=255)
                                                                                                    ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGCTCCTGCGATCCGGGAAAAACGCCCGCGGAGGGAGCGGCGCCGGGAAGGAGAAAACGGGAAAGCCGAAAAATAGCCCGGCGGGAGGCGGGGGGATAGTTaca                                               >  1:1277/1‑150 (MQ=255)
                                                                                                      cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGAGCCTGCGATCCGGCAAAAAAGCCCGCGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGGGGGGATAGTTtcaa                                              >  1:214350/1‑149 (MQ=255)
                                                                                                      cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGGAGTCTGCCGGGAGGCGTGGAGATAGTTacaa                                              >  2:1326949/1‑149 (MQ=255)
                                                                                                       aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGGTCGCCGCCGTGGAGGACAAACTGGGAAAGCTGAAAAGAAGTCTGCCGGAAGGGGTGGAGATAGTTacaaaa                                            >  1:80145/1‑148 (MQ=255)
                                                                                                         tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGAGGGGGGGGGGTCCTGCGCTCCCGCAAAAACGCCCCTGGAGGGATGGCCGCCGCGGAGGCCAAACTGGAAACGCCGTAAAGTAGTCTGCCCGGAGGAGTGGGGATAGGTtcaac                                             >  2:1386638/1‑147 (MQ=255)
                                                                                                             gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGATCCGGGAAAAAAGGCCGAGGAGGGATGGCCGCCGTGAAGGAAAAACTGGGAAAGCTGAAAAGGAGTCTGCCGGGAGGGGGGGGGATAGGTACAACATaag                                        >  1:101333/1‑146 (MQ=255)
                                                                                                              aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCg                                    <  2:580787/151‑1 (MQ=255)
                                                                                                              aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGa                                       >  2:633054/1‑148 (MQ=255)
                                                                                                               aCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCg                                    <  1:212759/150‑1 (MQ=255)
                                                                                                               aCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGAACCGGCAAAAACGGCCGAGAAGCGATCGCCGCCGTGCAGGACCAAATGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGGGGAGATATTTacacaaatccattgc                                   >  2:384409/1‑139 (MQ=255)
                                                                                                                  aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                                  <  2:104193/149‑1 (MQ=255)
                                                                                                                    aCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcca                              <  2:296696/151‑1 (MQ=255)
                                                                                                                      ggcgaaggcgaaGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagccagc                            >  2:451182/1‑151 (MQ=255)
                                                                                                                               gaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc                   <  2:212263/151‑1 (MQ=255)
                                                                                                                                     gcgggcggGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGcc              >  2:38726/1‑150 (MQ=255)
                                                                                                                                     gcgggcggGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCa             <  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)
                                                                                                                                                   |                                                                                                                                                    
CGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACCTC  >  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: