Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,568,915 G→A 6.3% intergenic (+52/‑275) 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,9150GA 6.3% 237.3 / 2.6 79intergenic (+52/‑275)pcaK/pcaF‑I4‑hydroxybenzoate transporter/beta‑ketoadipyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base A (2/3);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

AGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCAC  >  NC_002947/1568772‑1569062
                                                                                                                                                   |                                                                                                                                                    
aGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGggcggc                                                                                                                                                 <  7:147431/149‑1 (MQ=255)
 gCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGgcggcg                                                                                                                                                <  3:302902/149‑1 (MQ=255)
  taGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCCCCTGCGGGAGCGGCGt                                                                                                                                               <  3:167006/148‑1 (MQ=255)
    ggTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtg                                                                                                                                             <  8:359349/149‑1 (MQ=255)
      tGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCTGCTCCCGCAGGGGCGGCGTtgtg                                                                                                                                           <  1:362585/149‑1 (MQ=255)
      tGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgtg                                                                                                                                           <  6:53548/149‑1 (MQ=255)
        cTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgtg                                                                                                                                           <  4:303263/147‑1 (MQ=255)
        cTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTtgtg                                                                                                                                           >  3:303263/1‑147 (MQ=255)
             ggCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAAt                                                                                                                                    >  5:311044/1‑149 (MQ=255)
                    cTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGgcg                                                                                                                             <  5:259124/149‑1 (MQ=255)
                     tGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGgcg                                                                                                                             >  3:345162/1‑148 (MQ=255)
                       gTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGca                                                                                                                           >  3:83197/1‑148 (MQ=255)
                        tACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGca                                                                                                                           >  6:66074/1‑147 (MQ=255)
                          ccGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGt                                                                                                                       >  6:146515/1‑149 (MQ=255)
                                gcgcTGGCGACCGTGGGTGTGGTCGTGAATGGGGGTGTCAGCCACGCCGCCTACAGAACCTGAGTGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGc                                                                                                                   <  1:359733/149‑1 (MQ=255)
                                  gcTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTgc                                                                                                              >  3:21650/1‑149 (MQ=255)
                                   cTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACACAACCTGAGGGCGAGGGCCGCGTTGCGGCCCTATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTgc                                                                                                              <  1:354559/149‑1 (MQ=255)
                                    tGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTgcgc                                                                                                            >  8:283716/1‑149 (MQ=255)
                                      gCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTgcgc                                                                                                            <  7:264337/147‑1 (MQ=255)
                                      gCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGt                                                                                                          >  6:181302/1‑149 (MQ=255)
                                        gACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTgg                                                                                                        >  4:327603/1‑149 (MQ=255)
                                        gACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTgg                                                                                                        >  1:242832/1‑149 (MQ=255)
                                        gACCGGGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTgg                                                                                                        <  5:96435/149‑1 (MQ=255)
                                           cGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGcac                                                                                                                          <  6:208420/129‑1 (MQ=255)
                                           cGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGcac                                                                                                                          >  5:208420/1‑129 (MQ=255)
                                            gTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTg                                                                                                    <  2:242832/149‑1 (MQ=255)
                                              gggAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGa                                                                                                  >  6:421444/1‑149 (MQ=255)
                                               ggAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGAt                                                                                                 <  4:21650/149‑1 (MQ=255)
                                                gAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGa                                                                                                  >  6:234481/1‑147 (MQ=255)
                                                 aGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAg                                                                                               <  8:350403/149‑1 (MQ=255)
                                                    ggTCGTCAAGGGGCTCGTCAGACACGCAG‑‑‑ACGCAACCTGAGGGCCAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCACGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCg                                                                                            >  6:332855/1‑149 (MQ=255)
                                                          cAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGcaa                                                                                       >  6:341050/1‑148 (MQ=255)
                                                           aaGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGcaaca                                                                                     <  4:345162/149‑1 (MQ=255)
                                                                gCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATgg                                                                                <  5:421444/149‑1 (MQ=255)
                                                                  tGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGtt                                                                              >  5:91940/1‑149 (MQ=255)
                                                                    gTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTc                                                                             <  7:283716/148‑1 (MQ=255)
                                                                    gTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCg                                                                            <  4:404956/149‑1 (MQ=255)
                                                                    gTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCg                                                                            <  1:85434/149‑1 (MQ=255)
                                                                     tCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGt                                                                           >  6:196842/1‑149 (MQ=255)
                                                                          cACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGAGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATc                                                                      <  4:131320/149‑1 (MQ=255)
                                                                           aCGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCg                                                                     >  5:226280/1‑149 (MQ=255)
                                                                              ccG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGca                                                                   >  4:218580/1‑148 (MQ=255)
                                                                                      gCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATa                                                             >  7:89111/1‑149 (MQ=255)
                                                                                      gCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAAAGGGTCGTTAATCGCACAGa                                                               >  4:395715/1‑147 (MQ=255)
                                                                                         aCCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTc                                                          <  6:311044/149‑1 (MQ=255)
                                                                                         aCCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTc                                                          <  5:181302/149‑1 (MQ=255)
                                                                                              tgtgaatgggcCGCAACGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTa                                                     >  3:402229/8‑149 (MQ=255)
                                                                                                   gagggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGa                                                <  5:341050/149‑1 (MQ=255)
                                                                                                     gggcCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGAtt                                              >  8:7640/1‑149 (MQ=255)
                                                                                                        ccGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGgc                                                                                                                              <  8:76742/67‑1 (MQ=255)
                                                                                                        ccGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGgc                                                                                                                              >  7:76742/1‑67 (MQ=255)
                                                                                                        ccGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTa                                           <  3:327603/149‑1 (MQ=255)
                                                                                                         cgcgTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGt                                            >  1:328129/1‑147 (MQ=255)
                                                                                                          gcgTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAg                                         <  7:7640/149‑1 (MQ=255)
                                                                                                          gcgTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAg                                         <  5:196842/149‑1 (MQ=255)
                                                                                                                 ggCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGtaccta                                  <  4:204505/149‑1 (MQ=255)
                                                                                                                  gCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGAtt                                              <  8:344711/136‑1 (MQ=255)
                                                                                                                  gCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGAtt                                              >  7:344711/1‑136 (MQ=255)
                                                                                                                   cccATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGAGCACAGTTGGCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGtacctacc                                <  3:250808/149‑1 (MQ=255)
                                                                                                                    ccATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGtacctac                                 >  6:271130/1‑147 (MQ=255)
                                                                                                                     cATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAg                              >  3:145440/1‑149 (MQ=255)
                                                                                                                           gcagcacaaCGCCGCCCCTGCGGGAGCGGCGTTATGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTcc                         >  5:202197/8‑149 (MQ=38)
                                                                                                                           acaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGtctc                          <  6:226280/148‑1 (MQ=255)
                                                                                                                           acaGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGtctc                          <  7:304965/148‑1 (MQ=255)
                                                                                                                             aGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCtt                       <  5:393530/149‑1 (MQ=255)
                                                                                                                                 tAACGCCGCTCCCGCAGGGGCGGGGTTTTGCTGCAATCAGGGCGCTCATGTGACCTGCGCGGGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGTATTGTAAGGTACATACCAGTCTcc                         <  6:390437/143‑1 (MQ=255)
                                                                                                                                 tAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTcc                         >  5:390437/1‑143 (MQ=255)
                                                                                                                                 tAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATc                   >  3:34614/1‑149 (MQ=255)
                                                                                                                                  aaCGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCt                  >  7:211245/1‑149 (MQ=255)
                                                                                                                                  aaCGCCGCCCCTGCGGGAGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCg                                                                            >  8:13287/1‑91 (MQ=38)
                                                                                                                                  aaCGCCGCCCCTGCGGGAGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCg                                                                            <  7:13287/91‑1 (MQ=38)
                                                                                                                                    cgccgcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTgg                <  5:66074/149‑1 (MQ=255)
                                                                                                                                    cgccgcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTgg                <  5:234481/149‑1 (MQ=255)
                                                                                                                                     gccgcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTggg               >  6:386519/1‑149 (MQ=255)
                                                                                                                                     gccgcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGACAGTCGATTATCGGATTGTAAGGtacct                                   >  7:32147/1‑129 (MQ=255)
                                                                                                                                        gcTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTGTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGtt             >  1:84301/1‑149 (MQ=255)
                                                                                                                                           cccGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGAGAGTCGATTATCGGATTGTAAGCTACCTACAAGTCTCCTTAATCTGGGATTAcc         >  2:349565/1‑149 (MQ=255)
                                                                                                                                            ccGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCg        <  4:402229/149‑1 (MQ=255)
                                                                                                                                              gCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGaa      <  6:91940/149‑1 (MQ=255)
                                                                                                                                                aGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGc    <  4:83197/149‑1 (MQ=255)
                                                                                                                                                  gggCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCAc  >  1:420386/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCAGGTGCTGACGGCGTTGCTGGTACCGGCGGCGCTGGCGACCGTGGGAGTGGTCGTCAAGGGGCTGGTCAGCCACGCCG‑‑‑ACGCAACCTGAGGGCGAGGGCCGCGTTGCGGCCCATTCACAGCGTAACGCCGCTCCCGCAGGGGCGGCGTTGTGCTGCAATCAGGGCGCACAGTTGCCCTGCGCGTGGGGTGGATAGTCGTGCAACAAATGGTTCGTTAATCGCACAGATAGTCGATTATCGGATTGTAAGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCAC  >  NC_002947/1568772‑1569062

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