Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_008463114,7940AG40.9% 919.3 / inf 817K505R (AAG→AGG) vgrG1btype VI secretion system spike protein VgrG1b
Reads supporting (aligned to +/- strand):  ref base A (222/257);  new base G (137/195);  total (362/455)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 847 total aligned reads displayed.

CGAGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCGACC  >  NC_008463/114645‑114943
                                                                                                                                                     |                                                                                                                                                     
  aGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGagaaga                                                                                                                                                     >  1:8565480/1‑150 (MQ=34)
     aGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAAcc                                                                                                                                                  <  2:5151244/150‑1 (MQ=34)
      gACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGCAGGGCGCAGAGCAGTTGTGCATCCATGCCGAGAAGAAACa                                                                                                                                                 >  2:6757329/1‑150 (MQ=255)
       aCGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAg                                                                                                                                                >  1:6364958/1‑150 (MQ=34)
       aCGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAg                                                                                                                                                <  2:7243446/150‑1 (MQ=34)
        cGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAgg                                                                                                                                               <  2:3880812/150‑1 (MQ=34)
          gTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGAc                                                                                                                                             >  2:1707824/1‑150 (MQ=34)
                tACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAg                                                                                                                                       >  1:3088754/1‑150 (MQ=37)
                  cGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCTCGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGt                                                                                                                                     >  1:526007/1‑150 (MQ=37)
                  cGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGt                                                                                                                                     <  1:8862582/150‑1 (MQ=37)
                   gAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTc                                                                                                                                    <  2:8659207/150‑1 (MQ=38)
                       tGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCgaga                                                                                                                                <  2:9739485/150‑1 (MQ=39)
                       tGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCgaga                                                                                                                                <  2:5926910/150‑1 (MQ=39)
                          cGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAacg                                                                                                                             >  1:3727141/1‑150 (MQ=39)
                           ggCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAacga                                                                                                                            >  2:7049157/1‑150 (MQ=39)
                             cGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAacgacg                                                                                                                          <  1:6233666/150‑1 (MQ=39)
                               aaCGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACgag                                                                                                                        >  1:4684050/1‑150 (MQ=39)
                                   ctaCCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGcc                                                                                                                    <  1:8430812/148‑1 (MQ=11)
                                     aCCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCAc                                                                                                                  >  1:7972700/1‑150 (MQ=17)
                                       ccAGCGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACTAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTg                                                                                                                <  2:9353253/150‑1 (MQ=255)
                                       ccAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTg                                                                                                                >  1:4833831/1‑150 (MQ=39)
                                        cAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTgg                                                                                                               >  1:1093326/1‑150 (MQ=25)
                                        cAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTgg                                                                                                               <  2:8178864/150‑1 (MQ=40)
                                          gagCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGt                                                                                                             >  2:10423607/1‑150 (MQ=255)
                                           agCGGCACCAAGAGCCGTTCGAGCCAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTc                                                                                                            <  2:6224569/150‑1 (MQ=255)
                                           agCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTc                                                                                                            <  1:7217701/150‑1 (MQ=37)
                                            gCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCg                                                                                                           >  1:2908082/1‑150 (MQ=37)
                                            gCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCg                                                                                                           <  2:3445879/150‑1 (MQ=255)
                                            gCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCg                                                                                                           <  2:3259925/150‑1 (MQ=255)
                                            gCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCg                                                                                                           >  2:2731648/1‑150 (MQ=255)
                                                  ccAAGAGCCGCTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACg                                                                                                     >  1:214811/1‑150 (MQ=255)
                                                    cagagCCGTTCGAGCATGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGTGGACAAGTAGGGCGCCGAGCAGTTGTTCTTCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGAc                                                                                                   <  2:7912766/149‑1 (MQ=255)
                                                          cGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCTAGCAGTTGTTCATCCATGCCGAGAAGACCCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGAcc                                                                                             <  2:4140121/150‑1 (MQ=255)
                                                            ttCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCaa                                                                                           >  2:1980355/1‑150 (MQ=255)
                                                               gAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGAc                                                                                        <  1:6045348/150‑1 (MQ=37)
                                                                   aaGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATc                                                                                    >  1:10341515/1‑150 (MQ=255)
                                                                    aGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCg                                                                                   <  2:1134841/150‑1 (MQ=255)
                                                                    aGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCg                                                                                   >  2:10003591/1‑150 (MQ=255)
                                                                       gcggcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGAcc                                                                                >  2:10962611/1‑149 (MQ=255)
                                                                       gcggcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGAcc                                                                                <  2:8465890/150‑2 (MQ=255)
                                                                        cggcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTa                                                                               <  1:8647659/150‑1 (MQ=255)
                                                                         ggcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCac                                                                              <  1:9081918/150‑1 (MQ=255)
                                                                         ggcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTac                                                                              >  1:1573695/1‑150 (MQ=255)
                                                                          gcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTacg                                                                             >  1:5413571/1‑150 (MQ=255)
                                                                          gcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTacg                                                                             >  1:271395/1‑150 (MQ=255)
                                                                           cACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTacga                                                                            <  1:4485245/150‑1 (MQ=255)
                                                                            aCGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCacgac                                                                           <  2:3013498/150‑1 (MQ=255)
                                                                              gCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTAcgacga                                                                         <  2:4770525/150‑1 (MQ=255)
                                                                               ccGGCCAACTTCAACGAGATCCGCATGGAGGACCAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACgag                                                                        >  1:7624992/1‑150 (MQ=255)
                                                                                cGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACgaga                                                                       <  2:10699057/150‑1 (MQ=255)
                                                                                cGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACgaga                                                                       <  2:5389843/150‑1 (MQ=255)
                                                                                 ggCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGAc                                                                      >  1:9302535/1‑150 (MQ=255)
                                                                                  gCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACg                                                                     >  2:6592425/1‑150 (MQ=255)
                                                                                  gCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACg                                                                     <  2:9653026/150‑1 (MQ=255)
                                                                                   ccAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACgg                                                                    <  2:10095690/150‑1 (MQ=255)
                                                                                    cAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGt                                                                   <  2:4392414/150‑1 (MQ=255)
                                                                                     aaCTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGAAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTg                                                                  <  1:4045282/150‑1 (MQ=255)
                                                                                      aCTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGc                                                                 <  2:186302/150‑1 (MQ=255)
                                                                                      aCTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGc                                                                 >  1:7267860/1‑150 (MQ=255)
                                                                                       cTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCa                                                                >  2:735935/1‑150 (MQ=255)
                                                                                        ttCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCAt                                                               >  2:8997314/1‑150 (MQ=255)
                                                                                           aaCGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCAtgtg                                                            <  2:9183175/150‑1 (MQ=255)
                                                                                            aCGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGa                                                           <  2:3696852/150‑1 (MQ=255)
                                                                                             cGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGaa                                                          >  1:5022167/1‑150 (MQ=255)
                                                                                              gagaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAg                                                         <  2:9838769/150‑1 (MQ=255)
                                                                                               agaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGc                                                        >  2:2432791/1‑150 (MQ=255)
                                                                                               agaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGc                                                        <  1:5699027/150‑1 (MQ=255)
                                                                                                gaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCa                                                       <  2:1352478/150‑1 (MQ=255)
                                                                                                  tCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACg                                                     <  2:1564236/150‑1 (MQ=255)
                                                                                                   ccGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCCTCCATGCCGAGAGGAACCAGGACCTCGAGGTCGAGACCGACGAGAGCCACTGGGTCGGCCACGCCCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGa                                                    <  2:7695087/150‑1 (MQ=11)
                                                                                                    cGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGAt                                                   <  1:9793258/150‑1 (MQ=255)
                                                                                                    cGCAAGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGACAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAATCACGaa                                                   >  2:4529747/1‑149 (MQ=255)
                                                                                                      cATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCg                                                 >  1:8149290/1‑150 (MQ=255)
                                                                                                          gaggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCAcc                                             >  1:7009830/1‑150 (MQ=255)
                                                                                                            ggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCga                                           >  2:9522122/1‑150 (MQ=255)
                                                                                                            ggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCga                                           <  1:4126717/150‑1 (MQ=255)
                                                                                                              aCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCgaga                                         <  1:10663373/150‑1 (MQ=255)
                                                                                                                    agGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTgg                                   >  2:8301460/1‑150 (MQ=255)
                                                                                                                    agGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTgg                                   >  1:11169025/1‑150 (MQ=255)
                                                                                                                      ggCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGac                                 >  2:7553333/1‑150 (MQ=255)
                                                                                                                       gcgcCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGaca                                >  1:1698183/1‑150 (MQ=255)
                                                                                                                            gAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACaacaac                           <  1:3412205/150‑1 (MQ=255)
                                                                                                                            gAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACaacaac                           >  1:1399261/1‑150 (MQ=255)
                                                                                                                              gCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACga                         >  1:3278389/1‑150 (MQ=255)
                                                                                                                              gCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACga                         >  1:10205954/1‑150 (MQ=255)
                                                                                                                                aGTTGTTCATCCAAGCCGAGAAGAACCAGGACATCGACGTCGAGCACGACGAGCGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTCCGACGAGACGGTGCATGTGTATCACGATCGCACCGAGACCGTGGACAACAACgaga                       <  2:2240081/150‑1 (MQ=14)
                                                                                                                                 gttgttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGAc                      <  1:10503607/150‑1 (MQ=255)
                                                                                                                                   tgttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCa                    >  1:8322961/1‑150 (MQ=255)
                                                                                                                                    gttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCAt                   <  1:6525791/150‑1 (MQ=255)
                                                                                                                                    gttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCAt                   <  1:10008058/150‑1 (MQ=255)
                                                                                                                                        agccATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAAcc               <  2:2092151/148‑1 (MQ=255)
                                                                                                                                         tCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATCAccg              >  1:7403994/1‑149 (MQ=255)
                                                                                                                                            aTGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCg           <  1:6850881/150‑1 (MQ=255)
                                                                                                                                              gCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCCTAACCATCGGc         >  2:10281498/1‑150 (MQ=255)
                                                                                                                                              gCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGc         <  1:2485215/150‑1 (MQ=255)
                                                                                                                                               ccGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCg        >  1:1834036/1‑150 (MQ=255)
                                                                                                                                                 gagaAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTc      <  1:2348808/150‑1 (MQ=255)
                                                                                                                                                 gagaAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTc      >  1:8745749/1‑150 (MQ=255)
                                                                                                                                                  agaagaACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCg     <  2:4943762/150‑1 (MQ=255)
                                                                                                                                                   gaagaaCCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCGa    <  2:5683167/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGAGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCGACC  >  NC_008463/114645‑114943

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