Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_008463114,7940AG43.2% 898.7 / inf 755K505R (AAG→AGG) vgrG1btype VI secretion system spike protein VgrG1b
Reads supporting (aligned to +/- strand):  ref base A (209/219);  new base G (171/154);  total (382/373)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 772 total aligned reads displayed.

CCGAGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCGACC  >  NC_008463/114644‑114943
                                                                                                                                                      |                                                                                                                                                     
  gAGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGagaaga                                                                                                                                                     <  1:7703625/151‑1 (MQ=33)
   aGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAgaagaa                                                                                                                                                    <  1:4380395/151‑1 (MQ=33)
     cAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAAcc                                                                                                                                                  >  1:526386/1‑151 (MQ=33)
          ggTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGAc                                                                                                                                             >  1:3079799/1‑151 (MQ=33)
            tGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACAt                                                                                                                                           >  1:10082878/1‑151 (MQ=33)
             gCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATc                                                                                                                                          <  1:1941171/151‑1 (MQ=34)
              cccTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCg                                                                                                                                         <  2:898209/151‑1 (MQ=35)
                 tACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAgg                                                                                                                                      <  1:9938279/151‑1 (MQ=35)
                   cGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACACGAAGGGCGCCGAGCAGTTGTTCATCCCTGCCGAGACGAACCAGGACATCGAGGTc                                                                                                                                    >  2:7144191/1‑151 (MQ=255)
                          ccGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAacg                                                                                                                             <  1:7279300/151‑1 (MQ=38)
                            ggCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAacgac                                                                                                                           >  2:5841927/1‑151 (MQ=38)
                            ggCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAacgac                                                                                                                           >  1:4645359/1‑151 (MQ=38)
                             gCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAacgacg                                                                                                                          <  1:10305492/151‑1 (MQ=38)
                                aaCGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACgaga                                                                                                                       >  2:3380612/1‑151 (MQ=38)
                                 aCGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGagag                                                                                                                      >  1:9647046/1‑151 (MQ=38)
                                     taCCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCAc                                                                                                                  >  1:1511493/2‑151 (MQ=14)
                                            agCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCg                                                                                                           <  1:4795379/151‑1 (MQ=25)
                                              cGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGc                                                                                                         >  1:3174404/1‑151 (MQ=40)
                                               ggCACCTAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGTAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGcc                                                                                                        <  2:1439558/151‑1 (MQ=255)
                                               ggCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGcc                                                                                                        <  1:2713825/151‑1 (MQ=40)
                                                  aCCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACg                                                                                                     >  1:4110546/1‑151 (MQ=40)
                                                    cAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGAc                                                                                                   <  1:8442951/151‑1 (MQ=40)
                                                       gagCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCgg                                                                                                >  2:1275942/1‑151 (MQ=40)
                                                          ccGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAATAAGGGCGCGGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGACGGCCACGACCGGAcc                                                                                             <  1:166070/151‑1 (MQ=255)
                                                          ccGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGAcc                                                                                             >  1:3724840/1‑151 (MQ=25)
                                                          ccGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGAcc                                                                                             >  2:3841201/1‑151 (MQ=40)
                                                           cGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACCAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGGCGGCCACGACCGGACCa                                                                                            >  1:4238249/1‑151 (MQ=255)
                                                            gTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCaa                                                                                           >  1:9793509/1‑151 (MQ=25)
                                                             ttCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAg                                                                                          >  2:10248172/1‑151 (MQ=40)
                                                              tCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGa                                                                                         >  2:459747/1‑151 (MQ=25)
                                                                gAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGGCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGAcc                                                                                       >  1:5562451/1‑151 (MQ=255)
                                                                 aGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGACCGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCa                                                                                      >  2:1854253/1‑151 (MQ=25)
                                                                 aGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCa                                                                                      <  2:5425695/151‑1 (MQ=25)
                                                                  gCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCAt                                                                                     <  2:5514363/151‑1 (MQ=25)
                                                                  gCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCAGGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGGAGAACCAGGACATCGAGGTCGAGACCGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCAt                                                                                     <  2:6084640/151‑1 (MQ=255)
                                                                   cAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATc                                                                                    >  2:8449345/1‑151 (MQ=37)
                                                                   cAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATc                                                                                    <  1:6283104/151‑1 (MQ=255)
                                                                    aaGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCg                                                                                   >  2:7340241/1‑151 (MQ=38)
                                                                      gggCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGATGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGAc                                                                                 <  1:745990/151‑1 (MQ=255)
                                                                         cggcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTac                                                                              >  1:12026/1‑151 (MQ=255)
                                                                          ggcACGCCGTCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCacg                                                                             <  2:245448/151‑1 (MQ=255)
                                                                             aCGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCacgacg                                                                          <  2:497521/151‑1 (MQ=25)
                                                                             aCGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTacgacg                                                                          >  2:6942798/1‑151 (MQ=255)
                                                                                 cGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGAc                                                                      >  1:2404265/1‑151 (MQ=25)
                                                                                 cGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGAc                                                                      >  1:9128126/1‑151 (MQ=255)
                                                                                  ggCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACg                                                                     >  1:2830142/1‑151 (MQ=255)
                                                                                  ggCCAACTTCAAAGAGATCCGCATTGAGGACACGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAAACAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACg                                                                     <  2:747965/151‑1 (MQ=11)
                                                                                    ccAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGt                                                                   >  2:4454351/1‑151 (MQ=255)
                                                                                    ccAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCAACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGt                                                                   >  1:8528770/1‑151 (MQ=255)
                                                                                     cAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTg                                                                  >  1:3368614/1‑151 (MQ=25)
                                                                                      aaCTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGc                                                                 >  2:2991774/1‑151 (MQ=255)
                                                                                         ttCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGACGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCAtg                                                              >  2:1668654/1‑151 (MQ=255)
                                                                                          tCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCAtgt                                                             <  2:2913067/151‑1 (MQ=255)
                                                                                           cAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCAtgtg                                                            <  2:8315176/151‑1 (MQ=255)
                                                                                             aCGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGaa                                                          >  1:2170568/1‑151 (MQ=25)
                                                                                              cGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAg                                                         <  1:8826142/151‑1 (MQ=25)
                                                                                              cGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAg                                                         <  2:10567021/151‑1 (MQ=25)
                                                                                               gagaTCCGCATGGAGGACAAGAAGGGCGCCGAGCTGTAGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGc                                                        <  1:7905546/151‑1 (MQ=255)
                                                                                                 gaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGACAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCAc                                                      <  1:5757185/151‑1 (MQ=255)
                                                                                                  aTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACg                                                     <  1:6078194/151‑1 (MQ=255)
                                                                                                   tCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGa                                                    >  2:2084567/1‑151 (MQ=255)
                                                                                                    ccGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGAt                                                   <  1:5440276/151‑1 (MQ=25)
                                                                                                    ccGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGAt                                                   >  2:9334207/1‑151 (MQ=25)
                                                                                                     cGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATc                                                  <  1:8211929/151‑1 (MQ=255)
                                                                                                     cGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATc                                                  >  1:863252/1‑151 (MQ=255)
                                                                                                     cGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATc                                                  <  1:9483654/151‑1 (MQ=255)
                                                                                                      gCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCg                                                 >  2:4300323/1‑151 (MQ=25)
                                                                                                      gCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCg                                                 >  2:9481729/1‑151 (MQ=255)
                                                                                                          ggaggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCAcc                                             <  2:3587153/151‑1 (MQ=25)
                                                                                                             ggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCgag                                          <  1:10795668/151‑1 (MQ=25)
                                                                                                              gaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCgaga                                         <  2:7920992/151‑1 (MQ=25)
                                                                                                               aCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGAc                                        >  1:9221185/1‑151 (MQ=255)
                                                                                                                cAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGAcc                                       >  1:10572632/1‑151 (MQ=25)
                                                                                                                 aagaagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCg                                      <  2:6291431/151‑1 (MQ=25)
                                                                                                                 aagaagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCg                                      >  2:2587731/1‑151 (MQ=25)
                                                                                                                   gaagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTg                                    <  1:7056318/151‑1 (MQ=25)
                                                                                                                   gaagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTg                                    >  1:8043752/1‑151 (MQ=25)
                                                                                                                   gaagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCTAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTg                                    <  2:8809566/151‑1 (MQ=255)
                                                                                                                   gaagGGCGCCGAGCAGTTGTACATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTg                                    <  1:3936909/151‑1 (MQ=255)
                                                                                                                    aagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTgg                                   >  1:2562064/1‑151 (MQ=25)
                                                                                                                     agGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGa                                  >  2:9723637/1‑151 (MQ=25)
                                                                                                                           ccGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGAcaacaa                            <  2:9101807/151‑1 (MQ=25)
                                                                                                                           ccGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCTACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGAcaacaa                            >  1:3594229/1‑151 (MQ=255)
                                                                                                                            cGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGAGGACaacaac                           >  2:10413553/1‑151 (MQ=255)
                                                                                                                              aGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACga                         >  1:6567547/1‑151 (MQ=25)
                                                                                                                                 aGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGAc                      >  2:5242586/1‑151 (MQ=25)
                                                                                                                                 aGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGAc                      >  1:8852510/1‑151 (MQ=255)
                                                                                                                                  gttgttCATCCATGCCGAGAGGAACCAGGACCTCGAGGTCGTGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCTCCGAGACCGTGGACATCAACGAGAcc                     <  2:7719996/151‑1 (MQ=11)
                                                                                                                                  gttgttCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGAcc                     <  2:6671873/151‑1 (MQ=25)
                                                                                                                                  gttgttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGAcc                     <  2:10648576/151‑1 (MQ=255)
                                                                                                                                   ttgttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCa                    >  2:4800175/1‑151 (MQ=255)
                                                                                                                                   ttgttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCa                    <  2:1124059/151‑1 (MQ=255)
                                                                                                                                    tgttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCAt                   <  2:4133733/151‑1 (MQ=255)
                                                                                                                                    tgttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCAt                   <  2:651553/151‑1 (MQ=255)
                                                                                                                                      ttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATaa                 <  1:135648/151‑1 (MQ=255)
                                                                                                                                      ttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATaa                 <  2:8668817/151‑1 (MQ=255)
                                                                                                                                       tCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAAc                >  2:2345291/1‑151 (MQ=255)
                                                                                                                                          tCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATCAccgt             <  2:3242495/151‑3 (MQ=12)
                                                                                                                                             aTGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATCAccgtgca          >  1:5067733/1‑146 (MQ=2)
                                                                                                                                                     aagaaCCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCGAcc  <  1:6826404/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCGAGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCGACC  >  NC_008463/114644‑114943

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