Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_008463114,7940AG39.1% 1225.6 / inf 810K505R (AAG→AGG) vgrG1btype VI secretion system spike protein VgrG1b
Reads supporting (aligned to +/- strand):  ref base A (234/258);  new base G (156/160);  total (392/418)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 833 total aligned reads displayed.

CCGAGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCGACCG  >  NC_008463/114644‑114944
                                                                                                                                                      |                                                                                                                                                      
ccGAGCAGTCGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGATCAAGGGCGGCACGCCGGACAACTTCAACGAGATCCGCTTGGATGACAAGAAGGGCGCCGAGCATTTGTTCATCCATGCCGtgaa                                                                                                                                                        <  2:1526792/151‑1 (MQ=11)
ccGAGCAGACGGTGCCCTCCGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGagaa                                                                                                                                                        <  1:7662756/151‑1 (MQ=17)
   aGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAgaagaa                                                                                                                                                     >  1:6134380/1‑151 (MQ=33)
     cAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAAcc                                                                                                                                                   <  2:2059522/151‑1 (MQ=33)
      aGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCa                                                                                                                                                  >  1:1719723/1‑151 (MQ=33)
        aCGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAgg                                                                                                                                                <  1:7321967/151‑1 (MQ=33)
        aCGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAgg                                                                                                                                                <  2:10007840/151‑1 (MQ=33)
        aCGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAgg                                                                                                                                                <  2:10160462/151‑1 (MQ=33)
         cGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGa                                                                                                                                               <  1:10006475/151‑1 (MQ=33)
          ggTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGAc                                                                                                                                              <  2:5753731/151‑1 (MQ=33)
           gTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACa                                                                                                                                             >  2:4227475/1‑151 (MQ=33)
               ccTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGCCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGATGTTCATCCAAGCCGAGAAGAACCAGGACATCGa                                                                                                                                         >  2:7056445/1‑151 (MQ=255)
               ccTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTGCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAACGGCGCCGAGCAGATGATCATCCATGCCGAGAAGAACCAGGACCTCGa                                                                                                                                         >  2:1447349/1‑151 (MQ=255)
                  aCGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGt                                                                                                                                      >  2:3864185/1‑151 (MQ=35)
                    gAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCACGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCg                                                                                                                                    >  1:4306782/1‑151 (MQ=38)
                    gAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCg                                                                                                                                    >  1:3381901/1‑151 (MQ=38)
                     aGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCga                                                                                                                                   <  1:10413067/151‑1 (MQ=38)
                            ggCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAacgac                                                                                                                            >  1:7927952/1‑151 (MQ=38)
                             gCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAacgacg                                                                                                                           <  2:9236990/151‑1 (MQ=38)
                                       cccAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTg                                                                                                                 >  2:1825115/1‑151 (MQ=38)
                                        ccAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTgg                                                                                                                >  1:2363930/1‑151 (MQ=39)
                                        ccAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTgg                                                                                                                >  2:10335728/1‑151 (MQ=39)
                                         cAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTggg                                                                                                               <  1:3212904/151‑1 (MQ=40)
                                          agcgCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGt                                                                                                              <  2:3309197/148‑1 (MQ=25)
                                           gagCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTc                                                                                                             <  1:11157081/151‑1 (MQ=25)
                                           gagCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGACGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTc                                                                                                             >  1:1025786/1‑151 (MQ=38)
                                               ggCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGcc                                                                                                         >  2:7785364/1‑151 (MQ=25)
                                               ggCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGcc                                                                                                         >  1:5973267/1‑151 (MQ=25)
                                                gCACCAAGAGCCGATCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGATAGCCACTGGGTCGGCCa                                                                                                        >  2:11108428/1‑151 (MQ=255)
                                                  aCCAAGAGCCGTTCGAGCACGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACg                                                                                                      <  1:2120789/151‑1 (MQ=255)
                                                  aCCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACg                                                                                                      >  1:6653740/1‑151 (MQ=25)
                                                     aaGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGAcc                                                                                                   <  1:1854215/151‑1 (MQ=25)
                                                     aaGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGAcc                                                                                                   >  1:7818126/1‑151 (MQ=40)
                                                          ccGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCCACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGAcc                                                                                              <  2:3012122/151‑1 (MQ=25)
                                                            gTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCaa                                                                                            >  2:1686649/1‑151 (MQ=25)
                                                               cGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGAc                                                                                         >  2:2300956/1‑151 (MQ=40)
                                                                  gaaaGTGCGGCACGCCGGCCAACTTCAACGTGATCCGCATGGAGGACAAGAAGGGCGCCGTGCAGTTGTTCATCCATGCCGAGAGTAACAAGGACATCGAGGTCGTGAACGACGAGATCCACTGGGTCGGACACGACCGGACCAAGACCAt                                                                                      <  1:5449172/149‑1 (MQ=11)
                                                                  gCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCAt                                                                                      >  2:116862/1‑151 (MQ=40)
                                                                   cAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATc                                                                                     <  1:7103150/151‑1 (MQ=255)
                                                                    aaGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCg                                                                                    >  1:9010829/1‑151 (MQ=38)
                                                                    aaGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCg                                                                                    <  2:5595388/151‑1 (MQ=38)
                                                                    aaGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCg                                                                                    <  2:2557597/151‑1 (MQ=255)
                                                                     aGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGa                                                                                   <  2:4383459/151‑1 (MQ=38)
                                                                       ggcggcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACTAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGAcc                                                                                 <  2:8364924/151‑2 (MQ=255)
                                                                        gcggcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGTAGGGCGCCGAGCTGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGAcca                                                                                <  2:5941057/151‑3 (MQ=255)
                                                                         cggcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTac                                                                               <  1:3686792/151‑1 (MQ=255)
                                                                               gCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACgag                                                                         <  1:8747763/151‑1 (MQ=25)
                                                                               gCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACgag                                                                         >  2:7590223/1‑151 (MQ=255)
                                                                                ccGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCAAGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACgaga                                                                        <  2:10283932/151‑1 (MQ=255)
                                                                                 cGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGAc                                                                       >  1:6373884/1‑151 (MQ=255)
                                                                                 cGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGAc                                                                       <  2:630389/151‑1 (MQ=255)
                                                                                 cGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGAc                                                                       <  2:5116586/151‑1 (MQ=255)
                                                                                     cAACTTCAACGAGATCCGCATGGAGGACAAGACGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTg                                                                   <  1:8475898/151‑1 (MQ=255)
                                                                                          tCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCAtgt                                                              <  2:644260/151‑1 (MQ=25)
                                                                                            aaCGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGa                                                            <  2:6175146/151‑1 (MQ=25)
                                                                                            aaCGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGa                                                            >  2:8772727/1‑151 (MQ=25)
                                                                                             aCGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGaa                                                           <  2:5257902/151‑1 (MQ=255)
                                                                                               gagaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGc                                                         <  2:3688417/151‑1 (MQ=25)
                                                                                               gagaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGc                                                         >  2:7194339/1‑151 (MQ=25)
                                                                                               gagaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGc                                                         <  2:2799744/151‑1 (MQ=255)
                                                                                                agaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGCGAGGAACCAGGACATCGAGGTCGAGATCGACGAGAGCCACTGGGTCGGCCACGACCTGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCa                                                        <  2:7410839/151‑1 (MQ=255)
                                                                                                   tCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGa                                                     >  1:2567407/1‑151 (MQ=255)
                                                                                                    ccGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTACATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGAt                                                    <  2:4107708/151‑1 (MQ=255)
                                                                                                     cGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATc                                                   <  2:4560837/151‑1 (MQ=25)
                                                                                                         tGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCAc                                               >  1:3882515/1‑151 (MQ=25)
                                                                                                         tGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCAc                                               >  2:6778849/1‑151 (MQ=255)
                                                                                                            aggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCga                                            <  1:9863107/151‑1 (MQ=255)
                                                                                                               aCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGAc                                         <  2:7901054/151‑1 (MQ=25)
                                                                                                               aCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGAc                                         <  1:5080939/151‑1 (MQ=25)
                                                                                                               aCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGAc                                         <  2:9785227/151‑1 (MQ=255)
                                                                                                                cAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGAcc                                        <  1:5247411/151‑1 (MQ=25)
                                                                                                                 aagaagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCg                                       <  1:9158042/151‑1 (MQ=25)
                                                                                                                  agaagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGTCCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGt                                      <  1:10681386/151‑1 (MQ=255)
                                                                                                                    aagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTgg                                    <  1:5694372/151‑1 (MQ=255)
                                                                                                                    aagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTgg                                    <  1:8233316/151‑1 (MQ=255)
                                                                                                                      gggCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGac                                  >  1:4852359/1‑151 (MQ=25)
                                                                                                                      gggCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGac                                  >  1:7244113/1‑151 (MQ=255)
                                                                                                                      gggCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGac                                  >  1:4708582/1‑151 (MQ=255)
                                                                                                                        gcgcCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGacaa                                <  1:1100587/151‑1 (MQ=255)
                                                                                                                         cgcCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGacaac                               <  1:9452987/151‑1 (MQ=25)
                                                                                                                         cgcCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGacaac                               <  1:3338160/151‑1 (MQ=255)
                                                                                                                          gcCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGacaaca                              >  1:1273229/1‑151 (MQ=25)
                                                                                                                          gcCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGacaaca                              <  2:10893551/151‑1 (MQ=255)
                                                                                                                           ccGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGAcaacaa                             <  1:2917430/151‑1 (MQ=25)
                                                                                                                           ccGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGAcaacaa                             <  2:1326134/151‑1 (MQ=255)
                                                                                                                           ccGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGAcaacaa                             >  2:8248789/1‑151 (MQ=255)
                                                                                                                             gAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACg                           <  2:9053591/151‑1 (MQ=25)
                                                                                                                             gAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACg                           >  1:2734651/1‑151 (MQ=25)
                                                                                                                               gCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACgag                         <  2:9472986/151‑1 (MQ=255)
                                                                                                                                 aGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGAc                       <  1:9642906/151‑1 (MQ=255)
                                                                                                                                  gttgttCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGAcc                      >  1:7475787/1‑151 (MQ=25)
                                                                                                                                  gttgttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGCGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGAcc                      <  1:10491536/151‑1 (MQ=255)
                                                                                                                                       tCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCTGACCAAGACCATCGACTACTACCAGACGGTGCATGTGAAGAACGAGCTCACCGAGACCGTGGACGACAACGAGACCATAAc                 >  1:9329039/1‑151 (MQ=255)
                                                                                                                                        cATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATCAcc                >  2:9404703/1‑151 (MQ=16)
                                                                                                                                        cATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATAGCACCGAGACCGTGGACAACAACGAGACCATCAcc                <  1:408846/151‑1 (MQ=255)
                                                                                                                                         aTCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCa               >  1:1555396/1‑151 (MQ=255)
                                                                                                                                          tCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATCAccgt              >  2:3167724/1‑149 (MQ=12)
                                                                                                                                          tCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCAt              <  2:10491259/151‑1 (MQ=255)
                                                                                                                                            cATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATCAccgtgc            >  2:11003420/1‑147 (MQ=9)
                                                                                                                                            cATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCg            >  2:5436301/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CCGAGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCGACCG  >  NC_008463/114644‑114944

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: ‑