Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_008463114,7940AG42.0% 997.1 / inf 778K505R (AAG→AGG) vgrG1btype VI secretion system spike protein VgrG1b
Reads supporting (aligned to +/- strand):  ref base A (224/227);  new base G (150/176);  total (374/404)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 799 total aligned reads displayed.

CCGAGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCGACCG  >  NC_008463/114644‑114944
                                                                                                                                                      |                                                                                                                                                      
ccGAGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGagaa                                                                                                                                                        <  1:3682474/151‑1 (MQ=32)
    gCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAAc                                                                                                                                                    >  2:7407426/1‑151 (MQ=33)
     cAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAAcc                                                                                                                                                   <  1:10613508/151‑1 (MQ=33)
     cAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAAcc                                                                                                                                                   <  2:3791880/151‑1 (MQ=33)
       gACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAg                                                                                                                                                 <  1:8791217/151‑1 (MQ=33)
         cGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGa                                                                                                                                               <  2:5699293/151‑1 (MQ=33)
             gCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATc                                                                                                                                           <  2:7641055/151‑1 (MQ=34)
                 tACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAgg                                                                                                                                       >  2:5977601/1‑151 (MQ=35)
                  aCGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGt                                                                                                                                      >  1:8524656/1‑151 (MQ=35)
                        tGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGAAGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCgagac                                                                                                                                >  2:4438501/1‑150 (MQ=255)
                                    ccACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGGCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCa                                                                                                                    <  1:6541028/151‑1 (MQ=38)
                                        ccAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTgg                                                                                                                <  1:5177006/151‑1 (MQ=21)
                                         cAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTggg                                                                                                               >  1:3325374/1‑151 (MQ=25)
                                            agCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGATCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCg                                                                                                            >  2:1368447/1‑151 (MQ=25)
                                             gCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCgg                                                                                                           <  1:1690255/151‑1 (MQ=25)
                                               ggCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGcc                                                                                                         >  1:5659666/1‑151 (MQ=40)
                                                 cACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCAc                                                                                                       <  1:9475070/151‑1 (MQ=25)
                                                 cACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCAc                                                                                                       <  2:8695891/151‑1 (MQ=40)
                                                  aCCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACg                                                                                                      <  2:4545243/151‑1 (MQ=40)
                                                   ccAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGa                                                                                                     >  2:783877/1‑151 (MQ=25)
                                                       gagCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCgg                                                                                                 >  2:2278276/1‑151 (MQ=25)
                                                       gagCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCgg                                                                                                 >  1:6712193/1‑151 (MQ=40)
                                                       gagCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACAGCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCgg                                                                                                 >  2:3486211/1‑151 (MQ=255)
                                                        agCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGa                                                                                                >  1:6341854/1‑151 (MQ=40)
                                                        agCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCAAGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGa                                                                                                >  2:1758325/1‑151 (MQ=25)
                                                          ccGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGAcc                                                                                              >  2:9339807/1‑151 (MQ=25)
                                                          ccGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGAcc                                                                                              <  1:6163244/151‑1 (MQ=25)
                                                          ccGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGAcc                                                                                              >  2:6506831/1‑151 (MQ=40)
                                                           cGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCa                                                                                             <  2:9544377/151‑1 (MQ=25)
                                                           cGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCa                                                                                             >  1:4865298/1‑151 (MQ=40)
                                                           cGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCa                                                                                             <  2:3562488/151‑1 (MQ=40)
                                                            gTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCaa                                                                                            >  2:984742/1‑151 (MQ=40)
                                                              tCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGa                                                                                          <  1:5869780/151‑1 (MQ=40)
                                                                gAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGAcc                                                                                        <  1:8016837/151‑1 (MQ=25)
                                                                  gCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCAt                                                                                      <  2:8315613/151‑1 (MQ=40)
                                                                   cAAGGGCGGCACGCCTGCCATCTTCAACGCGATCCGCATGTCGGACAAGAAGGGCGCGGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATc                                                                                     <  2:10603022/151‑1 (MQ=11)
                                                                   cAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATc                                                                                     <  1:4080621/151‑1 (MQ=37)
                                                                   cAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATc                                                                                     <  1:1307481/151‑1 (MQ=37)
                                                                     aGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGa                                                                                   <  2:111252/151‑1 (MQ=38)
                                                                        gcggcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACACGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGAcca                                                                                <  2:200132/151‑3 (MQ=255)
                                                                        gcggcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGAcca                                                                                <  1:8190554/151‑3 (MQ=37)
                                                                           gcACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTacga                                                                             >  2:3313600/1‑151 (MQ=255)
                                                                             aCGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCacgacg                                                                           >  2:1024817/1‑151 (MQ=25)
                                                                             aCGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCacgacg                                                                           <  2:5188303/151‑1 (MQ=25)
                                                                             aCGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTacgacg                                                                           >  1:1144535/1‑151 (MQ=255)
                                                                                ccGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACgaga                                                                        >  1:7666634/1‑151 (MQ=255)
                                                                                 cGGCCAACTTCAACGAGCTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGAc                                                                       <  2:10504521/151‑1 (MQ=255)
                                                                                 cGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAATGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGAc                                                                       <  2:608633/151‑1 (MQ=255)
                                                                                 cGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGTACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGAc                                                                       <  2:8049144/151‑1 (MQ=255)
                                                                                 cGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGAc                                                                       >  2:8104333/1‑151 (MQ=25)
                                                                                 cGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGAc                                                                       >  2:8724038/1‑151 (MQ=25)
                                                                                  ggCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACg                                                                      <  1:3459982/151‑1 (MQ=25)
                                                                                   gCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACgg                                                                     <  1:7518173/151‑1 (MQ=25)
                                                                                    ccAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCCGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGCGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGt                                                                    <  2:9178059/151‑1 (MQ=255)
                                                                                     cAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTg                                                                   >  2:3927221/1‑151 (MQ=25)
                                                                                     cAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTg                                                                   >  1:10281330/1‑151 (MQ=25)
                                                                                         ttCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCAtg                                                               <  1:3936276/151‑1 (MQ=255)
                                                                                            aaCGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGa                                                            <  1:2674474/151‑1 (MQ=25)
                                                                                            aaCGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGa                                                            >  2:6244758/1‑151 (MQ=255)
                                                                                             aCGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGaa                                                           <  1:6893398/151‑1 (MQ=255)
                                                                                             aCGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGaa                                                           <  2:6838661/151‑1 (MQ=255)
                                                                                              cGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGATGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAg                                                          >  2:2482756/1‑151 (MQ=255)
                                                                                                cgaTCCGCAGGGCGGCCACGACGGGCGCCGAGCAGGGGTTCATCCATGCCGAGCAGAACCAGGACATCGAGGACGAGAACGCCGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCAGCGACTACGACGAGACGGTGCATGTGAAGCa                                                        <  2:10747989/150‑1 (MQ=14)
                                                                                                agaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCa                                                        <  2:490318/151‑1 (MQ=255)
                                                                                                 gaTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCAc                                                       >  2:7228115/1‑151 (MQ=25)
                                                                                                  aTCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACg                                                      >  2:5801764/1‑151 (MQ=255)
                                                                                                    ccGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGGTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGAt                                                    <  2:4266764/151‑1 (MQ=255)
                                                                                                     cGCATGGAGGACAAGAAGGGCGCCGAGCAGTGGTTCAACCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATc                                                   <  1:10017562/151‑1 (MQ=255)
                                                                                                      gCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGACCGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCCCGATCg                                                  >  2:5379230/1‑151 (MQ=255)
                                                                                                      gCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCg                                                  >  2:10084348/1‑151 (MQ=25)
                                                                                                          ggaggaCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCAcc                                              >  1:8954506/1‑151 (MQ=255)
                                                                                                              gaCAAGTAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCgaga                                          <  2:8897642/151‑1 (MQ=255)
                                                                                                               aCAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGAc                                         >  1:9261217/1‑151 (MQ=255)
                                                                                                                  agaagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGt                                      >  1:5453758/1‑151 (MQ=255)
                                                                                                                   gaagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTg                                     <  1:2032997/151‑1 (MQ=25)
                                                                                                                   gaagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTg                                     <  1:4236823/151‑1 (MQ=255)
                                                                                                                    aagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTgg                                    >  1:2886335/1‑151 (MQ=25)
                                                                                                                    aagGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTgg                                    >  1:7300767/1‑151 (MQ=255)
                                                                                                                        gcgcCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGacaa                                <  2:3101341/151‑1 (MQ=255)
                                                                                                                          gcCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGacaaca                              >  2:8507014/1‑151 (MQ=255)
                                                                                                                           ccGAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGAcaacaa                             <  2:700444/151‑1 (MQ=25)
                                                                                                                            cGAGCAGTTGTTCATCCATGCCGAGTGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACaacaac                            <  1:10431891/151‑1 (MQ=255)
                                                                                                                            cGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACaacaac                            <  1:2512058/151‑1 (MQ=255)
                                                                                                                             gAGCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACg                           <  2:5539545/151‑1 (MQ=25)
                                                                                                                             gAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACg                           >  1:8362091/1‑151 (MQ=255)
                                                                                                                              aGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACga                          >  1:3096000/1‑151 (MQ=255)
                                                                                                                               gCAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCTGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACgag                         >  2:8999525/1‑151 (MQ=255)
                                                                                                                                cAGTTGTTCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACgaga                        >  2:2855317/1‑151 (MQ=25)
                                                                                                                                   ttgttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCa                     >  1:9684604/1‑151 (MQ=255)
                                                                                                                                    tgttCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCAt                    >  2:5241555/1‑151 (MQ=25)
                                                                                                                                     gttCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATa                   >  1:1514636/1‑151 (MQ=255)
                                                                                                                                      ttCATCCATGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCAtca                  <  1:7823101/151‑3 (MQ=19)
                                                                                                                                        cATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAAcc                <  1:1831386/151‑1 (MQ=255)
                                                                                                                                            cATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCg            >  2:4689538/1‑151 (MQ=255)
                                                                                                                                             aggCCGAGAAGCACCAGGACACCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCCAGACCATCGACTACGACGAGACGGTGCATTTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCgg           <  2:6699769/149‑1 (MQ=255)
                                                                                                                                              tGCCGAGAGGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACCACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATCAccgtgcat          >  1:4569029/1‑145 (MQ=0)
                                                                                                                                                ccGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGt        >  1:9133362/1‑151 (MQ=255)
                                                                                                                                                ccGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGt        <  2:1637429/151‑1 (MQ=255)
                                                                                                                                                 cGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTc       <  1:2307252/151‑1 (MQ=255)
                                                                                                                                                 cGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTc       <  1:10765697/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CCGAGCAGACGGTGCCCTACGAGCTGCCGGCGAACGCCACCCAGAGCGGCACCAAGAGCCGTTCGAGCAAGGGCGGCACGCCGGCCAACTTCAACGAGATCCGCATGGAGGACAAGAAGGGCGCCGAGCAGTTGTTCATCCATGCCGAGAAGAACCAGGACATCGAGGTCGAGAACGACGAGAGCCACTGGGTCGGCCACGACCGGACCAAGACCATCGACTACGACGAGACGGTGCATGTGAAGCACGATCGCACCGAGACCGTGGACAACAACGAGACCATAACCATCGGCGTCGACCG  >  NC_008463/114644‑114944

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: