Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,568,919 G→A 5.0% intergenic (+14/‑313) pcaK → / → pcaF‑I 4‑hydroxybenzoate transporter/beta‑ketoadipyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,568,9190GA 5.0% 642.8 / 13.6 197intergenic (+14/‑313)pcaK/pcaF‑I4‑hydroxybenzoate transporter/beta‑ketoadipyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base G (71/116);  new base A (7/3);  total (78/119)
Fisher's exact test for biased strand distribution p-value = 5.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

GCCTGGAGTGGTGCAACGCTGCTGGGCCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTA  >  NC_002947/1568771‑1569068
                                                                                                                                                    |                                                                                                                                                     
   tGGAGTGGTGCAACGCTGCTGGGCCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGc                                                                                                                                                   <  1:849061/150‑1 (MQ=255)
     gAGTGGTGCAACGCTGCTGGGCCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCgg                                                                                                                                                 <  2:1005381/150‑1 (MQ=255)
     gAGTGGTGCAACGCTGCTGGGCCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCgg                                                                                                                                                 <  2:3295083/150‑1 (MQ=255)
          gtgCATCGCTGCTGGGCCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCTCTGGCGACCGTGTGAGTGGTCGTCTAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGAGAGGGCCGCATCGCGGccctc                                                                                                                                            <  2:4221389/150‑3 (MQ=255)
          gtgCAACGCTGCTGGGCCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCATCCTGAGGGCGAGGGCCGCAACGCGGccctc                                                                                                                                            <  2:3367447/150‑3 (MQ=255)
          gtgCAACGCTGCTGGGCCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCAt                                                                                                                                            >  2:2580547/1‑150 (MQ=255)
          gtgCAACGCTGCTGGGCCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCAt                                                                                                                                            >  1:3011802/1‑150 (MQ=255)
           tgCAACGCTGCTGGGCCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCAtt                                                                                                                                           >  1:49541/1‑150 (MQ=255)
            gCAACGCTGCTGGGCCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTc                                                                                                                                          >  2:3140522/1‑150 (MQ=255)
             cAACGCTGCTGGGCCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTca                                                                                                                                         >  2:2685275/1‑150 (MQ=255)
                      tgggcctgggcTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAAc                                                                                                                                <  1:4293562/150‑1 (MQ=255)
                        ggcctgggcTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAAcgc                                                                                                                              <  1:1472503/150‑1 (MQ=255)
                         gcctgggcTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAAcgcc                                                                                                                             <  1:2134240/150‑1 (MQ=255)
                           ctgggcTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAAcgccgc                                                                                                                           <  2:2454254/150‑1 (MQ=255)
                           ctgggcTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAAcgccgc                                                                                                                           <  2:1342292/150‑1 (MQ=255)
                           ctgggcTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCACAAcgccgc                                                                                                                           <  1:123517/150‑1 (MQ=255)
                             gggcTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTc                                                                                                                         <  1:4335539/150‑1 (MQ=255)
                              ggcTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTcc                                                                                                                        <  2:2178695/150‑1 (MQ=255)
                               gcTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTccc                                                                                                                       >  2:2442575/1‑150 (MQ=255)
                                 tGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGc                                                                                                                     >  2:2296069/1‑150 (MQ=255)
                                  ggAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCACGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCa                                                                                                                    <  1:3512157/150‑1 (MQ=255)
                                   gAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAg                                                                                                                   <  2:3854677/150‑1 (MQ=255)
                                    aGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAgg                                                                                                                  <  2:1096746/150‑1 (MQ=255)
                                    aGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAgg                                                                                                                  <  1:3507243/150‑1 (MQ=255)
                                    aGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAgg                                                                                                                  <  2:1472575/150‑1 (MQ=255)
                                     gCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAggg                                                                                                                 >  1:2180279/1‑150 (MQ=255)
                                     gCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAggg                                                                                                                 <  1:4244493/150‑1 (MQ=255)
                                     gCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAggg                                                                                                                 <  1:180402/150‑1 (MQ=255)
                                         cGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGggcgg                                                                                                             <  2:2318687/150‑1 (MQ=255)
                                         cGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGggcgg                                                                                                             <  1:1314480/150‑1 (MQ=255)
                                               ggTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgt                                                                                                       <  2:2649712/150‑1 (MQ=255)
                                                gTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgtg                                                                                                      >  2:2974191/1‑150 (MQ=255)
                                                 tGCCGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGTCGTTGtgc                                                                                                     <  1:2765305/150‑1 (MQ=255)
                                                   cTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCACAACGCCGCCCCTGCGGGAGCGGCGTTAcgctg                                                                                                   <  2:393792/150‑1 (MQ=38)
                                                    tGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCGCAGCACAACGCCGCCCCTGCGGGAGCGGCGTTAcgctgt                                                                                                  <  1:1779993/150‑2 (MQ=38)
                                                    tGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGtgctgc                                                                                                  >  2:319620/1‑150 (MQ=255)
                                                          cGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGtgctgtgaatgg                                                                                            <  2:4254680/150‑8 (MQ=255)
                                                           gTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAgg                                                                                           <  1:1939735/150‑1 (MQ=255)
                                                             tGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGgc                                                                                         <  1:994873/150‑1 (MQ=255)
                                                                tGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGca                                                                                      >  2:2042057/1‑150 (MQ=255)
                                                                tGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGca                                                                                      >  1:113617/1‑150 (MQ=255)
                                                                 ggTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGcac                                                                                     >  1:3995064/1‑150 (MQ=255)
                                                                  gTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGcaca                                                                                    <  2:1798542/150‑1 (MQ=255)
                                                                     ccGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGtt                                                                                 <  2:2507875/150‑1 (MQ=255)
                                                                     ccGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGtt                                                                                 >  1:3351596/1‑150 (MQ=255)
                                                                      cggcggCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTg                                                                                <  2:2234872/150‑1 (MQ=255)
                                                                      cggcggCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTg                                                                                <  2:4133367/150‑1 (MQ=255)
                                                                       ggcggcGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGc                                                                               <  1:1972978/150‑1 (MQ=255)
                                                                       ggcggcGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGc                                                                               <  2:2192514/150‑1 (MQ=255)
                                                                        gcggcgCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGcc                                                                              <  2:2202620/150‑1 (MQ=255)
                                                                         cggcgCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGccc                                                                             <  2:2617200/150‑1 (MQ=255)
                                                                           gcgcTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTg                                                                           >  1:134874/1‑150 (MQ=255)
                                                                             gcTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTgcg                                                                         >  2:2893092/1‑150 (MQ=255)
                                                                                      cGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTg                                                                >  2:446014/1‑150 (MQ=255)
                                                                                       gTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTgg                                                               >  1:2745664/1‑150 (MQ=255)
                                                                                        tGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGa                                                              >  1:896025/1‑150 (MQ=255)
                                                                                        tGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGa                                                              >  1:3345562/1‑150 (MQ=255)
                                                                                        tGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGa                                                              <  1:1438279/150‑1 (MQ=255)
                                                                                          ggAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAAAGCCGCTCCCGCCGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATa                                                            >  1:746341/1‑150 (MQ=255)
                                                                                                gtcgtcAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTg                                                      <  2:567472/150‑1 (MQ=255)
                                                                                                   gtcAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGcaa                                                   >  1:1046380/1‑150 (MQ=255)
                                                                                                   gtcAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGcaa                                                   <  2:2733646/150‑1 (MQ=255)
                                                                                                      aaGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGcaacaa                                                >  2:400873/1‑150 (MQ=255)
                                                                                                       aGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACaaa                                               <  2:860768/150‑1 (MQ=255)
                                                                                                            cTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGtt                                          <  1:4111925/150‑1 (MQ=255)
                                                                                                             tGGTCAGCCACGCCGCCGCACCCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCTCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTc                                         <  2:3051290/150‑1 (MQ=255)
                                                                                                             tGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTc                                         <  2:2988420/150‑1 (MQ=255)
                                                                                                             tGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTc                                         <  2:2816310/150‑1 (MQ=255)
                                                                                                              ggTCAGCCACGCCGCCGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCg                                        <  2:4318802/150‑1 (MQ=255)
                                                                                                              ggTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCg                                        <  2:2377201/150‑1 (MQ=255)
                                                                                                               gTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCTTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGt                                       <  2:1958414/150‑1 (MQ=255)
                                                                                                                tCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGtt                                      >  2:1191865/1‑150 (MQ=255)
                                                                                                                   gCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAAt                                   <  2:3927821/150‑1 (MQ=255)
                                                                                                                    ccACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATc                                  >  2:3531123/1‑150 (MQ=255)
                                                                                                                    ccACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATc                                  >  2:777730/1‑150 (MQ=255)
                                                                                                                    ccACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATc                                  <  1:455610/150‑1 (MQ=255)
                                                                                                                    ccACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATc                                  <  1:702113/150‑1 (MQ=255)
                                                                                                                    ccACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATc                                  >  1:965602/1‑150 (MQ=255)
                                                                                                                    ccACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGGGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATc                                  >  1:3754348/1‑150 (MQ=255)
                                                                                                                        gCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGcac                              >  1:3669540/1‑150 (MQ=255)
                                                                                                                              gCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAg                        >  2:161369/1‑150 (MQ=255)
                                                                                                                              gCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAg                        >  1:2717509/1‑150 (MQ=255)
                                                                                                                              gCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAg                        >  1:2945258/1‑150 (MQ=255)
                                                                                                                               cAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGAGCTGCAATCAGCGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGt                       <  2:500854/150‑1 (MQ=255)
                                                                                                                                aaCCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTc                      <  1:3883336/150‑1 (MQ=255)
                                                                                                                                aaCCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTc                      <  1:2492260/150‑1 (MQ=255)
                                                                                                                                    tGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATAGCACAGATAGTCGAtt                  >  1:2016227/1‑150 (MQ=255)
                                                                                                                                      tgtgaatgggcCGCAACGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTAt                >  1:2349818/8‑150 (MQ=255)
                                                                                                                                      agggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTAt                <  2:686101/150‑1 (MQ=255)
                                                                                                                                        ggcgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCg              <  1:2805362/150‑1 (MQ=255)
                                                                                                                                          cgagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGa            >  2:1161636/1‑150 (MQ=255)
                                                                                                                                          aatgggcCGCAACGCGGCCCTTTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGCTAATCGCACAGATAGTCGATTATCGGa            >  1:3125071/4‑150 (MQ=255)
                                                                                                                                          aatgggcCGCAACGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGa            >  2:3693270/4‑150 (MQ=255)
                                                                                                                                            tgggcCGCAACGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGAtt          >  1:1938347/2‑150 (MQ=255)
                                                                                                                                            agggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGAtt          <  1:150068/150‑1 (MQ=255)
                                                                                                                                             gggcCGCAACGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGAGTg         >  1:4402516/1‑150 (MQ=255)
                                                                                                                                              ggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATTAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGt        >  1:406496/1‑150 (MQ=255)
                                                                                                                                                ccGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTaa      <  1:2399451/150‑1 (MQ=255)
                                                                                                                                                ccGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTaa      <  1:2442575/150‑1 (MQ=255)
                                                                                                                                                    gTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGta  <  1:2580630/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCTGGAGTGGTGCAACGCTGCTGGGCCTGGGCTGGAGCTTCGAGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCGACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTA  >  NC_002947/1568771‑1569068

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: