Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,590,2130TC23.3% 1516.9 / inf 554T254T (ACT→ACCPA14_RS21015site‑specific integrase
Reads supporting (aligned to +/- strand):  ref base T (215/209);  new base C (66/63);  total (282/272)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 561 total aligned reads displayed.

GCAGGAGGCCGATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGCCGACTACGACCACATGCGTGCCGTGCTG  >  NC_008463/4590063‑4590363
                                                                                                                                                      |                                                                                                                                                      
 cAGGAGGCCGATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTc                                                                                                                                                       <  1:7846789/151‑1 (MQ=33)
             agagCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTcc                                                                                                                                           <  1:2713500/151‑1 (MQ=33)
               agCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCgg                                                                                                                                         <  1:136297/151‑1 (MQ=34)
                 ccGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGca                                                                                                                                       >  2:2976906/1‑151 (MQ=35)
                   ggATTCGCCGATGACGACCAGTTGTTCACCGTCAACCGGTTCTCACCGCACTTCGAGAGCACGGTGTTGAACTCCGACCAGGTCGAAGCCATGTACCGGCAGTTGACCGAGAAGGGTGGGGGGCGGATGACCCCGCACCGTTTCCGGcaca                                                                                                                                     >  2:4405115/1‑151 (MQ=11)
                     aTTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACAcc                                                                                                                                   <  1:1210501/151‑1 (MQ=35)
                        cGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTg                                                                                                                                >  2:10774077/1‑151 (MQ=35)
                        cGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTg                                                                                                                                <  2:1339678/151‑1 (MQ=35)
                        cGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTg                                                                                                                                <  1:3960714/151‑1 (MQ=35)
                        cGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTg                                                                                                                                >  1:11212164/1‑151 (MQ=35)
                        cGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTg                                                                                                                                >  2:10673600/1‑151 (MQ=14)
                        cGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTg                                                                                                                                <  1:6961604/151‑1 (MQ=14)
                          ccGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGc                                                                                                                              >  1:3022968/1‑151 (MQ=14)
                          ccGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGGTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGc                                                                                                                              >  2:6539731/1‑151 (MQ=14)
                           cGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGcc                                                                                                                             >  1:1556331/1‑151 (MQ=35)
                           cGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGcc                                                                                                                             >  2:2791517/1‑151 (MQ=35)
                            gATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGACGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCCCACCCTGGCCa                                                                                                                            >  2:8830010/1‑151 (MQ=25)
                             aTGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCAc                                                                                                                           >  2:1958321/1‑151 (MQ=35)
                              tGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCAcc                                                                                                                          <  2:3398865/151‑1 (MQ=35)
                               gacgacCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCg                                                                                                                         <  1:2046858/151‑1 (MQ=35)
                               gacgacCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCg                                                                                                                         <  1:9219313/151‑1 (MQ=35)
                                 cgacCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGAc                                                                                                                       <  1:7609011/151‑1 (MQ=14)
                                  gacCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACt                                                                                                                      <  1:10916694/151‑1 (MQ=14)
                                   acCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACtt                                                                                                                     >  1:7526743/1‑151 (MQ=35)
                                   acCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACtt                                                                                                                     <  2:10967154/151‑1 (MQ=35)
                                   acCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACtt                                                                                                                     <  1:10212400/151‑1 (MQ=14)
                                       gttgttCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTtgatg                                                                                                                 >  1:10586134/1‑151 (MQ=35)
                                        ttgttCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTtgatga                                                                                                                >  2:7673209/1‑151 (MQ=14)
                                           ttCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAgg                                                                                                             >  2:8990605/1‑151 (MQ=14)
                                           ttCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAgg                                                                                                             <  1:5737865/151‑1 (MQ=14)
                                            tCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGc                                                                                                            <  1:2202223/151‑1 (MQ=35)
                                            tCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGc                                                                                                            <  2:6346590/151‑1 (MQ=35)
                                              aaCGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCAc                                                                                                          >  1:7974609/1‑151 (MQ=14)
                                               aCGTCAACGGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCAcc                                                                                                         >  2:167202/1‑151 (MQ=25)
                                               aCGTCAACCGGTTCTCACCGCACTACAAGTGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGACGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGGCCACCATGGCCACCGACTTGATGAAGGCAcc                                                                                                         >  2:10587394/1‑151 (MQ=255)
                                                cGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCAccc                                                                                                        >  1:7738357/1‑151 (MQ=14)
                                                 gTCAGCCGGTTCTCTCCGCACTACAAGAGCAAGGTGAGGAACTGCGACCAGGTCGAAGCCTTGTACCGGTAGTTGACCGTGTGGGGTGGGGTGCGGATGAGTCGGCAGCGGTTCGGGCAGACCCTGGCCACCGGCGGGGGGATGGCACCCg                                                                                                       >  2:7335603/1‑151 (MQ=2)
                                                 gTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCg                                                                                                       >  1:8462235/1‑151 (MQ=35)
                                                  tCAACCGGTTCTCACCGCACTCCAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTAGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGa                                                                                                      <  2:5385701/151‑1 (MQ=25)
                                                  tCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGa                                                                                                      <  2:9406454/151‑1 (MQ=35)
                                                   cAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAg                                                                                                     <  1:8999828/151‑1 (MQ=35)
                                                   cAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAg                                                                                                     <  1:4246763/151‑1 (MQ=14)
                                                    aaCCGGTTCTCACCGCACGACAAGAGCAAGGTGATGAACTCCGACCCGGTCGCTGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGc                                                                                                    <  2:6851049/151‑1 (MQ=255)
                                                     aCCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCg                                                                                                   <  2:8235550/151‑1 (MQ=14)
                                                         gTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAAc                                                                                               <  2:1872942/151‑1 (MQ=35)
                                                          ttCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACa                                                                                              >  2:4813280/1‑151 (MQ=14)
                                                           tctcACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACAt                                                                                             <  1:4530608/151‑1 (MQ=14)
                                                             tcACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATcc                                                                                           <  1:6723750/151‑1 (MQ=37)
                                                              cACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCa                                                                                          <  2:4162965/151‑1 (MQ=38)
                                                               aCCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCAc                                                                                         >  2:7544062/1‑151 (MQ=38)
                                                                  gCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTc                                                                                      >  1:3276346/1‑151 (MQ=38)
                                                                    aCTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCAc                                                                                    >  1:3892799/1‑151 (MQ=38)
                                                                    aCTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCAc                                                                                    >  1:1683794/1‑151 (MQ=38)
                                                                    aCTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCAc                                                                                    <  1:1889157/151‑1 (MQ=38)
                                                                    aCTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCAc                                                                                    <  1:3413356/151‑1 (MQ=38)
                                                                    aCTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACCCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATACACCTCAc                                                                                    >  2:3096530/1‑151 (MQ=14)
                                                                       aCAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGaa                                                                                 >  2:916347/1‑151 (MQ=38)
                                                                       aCAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGaa                                                                                 >  2:3850767/1‑151 (MQ=38)
                                                                        cAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAg                                                                                >  2:2503938/1‑151 (MQ=38)
                                                                        cAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAg                                                                                <  2:2227590/151‑1 (MQ=38)
                                                                            agCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGcc                                                                            <  1:6120215/151‑1 (MQ=38)
                                                                               aaGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGAACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGc                                                                         >  2:5565400/1‑151 (MQ=38)
                                                                                aGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCt                                                                        <  1:9814496/151‑1 (MQ=38)
                                                                                aGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCt                                                                        <  1:10052072/151‑1 (MQ=38)
                                                                                    gatgaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAAc                                                                    <  1:8139307/151‑1 (MQ=38)
                                                                                      tgaACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCa                                                                  >  2:3649508/1‑151 (MQ=38)
                                                                                        aaCTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACt                                                                <  1:5509672/151‑1 (MQ=38)
                                                                                         aCTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTc                                                               <  1:6319806/151‑1 (MQ=38)
                                                                                           tCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGa                                                             <  2:7450029/151‑1 (MQ=38)
                                                                                            ccGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGaa                                                            >  2:5884078/1‑151 (MQ=38)
                                                                                            ccGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGaa                                                            >  1:8859470/1‑151 (MQ=38)
                                                                                            ccGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGaa                                                            >  2:1515613/1‑151 (MQ=38)
                                                                                              gACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAata                                                          >  2:10173989/1‑151 (MQ=38)
                                                                                               aCCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAatat                                                         <  1:9918388/151‑1 (MQ=38)
                                                                                                  aGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCa                                                      >  2:11108919/1‑151 (MQ=38)
                                                                                                     tCGAAGCCATGTACCGTAAGTTGACCGAGAAGGTTGGGGTGCGGTTGAATCCGCACCGTTTCCGGCTCACCCTGGCCACCGACTTGATGAAGGCAGCCGAGCGGAACATCCACCTCACGATGTGCCTGCTCAACCACTCGAATATCCAGac                                                   <  2:3946428/151‑1 (MQ=11)
                                                                                                     tCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGac                                                   <  1:4726343/151‑1 (MQ=38)
                                                                                                       gtaGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGacca                                                 <  2:1568795/149‑1 (MQ=38)
                                                                                                       gAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGacca                                                 <  1:3722715/151‑1 (MQ=38)
                                                                                                          gCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGAccacca                                              >  2:10356028/1‑151 (MQ=39)
                                                                                                           ccATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCAt                                             <  1:5953643/151‑1 (MQ=40)
                                                                                                              tGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAg                                          <  2:5007064/151‑1 (MQ=40)
                                                                                                              tGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAg                                          <  1:730135/151‑1 (MQ=40)
                                                                                                                 aCCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCtc                                       >  1:2657099/1‑150 (MQ=40)
                                                                                                                    ggAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCTGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACAt                                    <  2:3246007/151‑1 (MQ=255)
                                                                                                                     gAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATc                                   <  1:8219552/151‑1 (MQ=40)
                                                                                                                     gAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATc                                   >  2:1718882/1‑151 (MQ=40)
                                                                                                                       aGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGa                                 >  2:9734234/1‑151 (MQ=40)
                                                                                                                          tGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGc                              >  1:5212312/1‑151 (MQ=40)
                                                                                                                           gACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGcc                             >  1:5345792/1‑151 (MQ=40)
                                                                                                                             ccGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGCCGa                           <  1:4109390/151‑1 (MQ=40)
                                                                                                                               gagaAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGAATTGATGAAGGCACCCGAGCGGTACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGCCGACt                         <  1:1956689/151‑1 (MQ=255)
                                                                                                                                agaAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGCCGACTa                        <  1:1086023/151‑1 (MQ=40)
                                                                                                                                 gaAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGCCGACTAc                       >  1:8577301/1‑151 (MQ=40)
                                                                                                                                   aGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGATCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGCCGACTACGa                     >  2:858757/1‑151 (MQ=255)
                                                                                                                                       tGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGCCGACTACGACcac                 >  1:2330620/1‑151 (MQ=255)
                                                                                                                                        ggggTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGCCGACTACGACcaca                <  1:5042229/151‑1 (MQ=255)
                                                                                                                                             gCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGCCGACTACGACCACATGcgt           <  2:8128957/151‑1 (MQ=255)
                                                                                                                                               ggATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGCCGACTACGACCACAAGcgtgc         >  2:11182502/1‑151 (MQ=255)
                                                                                                                                                    aCTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGCCGACTACGACCACATGcgtgccgtgc    <  1:9588166/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GCAGGAGGCCGATAGAGCCGGATTCGCCGATGACGACCAGTTGTTCAACGTCAACCGGTTCTCACCGCACTACAAGAGCAAGGTGATGAACTCCGACCAGGTCGAAGCCATGTACCGGAAGTTGACCGAGAAGGTTGGGGTGCGGATGACTCCGCACCGTTTCCGGCACACCCTGGCCACCGACTTGATGAAGGCACCCGAGCGGAACATCCACCTCACGAAGTGCCTGCTCAACCACTCGAATATCCAGACCACCATGAGCTACATCGAGGCCGACTACGACCACATGCGTGCCGTGCTG  >  NC_008463/4590063‑4590363

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