Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,855,2560AG31.7% 610.1 / inf 364N2230S (AAC→AGC) PA14_RS13370hemagglutinin repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (137/111);  new base G (58/57);  total (195/169)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.28e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 369 total aligned reads displayed.

ACCTTCGCGCCGACACGCTCGACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAACGGCGACGGCGGGCGGCTCGACAGCCA  >  NC_008463/2855106‑2855406
                                                                                                                                                      |                                                                                                                                                      
  cTTCGCGCCGACACGCTCGACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAAcc                                                                                                                                                      >  1:4683170/1‑151 (MQ=35)
     cgcgCCGACACGCTCGACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACcgcg                                                                                                                                                   >  2:4392176/1‑151 (MQ=33)
          cGACACGCTCGACAACAGCCAGTCCGGGCTGATCGCCGCGACTGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCggcggc                                                                                                                                              <  1:1026906/151‑1 (MQ=14)
           gACACGCTCGACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCa                                                                                                                                             <  2:5753466/151‑1 (MQ=31)
             cacGCTCGACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAg                                                                                                                                           >  2:7898294/1‑151 (MQ=31)
               cGCTCGACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGt                                                                                                                                         >  2:5253401/1‑151 (MQ=31)
                tcTCGACAACAGCCAGTCCGGGCTGATCGCCGCGAAAGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGACGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGTCGGCAAGGTc                                                                                                                                        >  2:9070550/2‑151 (MQ=255)
                gCTCGACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTc                                                                                                                                        <  1:1415712/151‑1 (MQ=14)
                gCTCGACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCAGggcggcgcgctg                                                                                                                                        <  1:3281683/151‑7 (MQ=9)
                   cGACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATc                                                                                                                                     <  2:9765081/151‑1 (MQ=14)
                   cGACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATc                                                                                                                                     <  1:6904465/151‑1 (MQ=31)
                    gACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTCTCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGAATTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCACCTGGACAACCGCGGCGGCAAGGTCATCg                                                                                                                                    >  2:7094032/1‑151 (MQ=9)
                    gACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCg                                                                                                                                    <  1:2666063/151‑1 (MQ=14)
                    gACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCg                                                                                                                                    >  1:4457363/1‑151 (MQ=14)
                     acaacaGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCgg                                                                                                                                   <  1:2018385/151‑1 (MQ=14)
                     acaacaGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGc                                                                                                                                  <  1:10632047/151‑1 (MQ=14)
                      caacaGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGc                                                                                                                                  >  1:9179105/1‑151 (MQ=31)
                      caacaGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGc                                                                                                                                  >  2:8525343/1‑151 (MQ=31)
                      caacaGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGc                                                                                                                                  <  1:9456829/151‑1 (MQ=31)
                       aacaGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCg                                                                                                                                 >  1:2565126/1‑151 (MQ=14)
                       aacaGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCg                                                                                                                                 <  1:2992133/151‑1 (MQ=31)
                       aacaGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCg                                                                                                                                 >  1:4890351/1‑151 (MQ=31)
                             cAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCg                                                                                                                           <  1:5834105/151‑1 (MQ=31)
                              aGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCgg                                                                                                                          <  2:3855145/151‑1 (MQ=14)
                                tCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGcc                                                                                                                        <  1:4636779/151‑1 (MQ=31)
                                tCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGcc                                                                                                                        <  2:2280768/151‑1 (MQ=31)
                                 ccGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCt                                                                                                                       <  2:4963533/151‑1 (MQ=31)
                                  cGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTg                                                                                                                      >  1:1369989/1‑151 (MQ=14)
                                   gggCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTgc                                                                                                                     <  1:1957820/151‑1 (MQ=31)
                                    ggCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTgcg                                                                                                                    <  2:1914324/151‑1 (MQ=14)
                                    ggCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTgcg                                                                                                                    <  2:10842726/151‑1 (MQ=14)
                                     gCTGATCGCCGCGAATGGCGGCATCGCTATCGTGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTgcgc                                                                                                                   >  1:5224681/1‑151 (MQ=14)
                                     gCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTgcgc                                                                                                                   >  2:7266773/1‑151 (MQ=31)
                                      cTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGcc                                                                                                                  >  1:168981/1‑151 (MQ=31)
                                      cTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGcc                                                                                                                  >  2:3664135/1‑151 (MQ=31)
                                       tGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCt                                                                                                                 <  1:7403642/151‑1 (MQ=31)
                                            gccgcGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCg                                                                                                            <  1:6317885/151‑1 (MQ=14)
                                              cgcgAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTg                                                                                                          >  2:9195399/1‑151 (MQ=31)
                                                  aaTGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAgc                                                                                                      >  1:3875075/1‑151 (MQ=31)
                                                     ggcggcATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGcc                                                                                                   >  2:8598565/1‑151 (MQ=31)
                                                      gcggcATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCACGGTCAGCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCct                                                                                                  >  2:10113731/1‑151 (MQ=14)
                                                        ggcATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCctgc                                                                                                >  2:2564598/1‑151 (MQ=31)
                                                         gcATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCctgct                                                                                               <  1:23629/151‑1 (MQ=11)
                                                                                        ccGCGCCGGCGAGAGTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGt                                                                <  2:5104688/151‑1 (MQ=11)
                                                                                           cgcCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCAACGTGCAGCGCCTGCTGAACCAGGCCGTGGGGGTGCTGGCCGGGCGTGAc                                                             >  2:4561176/1‑151 (MQ=9)
                                                                                             ccGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACgg                                                           >  2:1726097/1‑151 (MQ=31)
                                                                                             ccGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACgg                                                           >  1:2110663/1‑151 (MQ=31)
                                                                                              cGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGCCAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCCCTGGGGTGCTGGCCGGGCGTGACGGc                                                          >  1:1074560/1‑151 (MQ=9)
                                                                                              cGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGc                                                          <  1:9396096/151‑1 (MQ=31)
                                                                                               ggCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGcc                                                         >  1:7249267/1‑151 (MQ=31)
                                                                                                gCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCt                                                        <  1:2719278/151‑1 (MQ=31)
                                                                                                 cGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTg                                                       >  1:7755347/1‑151 (MQ=31)
                                                                                                 cGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTg                                                       >  2:5886709/1‑151 (MQ=31)
                                                                                                  gagaTTTCCAGCACCTCGATGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGa                                                      >  1:7095097/1‑151 (MQ=14)
                                                                                                   agaTTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAg                                                     >  1:816367/1‑151 (MQ=14)
                                                                                                   agaTTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAg                                                     >  1:6231801/1‑151 (MQ=14)
                                                                                                     aTTTCCAGCACCTCGCAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGcc                                                   <  1:764516/151‑1 (MQ=18)
                                                                                                      tttCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCt                                                  >  2:8080421/1‑151 (MQ=33)
                                                                                                        tCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTgg                                                <  1:9877916/151‑1 (MQ=33)
                                                                                                         ccAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGa                                               >  2:4822280/1‑151 (MQ=18)
                                                                                                          cAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGAc                                              <  2:10752734/151‑1 (MQ=33)
                                                                                                          cAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGAc                                              <  2:10300376/151‑1 (MQ=33)
                                                                                                          cAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGAc                                              >  1:3185327/1‑151 (MQ=33)
                                                                                                           aGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGAcg                                             <  1:3780224/151‑1 (MQ=18)
                                                                                                             cACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGAcggc                                           >  1:9686680/1‑151 (MQ=18)
                                                                                                              aCCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGAcggcg                                          <  2:6349074/151‑1 (MQ=18)
                                                                                                              aCCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGAcggcg                                          >  2:6654762/1‑151 (MQ=33)
                                                                                                               ccTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGAcggcgg                                         >  1:4772022/1‑151 (MQ=33)
                                                                                                                cTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACggcggc                                        <  2:8261902/151‑1 (MQ=33)
                                                                                                                cTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACggcggc                                        <  2:4111926/151‑1 (MQ=33)
                                                                                                                cTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACggcggc                                        <  1:8758845/151‑1 (MQ=33)
                                                                                                                   gAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGaa                                     <  2:3011768/151‑1 (MQ=33)
                                                                                                                     aGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACt                                   <  2:9463981/151‑1 (MQ=33)
                                                                                                                      ggTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTg                                  <  1:435592/151‑1 (MQ=33)
                                                                                                                       gTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGt                                 >  2:2752773/1‑151 (MQ=33)
                                                                                                                           ccGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCaa                             >  2:3455101/1‑151 (MQ=18)
                                                                                                                           ccGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCaa                             >  1:5919993/1‑151 (MQ=33)
                                                                                                                            cGTGAACGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAac                            >  2:3794640/1‑151 (MQ=18)
                                                                                                                             gTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacg                           >  1:1242716/1‑151 (MQ=33)
                                                                                                                               gAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggc                         <  1:7996150/151‑1 (MQ=33)
                                                                                                                                aaCGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggcg                        <  1:8624249/151‑1 (MQ=18)
                                                                                                                                aaCGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggcg                        <  1:1924907/151‑1 (MQ=18)
                                                                                                                                aaCGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggcg                        >  2:5705236/1‑151 (MQ=33)
                                                                                                                                 aCGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggcga                       >  1:5490982/1‑151 (MQ=33)
                                                                                                                                  cGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggcgac                      <  2:5388463/151‑1 (MQ=18)
                                                                                                                                  cGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggcgac                      <  1:5072876/151‑1 (MQ=18)
                                                                                                                                  cGCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggcgac                      <  2:2048937/151‑1 (MQ=18)
                                                                                                                                  cGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggcgac                      <  1:10098408/151‑1 (MQ=33)
                                                                                                                                   gCCCGCGAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggcgacg                     <  2:4734062/151‑1 (MQ=18)
                                                                                                                                    cccGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggcgacgg                    <  1:2484267/151‑1 (MQ=33)
                                                                                                                                    cccGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggcgacgg                    <  2:8763740/151‑1 (MQ=33)
                                                                                                                                     ccGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAacggcgacggc                   <  2:4469145/151‑1 (MQ=33)
                                                                                                                                       gcgAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAACGGCGAcggcgg                 <  1:655969/151‑1 (MQ=33)
                                                                                                                                       gcgAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAACGGCGAcggcgg                 <  1:6749345/151‑1 (MQ=33)
                                                                                                                                         gAGCAACTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAACGGCGACggcgggc               <  1:5760677/151‑1 (MQ=25)
                                                                                                                                         gAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAACGGCGACggcgggc               >  2:3515669/1‑151 (MQ=35)
                                                                                                                                             aaCTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAACGGCGACGGCGGGCGGCt           >  2:1367969/1‑151 (MQ=25)
                                                                                                                                              aCTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAACGGCGACGGCGGGCGGCTc          >  1:1502014/1‑151 (MQ=35)
                                                                                                                                               cTGGACAGCCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAACGGCGACGGCGGGCGGCTCg         >  2:6093160/1‑151 (MQ=25)
                                                                                                                                                 ggACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAACGGCGACGGCGGGCGGCTCGAc       >  1:7883751/1‑151 (MQ=35)
                                                                                                                                                      |                                                                                                                                                      
ACCTTCGCGCCGACACGCTCGACAACAGCCAGTCCGGGCTGATCGCCGCGAATGGCGGCATCGCTATCGAGGCGCGGCAGGTCGACAACCGCGCCGGCGAGATTTCCAGCACCTCGAAGGTCGCCGTGAACGCCCGCGAGCAACTGGACAACCGCGGCGGCAAGGTCATCGGCGACAGCGGCCTGCGCCTCACCGTGCAGCGCCTGCTGAACCAGGCCAAGGGGGTGCTGGCCGGGCGTGACGGCCTGAGCCTGGACGGCGGCGAACTGTTCAACGGCGACGGCGGGCGGCTCGACAGCCA  >  NC_008463/2855106‑2855406

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: