Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,457,811 A→T 17.5% T6S (ACT→TCT)  paaX → transcriptional repressor of phenylacetic acid degradation paa operon, phenylacetyl‑CoA inducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,457,8110AT17.5% 790.6 / 152.5 301T6S (ACT→TCT) paaXtranscriptional repressor of phenylacetic acid degradation paa operon, phenylacetyl‑CoA inducer
Reads supporting (aligned to +/- strand):  ref base A (145/103);  new base T (24/29);  total (169/132)
Fisher's exact test for biased strand distribution p-value = 9.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 304 total aligned reads displayed.

GGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAGGGCTGGGATTTGGTGAGCGT  >  CP009273/1457662‑1457960
                                                                                                                                                     |                                                                                                                                                     
 gtgtTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATAc                                                                                                                                                      >  1:3034869/1‑150 (MQ=255)
 gtgtTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATAc                                                                                                                                                      >  1:3750222/1‑150 (MQ=255)
  tgtTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACt                                                                                                                                                     <  2:3353615/150‑1 (MQ=255)
    tttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCttt                                                                                                                                                   <  1:4008722/150‑1 (MQ=255)
     ttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACtttt                                                                                                                                                  >  1:998080/1‑150 (MQ=255)
         tCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATca                                                                                                                                              <  2:2079039/150‑2 (MQ=255)
         tCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCAAAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATcc                                                                                                                                              >  2:2417310/1‑150 (MQ=255)
         tCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCACTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGCAGCGAATGAGTAAACTTGCTACTTTTATcc                                                                                                                                              >  2:3113794/1‑150 (MQ=255)
          cTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCa                                                                                                                                             <  1:49643/150‑1 (MQ=255)
          cTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCa                                                                                                                                             >  1:4241906/1‑150 (MQ=255)
           tGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCaa                                                                                                                                            >  1:3328721/1‑150 (MQ=255)
            gcgcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAAc                                                                                                                                           >  2:3450899/1‑150 (MQ=255)
               cAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATg                                                                                                                                        >  1:2768578/1‑150 (MQ=255)
               cAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATg                                                                                                                                        <  2:1418320/150‑1 (MQ=255)
                aaTATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGc                                                                                                                                       >  1:2187673/1‑150 (MQ=255)
                aaTATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGc                                                                                                                                       >  1:33687/1‑150 (MQ=255)
                  tatTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTg                                                                                                                                     >  2:2999817/1‑150 (MQ=255)
                   atTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGt                                                                                                                                    >  1:3554088/1‑150 (MQ=255)
                   atTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGt                                                                                                                                    <  2:2562625/150‑1 (MQ=255)
                   atTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGt                                                                                                                                    >  2:645863/1‑150 (MQ=255)
                     tgttgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTaa                                                                                                                                  <  1:3719329/150‑1 (MQ=255)
                     tggtgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTaa                                                                                                                                  <  1:1707525/147‑1 (MQ=255)
                       ttgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAAc                                                                                                                                <  2:442811/150‑1 (MQ=255)
                           tGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTAAACGCTg                                                                                                                            <  1:1499712/150‑1 (MQ=255)
                           tGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTg                                                                                                                            >  2:3557686/1‑150 (MQ=255)
                            gCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGt                                                                                                                           >  2:2345937/1‑150 (MQ=255)
                             ccTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTAAACGCTGtt                                                                                                                          <  2:301907/150‑1 (MQ=255)
                               tGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTcc                                                                                                                        >  2:2688262/1‑150 (MQ=255)
                                   gTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTc                                                                                                                    >  1:3035026/1‑150 (MQ=255)
                                   gTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTc                                                                                                                    >  1:2824335/1‑150 (MQ=255)
                                     aTCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAg                                                                                                                  >  1:3253322/1‑150 (MQ=255)
                                       cGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTg                                                                                                                <  1:2680044/150‑1 (MQ=255)
                                        gTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTgg                                                                                                               >  2:2932163/1‑150 (MQ=255)
                                         tCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGc                                                                                                              >  2:3633445/1‑150 (MQ=255)
                                           tGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGcac                                                                                                            >  2:3190202/1‑150 (MQ=255)
                                              cccTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATc                                                                                                         >  2:865146/1‑150 (MQ=255)
                                                   gtggGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGa                                                                                                    <  1:145993/150‑1 (MQ=255)
                                                    tggGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAt                                                                                                   <  2:519093/150‑1 (MQ=255)
                                                     ggggTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAtc                                                                                                  <  1:236944/150‑1 (MQ=255)
                                                     ggggTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAtc                                                                                                  <  2:3995958/150‑1 (MQ=255)
                                                      gggTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAtct                                                                                                 <  2:4345051/150‑1 (MQ=255)
                                                       ggTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAtctc                                                                                                >  2:4256016/1‑150 (MQ=255)
                                                          aaaGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATctcctc                                                                                             >  2:2124245/1‑150 (MQ=255)
                                                          aaaGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCTGTCAGTGGCACATCTTTGATctcctc                                                                                             >  1:2041464/1‑150 (MQ=255)
                                                           aaGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCctct                                                                                            >  1:2336117/1‑150 (MQ=255)
                                                           aaGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCctct                                                                                            <  2:1547725/150‑1 (MQ=255)
                                                            aGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCctctc                                                                                           >  1:4340375/1‑150 (MQ=255)
                                                                ccAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTa                                                                                       <  1:750487/150‑1 (MQ=255)
                                                                 cAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTAt                                                                                      >  2:306687/1‑150 (MQ=255)
                                                                 cAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTAt                                                                                      >  2:1607059/1‑150 (MQ=255)
                                                                  aGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATg                                                                                     <  1:4217418/150‑1 (MQ=255)
                                                                   gggCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATgg                                                                                    >  2:1944925/1‑150 (MQ=255)
                                                                   gggCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATgg                                                                                    >  2:3915620/1‑150 (MQ=255)
                                                                    ggCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGt                                                                                   <  2:1202085/150‑1 (MQ=255)
                                                                       cAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGAt                                                                                >  1:2941416/1‑150 (MQ=255)
                                                                        aGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGAtt                                                                               <  2:1050076/150‑1 (MQ=255)
                                                                           aGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGc                                                                            <  2:1729557/150‑1 (MQ=255)
                                                                           aGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGc                                                                            >  1:3542876/1‑150 (MQ=255)
                                                                            gTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCt                                                                           <  1:3913681/150‑1 (MQ=255)
                                                                               gATACGACCTGTGCTATGATTCATAAATCACAACAATAACCACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCACCATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTc                                                                        <  2:3738819/150‑1 (MQ=255)
                                                                                aTACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTcc                                                                       <  2:318425/150‑1 (MQ=255)
                                                                                  aCGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCa                                                                     <  2:1587288/150‑1 (MQ=255)
                                                                                    gACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATc                                                                   <  2:4625842/150‑1 (MQ=255)
                                                                                     aCCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCg                                                                  >  1:662626/1‑150 (MQ=255)
                                                                                     aCCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCg                                                                  >  2:4058867/1‑150 (MQ=255)
                                                                                      ccTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCgt                                                                 <  1:2891959/150‑1 (MQ=255)
                                                                                      ccTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCgt                                                                 <  1:1150312/150‑1 (MQ=255)
                                                                                       cTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCgtg                                                                <  2:2187673/150‑1 (MQ=255)
                                                                                       cTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCgtg                                                                <  2:33781/150‑1 (MQ=255)
                                                                                            cTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCgtggtggt                                                           >  1:449354/1‑150 (MQ=255)
                                                                                              aagATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTACACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATATTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGa                                                         <  1:1628125/148‑1 (MQ=255)
                                                                                                gATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGaaa                                                       >  2:1848365/1‑150 (MQ=255)
                                                                                                 aTTCATAACTCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAAt                                                      <  1:568277/150‑1 (MQ=255)
                                                                                                    cATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTg                                                   <  1:1640726/150‑1 (MQ=255)
                                                                                                     aTAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTgg                                                  >  1:502662/1‑150 (MQ=255)
                                                                                                        aaTCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTg                                               >  2:2837219/1‑150 (MQ=255)
                                                                                                         aTCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTgg                                              <  2:2374555/150‑1 (MQ=255)
                                                                                                          tCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTggg                                             >  2:46592/1‑150 (MQ=255)
                                                                                                           cacaACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGt                                            >  1:1005068/1‑150 (MQ=255)
                                                                                                           cacaACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCAGTTCCGGTCAGTGTCACATCTTTGATCTCCTCTCTGTATGTTGATTCGCTTTCCCATCGTTGTGGTGAAATCTGGATGGGt                                            >  2:208211/1‑150 (MQ=255)
                                                                                                                 aaTAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTg                                      >  1:2264564/1‑150 (MQ=255)
                                                                                                                 aaTAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTg                                      >  1:1918351/1‑150 (MQ=255)
                                                                                                                   tAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGgc                                    >  2:1429897/1‑150 (MQ=255)
                                                                                                                        acaGACTGAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCtt                               <  1:1200247/150‑1 (MQ=255)
                                                                                                                         caGACTGAAGCGAATGATTAAACTTGATTCTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCttt                              >  2:1522511/1‑150 (MQ=255)
                                                                                                                           gACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGc                            <  2:2755434/150‑1 (MQ=255)
                                                                                                                              tGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTgg                         >  1:606941/1‑150 (MQ=255)
                                                                                                                               gAAGCGAATGAGTAAACTTGATTCTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCACATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGa                        <  2:2512287/150‑1 (MQ=255)
                                                                                                                                caGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGaa                       <  1:1532403/149‑1 (MQ=255)
                                                                                                                                  gCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAgg                     <  1:2573345/150‑1 (MQ=255)
                                                                                                                                   cGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAggg                    >  2:1124771/1‑150 (MQ=255)
                                                                                                                                    gAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAGGGc                   >  2:3492889/1‑150 (MQ=255)
                                                                                                                                    gAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAGGGc                   >  2:1784161/1‑150 (MQ=255)
                                                                                                                                    gAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAGGGc                   >  2:2496896/1‑150 (MQ=255)
                                                                                                                                             aaCTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAGGGCTGGGATTTg          <  2:4184667/150‑1 (MQ=255)
                                                                                                                                               cTTGATTCTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAGGGCTGGGATTTGGt        >  1:938204/1‑150 (MQ=255)
                                                                                                                                               cTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAGGGCTGGGATTTGGt        <  1:3450899/150‑1 (MQ=255)
                                                                                                                                               cTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAGGGCTGGGATTTGGt        >  2:780876/1‑150 (MQ=255)
                                                                                                                                                ttGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAGGGCTGGGATTTGGTg       >  2:3259984/1‑150 (MQ=255)
                                                                                                                                                 tGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAGGGCTGGGATTTGGTGa      <  1:3622042/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTCTGTATGGTGATTCGCTTTCCCATCGTGGTGGTGAAATCTGGTTGGGTAGTCTGGCTGCTTTGCTGGAAGGGCTGGGATTTGGTGAGCGT  >  CP009273/1457662‑1457960

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: