Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,457,725 G→T 95.0% intergenic (+30/‑71) paaK → / → paaX phenylacetyl‑CoA ligase/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,7250GT95.0% 737.0 / 24.9 249intergenic (+30/‑71)paaK/paaXphenylacetyl‑CoA ligase/transcriptional repressor of phenylacetic acid degradation paa operon, phenylacetyl‑CoA inducer
Reads supporting (aligned to +/- strand):  ref base G (6/7);  new base T (108/128);  total (114/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

TTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTC  >  CP009273/1457580‑1457869
                                                                                                                                                 |                                                                                                                                                
ttAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTcc                                                                                                                                                >  2:210582/1‑148 (MQ=255)
    gTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAggg                                                                                                                                            >  2:188412/1‑148 (MQ=255)
     tCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATGTTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGc                                                                                                                                           <  2:559484/148‑1 (MQ=255)
      cGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGcc                                                                                                                                          <  2:3234076/148‑1 (MQ=255)
      cGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGcc                                                                                                                                          >  2:3511686/1‑148 (MQ=255)
      cGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTAACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGcc                                                                                                                                          >  2:1948538/1‑148 (MQ=255)
       gATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCa                                                                                                                                         <  2:1606344/148‑1 (MQ=255)
        aTGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAg                                                                                                                                        >  1:3541079/1‑148 (MQ=255)
        aTGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAg                                                                                                                                        >  2:895853/1‑148 (MQ=255)
        aTGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAg                                                                                                                                        >  1:1419363/1‑148 (MQ=255)
         tggtggGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGa                                                                                                                                       >  1:2405641/1‑148 (MQ=255)
         tggtggGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGa                                                                                                                                       <  1:1704173/148‑1 (MQ=255)
         tggtggGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGa                                                                                                                                       >  2:1295829/1‑148 (MQ=255)
           gtggGGATCTCTACCGAGGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAg                                                                                                                                     <  2:2699070/148‑1 (MQ=255)
            tggGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATAGTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGt                                                                                                                                    >  2:209105/1‑148 (MQ=255)
             ggggATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTAACATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTc                                                                                                                                   >  1:1993962/1‑148 (MQ=255)
              gggATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCg                                                                                                                                  <  1:276377/148‑1 (MQ=255)
              gggATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCAGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCTTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCg                                                                                                                                  >  2:3281080/1‑148 (MQ=255)
               ggATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGa                                                                                                                                 >  2:3132645/1‑148 (MQ=255)
                gATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGAt                                                                                                                                >  2:1818068/1‑148 (MQ=255)
                  tctcTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAGGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATAc                                                                                                                              >  1:1619144/1‑148 (MQ=255)
                  tctcTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATAc                                                                                                                              >  1:3377823/1‑148 (MQ=255)
                      tACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGAcc                                                                                                                          <  2:720146/148‑1 (MQ=255)
                         cGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCtgt                                                                                                                       >  2:2548514/1‑148 (MQ=255)
                            tgtgATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCt                                                                                                                    >  1:68741/1‑148 (MQ=255)
                              tgatgaTCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTAt                                                                                                                  >  1:3419409/1‑148 (MQ=255)
                               gatgatCGTTAACTGTGGCAGTATCCCGCGTTCAGATGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCTTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATg                                                                                                                 <  1:2207157/148‑1 (MQ=255)
                               gatgatCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATg                                                                                                                 <  2:219269/148‑1 (MQ=255)
                                atgatCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGa                                                                                                                <  2:1708397/148‑1 (MQ=255)
                                 tgatCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGAt                                                                                                               >  1:3272850/1‑148 (MQ=255)
                                  gatCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGAtt                                                                                                              <  1:2915403/148‑1 (MQ=255)
                                   atCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTc                                                                                                             <  1:1537400/148‑1 (MQ=255)
                                     cGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCAt                                                                                                           >  2:83946/1‑148 (MQ=255)
                                     cGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCAt                                                                                                           >  1:2458825/1‑148 (MQ=255)
                                          aCTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATc                                                                                                      <  1:2035281/148‑1 (MQ=255)
                                          aCTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATc                                                                                                      <  2:607902/148‑1 (MQ=255)
                                            tgtgGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATcac                                                                                                    <  1:1740869/148‑1 (MQ=255)
                                               ggCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCacaac                                                                                                 >  1:1946888/1‑148 (MQ=255)
                                                gCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTAGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCacaaca                                                                                                <  1:851026/148‑1 (MQ=255)
                                                 cAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCAcaacaa                                                                                               >  2:2476783/1‑148 (MQ=255)
                                                    tATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaa                                                                                            >  2:2125399/1‑148 (MQ=255)
                                                    tATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaa                                                                                            <  2:503368/148‑1 (MQ=255)
                                                    tATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaa                                                                                            >  2:852386/1‑148 (MQ=255)
                                                     aTCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaac                                                                                           <  1:102457/148‑1 (MQ=255)
                                                     aTCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaac                                                                                           <  2:2444428/148‑1 (MQ=255)
                                                      tCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaaca                                                                                          <  1:2265885/148‑1 (MQ=255)
                                                          gcgTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAg                                                                                      >  2:3432604/1‑148 (MQ=255)
                                                             ttCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACt                                                                                   >  1:2838879/1‑148 (MQ=255)
                                                             ttCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACt                                                                                   >  1:1142692/1‑148 (MQ=255)
                                                               cAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGa                                                                                 <  1:2281930/148‑1 (MQ=255)
                                                               cAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGa                                                                                 <  2:3292755/148‑1 (MQ=255)
                                                               cAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGGGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGa                                                                                 <  1:39008/148‑1 (MQ=255)
                                                                aGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGaa                                                                                <  1:307157/148‑1 (MQ=255)
                                                                 gAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAg                                                                               >  2:333797/1‑148 (MQ=255)
                                                                  aaggcaaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGc                                                                              <  2:1791343/148‑1 (MQ=255)
                                                                  aaggcaaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGc                                                                              >  2:3007754/1‑148 (MQ=255)
                                                                  aaggcaaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGc                                                                              >  2:1111681/1‑148 (MQ=255)
                                                                    ggcaaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGa                                                                            <  2:946936/148‑1 (MQ=255)
                                                                    ggaaaggcGTGTCGGGTGTTTGATCTGCGCATTATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGa                                                                            <  1:2224965/145‑1 (MQ=255)
                                                                     gcaaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGaa                                                                           <  1:1089855/148‑1 (MQ=255)
                                                                     gcaaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGaa                                                                           <  1:2862963/148‑1 (MQ=255)
                                                                       aaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATg                                                                         <  2:594605/148‑1 (MQ=255)
                                                                        aggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGa                                                                        >  2:122038/1‑148 (MQ=255)
                                                                        aggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGa                                                                        <  2:2465526/148‑1 (MQ=255)
                                                                          gcGTGTCGGGTGTTTGATCTGCGCAATATTGGTGGGGCCGGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGt                                                                      <  1:2701032/148‑1 (MQ=255)
                                                                               tCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACt                                                                 <  1:454121/148‑1 (MQ=255)
                                                                                  ggTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCTGGGACAGAAGTATATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGa                                                              <  1:1343685/148‑1 (MQ=255)
                                                                                  ggTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGa                                                              >  2:1165458/1‑148 (MQ=255)
                                                                                    tgtTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATa                                                            >  1:1219271/1‑148 (MQ=255)
                                                                                    tgtTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATa                                                            <  1:684106/148‑1 (MQ=255)
                                                                                      tttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGGGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACACTAACAACAGACTGAAGCGAATGAGTAAACTTGATACt                                                          >  2:1350392/1‑148 (MQ=255)
                                                                                       ttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAATCTTGATACtt                                                         >  2:2201846/1‑148 (MQ=255)
                                                                                       ttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACtt                                                         <  1:3303435/148‑1 (MQ=255)
                                                                                       ttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACtt                                                         <  1:1097233/148‑1 (MQ=255)
                                                                                          aTCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTa                                                      >  2:1226553/1‑148 (MQ=255)
                                                                                           tCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTAt                                                     <  2:144754/148‑1 (MQ=255)
                                                                                           tCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTAt                                                     <  1:1818103/148‑1 (MQ=255)
                                                                                              gcgcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCa                                                  >  1:559289/1‑148 (MQ=255)
                                                                                              gcgcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCa                                                  >  1:1325057/1‑148 (MQ=255)
                                                                                                 cAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACa                                               >  2:1908099/1‑148 (MQ=255)
                                                                                                   atatTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATg                                             <  1:209105/148‑1 (MQ=255)
                                                                                                      ttgttgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTg                                          <  2:2809539/148‑1 (MQ=255)
                                                                                                        gttgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTa                                        <  1:127973/148‑1 (MQ=255)
                                                                                                          tgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTaaa                                      >  2:724651/1‑148 (MQ=255)
                                                                                                            gTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACg                                    <  2:2029583/148‑1 (MQ=255)
                                                                                                             tGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGc                                   <  2:1564496/148‑1 (MQ=255)
                                                                                                             tGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGc                                   <  2:2337693/148‑1 (MQ=255)
                                                                                                                  gACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGCTACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTc                              <  1:3331193/148‑1 (MQ=255)
                                                                                                                  gACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTc                              <  1:1485098/148‑1 (MQ=255)
                                                                                                                  gACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTc                              <  1:1884782/148‑1 (MQ=255)
                                                                                                                    cGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCg                            >  2:2371807/1‑148 (MQ=255)
                                                                                                                         cGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAg                       >  2:3315598/1‑148 (MQ=255)
                                                                                                                          gTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGt                      <  1:1133636/148‑1 (MQ=255)
                                                                                                                                   tggtggGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCt             <  1:3218891/148‑1 (MQ=255)
                                                                                                                                   tggtggGGTAAAGCGCCTGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCt             <  1:3007292/148‑1 (MQ=255)
                                                                                                                                      tggGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTg          <  2:2900742/148‑1 (MQ=255)
                                                                                                                                      tggGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTg          <  2:2626842/148‑1 (MQ=255)
                                                                                                                                        gggTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAt        <  2:3199955/148‑1 (MQ=255)
                                                                                                                                           tAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAtctc     <  2:2094529/148‑1 (MQ=255)
                                                                                                                                            aaaGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATctcc    <  1:290085/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTC  >  CP009273/1457580‑1457869

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: