Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_008463114,7940AG43.3% 805.1 / inf 875K505R (AAG→AGG) vgrG1btype VI secretion system spike protein VgrG1b
Reads supporting (aligned to +/- strand):  ref base A (244/249);  new base G (182/195);  total (428/447)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.30e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 888 total aligned reads displayed.

CGAGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCGACC  >  NC_008463/114645‑114943
                                                                                                                                                     |                                                                                                                                                     
cGAGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGagaa                                                                                                                                                       >  1:872494/1‑150 (MQ=32)
      gACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGACGAACCa                                                                                                                                                 >  2:6481806/1‑150 (MQ=17)
        cGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAgg                                                                                                                                               <  2:1760373/150‑1 (MQ=34)
        cGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAgg                                                                                                                                               <  2:7089778/150‑1 (MQ=34)
          gTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCTGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGAc                                                                                                                                             <  2:1154139/150‑1 (MQ=21)
            gCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACAt                                                                                                                                           <  1:8384326/150‑1 (MQ=34)
             cccTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATc                                                                                                                                          >  2:9436709/1‑150 (MQ=35)
             cccTACGAGCTGCCGGCGAACGCCACCCAGAGCGCCACCAACACCCGTTCGAGCAAGCGCGTCACGCCCGCCCACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTAGTTCATCCATGCCGAGAAGAACCCGGACATc                                                                                                                                          <  2:3824521/150‑1 (MQ=11)
                  cGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGt                                                                                                                                     <  2:4827114/150‑1 (MQ=37)
                  cGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGt                                                                                                                                     >  1:3146587/1‑150 (MQ=37)
                   gAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTc                                                                                                                                    >  1:712438/1‑150 (MQ=38)
                      cTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCgag                                                                                                                                 <  1:273530/150‑1 (MQ=39)
                        gCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGaa                                                                                                                               <  2:7902971/150‑1 (MQ=39)
                        gCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGaa                                                                                                                               >  1:4327729/1‑150 (MQ=39)
                           ggCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAacga                                                                                                                            >  1:10019018/1‑150 (MQ=39)
                           ggCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAacga                                                                                                                            <  1:7690654/150‑1 (MQ=39)
                                   ccACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGcc                                                                                                                    >  1:8240381/1‑150 (MQ=38)
                                   ccACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGcc                                                                                                                    <  2:5851173/150‑1 (MQ=38)
                                     aCCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACACGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCAc                                                                                                                  <  2:9688522/150‑1 (MQ=38)
                                     aCCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCAc                                                                                                                  <  2:10379164/150‑1 (MQ=17)
                                     aCCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCAc                                                                                                                  <  1:4162705/150‑1 (MQ=17)
                                     aCCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCAc                                                                                                                  <  1:6432801/150‑1 (MQ=38)
                                     aCCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCAc                                                                                                                  >  1:3530804/1‑150 (MQ=38)
                                      cccAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACTAGAGCCACt                                                                                                                 >  2:5075815/1‑150 (MQ=18)
                                      cccAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACt                                                                                                                 >  1:1693647/1‑150 (MQ=18)
                                       ccAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTg                                                                                                                <  1:5016631/150‑1 (MQ=21)
                                        cAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTgg                                                                                                               >  2:3464524/1‑150 (MQ=25)
                                        cAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACGTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTgg                                                                                                               <  1:5293270/150‑1 (MQ=255)
                                         agagCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCTGCCAACTTCAACGAGATCCGCATGGAGGACAAGTAGGGCGCCGAGCAGTTGTTCATCCAGGCCGAGAAGAACCCGGACATCGAGGTCGAGAACGACGAGAGCCACTggg                                                                                                              <  1:3715615/150‑1 (MQ=255)
                                         agagCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTggg                                                                                                              >  1:3257599/1‑150 (MQ=25)
                                         agagCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTggg                                                                                                              >  2:1989278/1‑150 (MQ=25)
                                          gagCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGt                                                                                                             >  1:10968504/1‑150 (MQ=255)
                                            gCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCg                                                                                                           >  2:3253320/1‑150 (MQ=255)
                                              ggCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGc                                                                                                         >  1:10095547/1‑150 (MQ=37)
                                              ggCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGc                                                                                                         <  1:301864/150‑1 (MQ=37)
                                              ggCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGc                                                                                                         >  1:2306688/1‑150 (MQ=255)
                                                cACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCa                                                                                                       >  2:5183120/1‑150 (MQ=37)
                                                cACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCa                                                                                                       <  2:4951001/150‑1 (MQ=37)
                                                 aCCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCAc                                                                                                      >  1:3920450/1‑150 (MQ=37)
                                                      gagCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCg                                                                                                 <  2:221110/150‑1 (MQ=37)
                                                      gagCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCg                                                                                                 <  1:6979705/150‑1 (MQ=37)
                                                         ccGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGAc                                                                                              >  1:5167608/1‑150 (MQ=255)
                                                          cGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGAcc                                                                                             >  1:7910279/1‑150 (MQ=37)
                                                              cGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGa                                                                                         >  1:4544463/1‑150 (MQ=37)
                                                               gAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGAc                                                                                        >  1:5842836/1‑150 (MQ=37)
                                                                aGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGAcc                                                                                       >  2:10574815/1‑150 (MQ=37)
                                                                 gCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCa                                                                                      <  2:6599904/150‑1 (MQ=37)
                                                                  cAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCAt                                                                                     <  2:2595684/150‑1 (MQ=37)
                                                                  cAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCAt                                                                                     <  1:10471545/150‑1 (MQ=37)
                                                                    aGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCg                                                                                   >  2:6907623/1‑150 (MQ=38)
                                                                     gggCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGa                                                                                  <  2:4172387/150‑1 (MQ=255)
                                                                       gcggcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACt                                                                                <  1:2123661/150‑1 (MQ=255)
                                                                        cggcACGCCGGCCACCTTCACCGAGATCCGCAGGGAGGACACGAAGGGCGCCGCGCAGTTGTTCCTCCATGCCGAGAGGAACCAGGACATCGAGGTCGCGCACGACGCGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGAcca                                                                               <  2:9089409/150‑3 (MQ=11)
                                                                           cACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTacga                                                                            >  2:4512051/1‑150 (MQ=255)
                                                                              gcagtccAACTTCAACGAGATCCGCATGGAGGAAAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGAAATCGAGGTCGAGATCGACGAGAGCCACTGGGTCGGCCACGACCGGACCATGACCATCGACCAcgacga                                                                         <  2:10036308/145‑1 (MQ=11)
                                                                              gCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCAcgacga                                                                         >  1:3464115/1‑150 (MQ=255)
                                                                                cGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACgaga                                                                       >  1:7257793/1‑150 (MQ=255)
                                                                                   ccAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACgg                                                                    <  1:7426543/150‑1 (MQ=255)
                                                                                    cAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGt                                                                   >  2:6040850/1‑150 (MQ=255)
                                                                                     aaCTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTg                                                                  >  2:10949541/1‑150 (MQ=255)
                                                                                     aaCTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTg                                                                  >  2:8188181/1‑150 (MQ=255)
                                                                                      aCTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGc                                                                 >  1:1388726/1‑150 (MQ=255)
                                                                                         tCACCGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCAtg                                                              <  2:9871722/150‑1 (MQ=255)
                                                                                         tCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCAtg                                                              >  2:3684406/1‑150 (MQ=255)
                                                                                          cAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCAtgt                                                             >  2:2925744/1‑150 (MQ=255)
                                                                                            aCGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGa                                                           >  1:9435168/1‑150 (MQ=255)
                                                                                              gagaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAg                                                         >  1:9008186/1‑150 (MQ=255)
                                                                                                gaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCa                                                       <  1:5565065/150‑1 (MQ=255)
                                                                                                 aTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCAc                                                      >  2:10211261/1‑150 (MQ=255)
                                                                                                   ccGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCCGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGa                                                    <  2:7645154/150‑1 (MQ=255)
                                                                                                        tGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCa                                               <  1:1129799/150‑1 (MQ=255)
                                                                                                          gaggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCAcc                                             <  1:2423997/150‑1 (MQ=255)
                                                                                                          gaggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCAcc                                             >  2:8491476/1‑150 (MQ=255)
                                                                                                           aggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGGCGCGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCg                                            <  2:1332353/150‑1 (MQ=255)
                                                                                                            ggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCga                                           >  2:8388282/1‑150 (MQ=255)
                                                                                                            ggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCga                                           >  2:5434253/1‑150 (MQ=255)
                                                                                                            ggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCga                                           >  1:10688495/1‑150 (MQ=255)
                                                                                                             gaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCgag                                          >  1:600980/1‑150 (MQ=255)
                                                                                                                  gaagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGt                                     >  2:1574615/1‑150 (MQ=255)
                                                                                                                   aagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTg                                    >  2:6285644/1‑150 (MQ=255)
                                                                                                                     gggCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGa                                  <  1:9714308/150‑1 (MQ=255)
                                                                                                                       gcgcCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGaca                                >  2:3020731/1‑150 (MQ=255)
                                                                                                                        cgcCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGACGGCCACGACCGGACCAACACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGacaa                               <  2:3321169/150‑1 (MQ=255)
                                                                                                                         gcCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGacaac                              >  1:6681581/1‑150 (MQ=255)
                                                                                                                           cGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGAcaacaa                            <  1:6216117/150‑1 (MQ=255)
                                                                                                                            gAGCAGTTGTTCATCCATGCCGTGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACaacaac                           >  2:3943628/1‑150 (MQ=255)
                                                                                                                            gAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACaacaac                           >  2:8810920/1‑150 (MQ=255)
                                                                                                                             aGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACg                          >  2:9119289/1‑150 (MQ=255)
                                                                                                                               cAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACgag                        <  2:6750972/150‑1 (MQ=255)
                                                                                                                               cAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACgag                        <  2:462099/150‑1 (MQ=255)
                                                                                                                                aGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACgaga                       >  1:4824794/1‑150 (MQ=255)
                                                                                                                                 gttgttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGAc                      >  1:7999677/1‑150 (MQ=255)
                                                                                                                                      tCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCAtca                 >  1:92160/1‑148 (MQ=255)
                                                                                                                                      tCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATaa                 >  1:2012094/1‑150 (MQ=255)
                                                                                                                                      tCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATaa                 <  2:274439/150‑1 (MQ=255)
                                                                                                                                      tCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATaa                 <  2:2437955/150‑1 (MQ=255)
                                                                                                                                      tCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATaa                 <  2:10710105/150‑1 (MQ=255)
                                                                                                                                      tCATCCATGCCGAGAACAACCAGGACATCGAGGTCGAGAACGACGAGAGCCTCTGGGTCGGCCACGACCGGACCAAGACCAGCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATaa                 <  2:7333085/150‑1 (MQ=255)
                                                                                                                                              gCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGc         >  1:988275/1‑150 (MQ=255)
                                                                                                                                                  agaagaACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCg     >  1:781803/1‑150 (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: ‑