Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,457,725 G→T 100% 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,7250GT96.7% 682.4 / 11.4 222intergenic (+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 (5/3);  new base T (87/127);  total (92/130)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 222 total aligned reads displayed.

TTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTC  >  CP009273/1457580‑1457869
                                                                                                                                                 |                                                                                                                                                
ttAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTcc                                                                                                                                                >  1:279765/1‑148 (MQ=255)
 tAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCa                                                                                                                                               <  1:21727/148‑1 (MQ=255)
   aGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAgg                                                                                                                                             <  2:2900463/148‑1 (MQ=255)
    gTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAggg                                                                                                                                            >  1:3504788/1‑148 (MQ=255)
     tCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGc                                                                                                                                           <  2:2651836/148‑1 (MQ=255)
        aTGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAg                                                                                                                                        >  1:691532/1‑148 (MQ=255)
                 aTCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATa                                                                                                                               <  2:2620110/148‑1 (MQ=255)
                 aTCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATa                                                                                                                               <  2:1716068/148‑1 (MQ=255)
                 aTCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATa                                                                                                                               <  2:2697854/148‑1 (MQ=255)
                 aTCTCTACCGAGGTGATGATCGTTATCTGTGGCAGTATCCCGCGGTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTGTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATa                                                                                                                               <  2:2697131/148‑1 (MQ=255)
                  tctcTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATAc                                                                                                                              >  1:2081629/1‑148 (MQ=255)
                    tctACCGATGTGATGATCGTTAACTGTGTCAGTATCCCGCGTTCAGTAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCACTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGa                                                                                                                            <  2:130079/148‑1 (MQ=255)
                      tACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGAcc                                                                                                                          <  2:3465130/148‑1 (MQ=255)
                      tACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGAcc                                                                                                                          <  2:1727534/148‑1 (MQ=255)
                      tACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATGGTGGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGAcc                                                                                                                          <  1:943446/148‑1 (MQ=255)
                         cGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCtgt                                                                                                                       >  1:3006628/1‑148 (MQ=255)
                          gATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCtgtg                                                                                                                      <  2:3205516/148‑1 (MQ=255)
                          gATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCtgtg                                                                                                                      >  1:2899793/1‑148 (MQ=255)
                             gtgATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTa                                                                                                                   <  2:1987353/148‑1 (MQ=255)
                             gtgATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTa                                                                                                                   <  2:1146515/148‑1 (MQ=255)
                             gtgATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTa                                                                                                                   <  1:671036/148‑1 (MQ=255)
                               gatgatCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATg                                                                                                                 <  2:2157551/148‑1 (MQ=255)
                                atgatCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGa                                                                                                                <  2:798237/148‑1 (MQ=255)
                                 tgatCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGAt                                                                                                               <  1:3411741/148‑1 (MQ=255)
                                   atCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTc                                                                                                             <  1:352203/148‑1 (MQ=255)
                                   atCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTc                                                                                                             <  1:2697300/148‑1 (MQ=255)
                                       ttAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATaa                                                                                                         >  1:1690589/1‑148 (MQ=255)
                                       ttAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATaa                                                                                                         >  1:1690747/1‑148 (MQ=255)
                                       ttAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATaa                                                                                                         >  1:935004/1‑148 (MQ=255)
                                           cTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATCCGACCTGTGCTATGATTCATAAATca                                                                                                     >  2:336062/1‑148 (MQ=255)
                                             gtgGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCGAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATcaca                                                                                                   <  1:1242059/148‑1 (MQ=255)
                                             gtgGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATcaca                                                                                                   >  2:1660043/1‑148 (MQ=255)
                                              tgGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCacaa                                                                                                  >  2:577581/1‑148 (MQ=255)
                                                   gTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAATGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATa                                                                                             <  1:2245542/148‑1 (MQ=255)
                                                   gTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATa                                                                                             <  2:2799511/148‑1 (MQ=255)
                                                   gTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATa                                                                                             >  1:751698/1‑148 (MQ=255)
                                                      tCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaaca                                                                                          <  2:2504561/148‑1 (MQ=255)
                                                      tCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaaca                                                                                          <  1:3157690/148‑1 (MQ=255)
                                                       cccGCGGTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaacaa                                                                                         <  2:633140/148‑1 (MQ=255)
                                                          gcgTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAg                                                                                      <  2:2274008/148‑1 (MQ=255)
                                                          gcgTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAg                                                                                      <  2:2059109/148‑1 (MQ=255)
                                                          gcgTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAg                                                                                      <  2:2753888/148‑1 (MQ=255)
                                                           cgTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGa                                                                                     <  2:668188/148‑1 (MQ=255)
                                                           cgTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGa                                                                                     <  2:1275589/148‑1 (MQ=255)
                                                              tCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTg                                                                                  <  1:1452738/148‑1 (MQ=255)
                                                               cAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGa                                                                                 <  2:2750044/148‑1 (MQ=255)
                                                                 gAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGNTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAg                                                                               >  1:1780337/1‑148 (MQ=255)
                                                                   aggcaaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCg                                                                             <  2:503343/148‑1 (MQ=255)
                                                                         ggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAg                                                                       <  1:2455609/148‑1 (MQ=255)
                                                                          gcGTGTCGGGTGTTTCATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGt                                                                      >  2:930554/1‑148 (MQ=255)
                                                                            gtgtCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTaa                                                                    <  2:3074798/148‑1 (MQ=255)
                                                                            gtgtCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTaa                                                                    <  2:1872857/148‑1 (MQ=255)
                                                                             tgtCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTaaa                                                                   >  1:2275962/1‑148 (MQ=255)
                                                                             tgtCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTaaa                                                                   >  1:2761287/1‑148 (MQ=255)
                                                                                 gggTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTg                                                               <  2:238264/148‑1 (MQ=255)
                                                                                 gggTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTg                                                               <  2:2864258/148‑1 (MQ=255)
                                                                                  ggTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGa                                                              <  1:2897096/148‑1 (MQ=255)
                                                                                  ggTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGa                                                              <  1:1936599/148‑1 (MQ=255)
                                                                                   gtgtTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGAt                                                             >  1:2055625/1‑148 (MQ=255)
                                                                                     gtTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATAc                                                           >  1:285390/1‑148 (MQ=255)
                                                                                      tttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACt                                                          <  2:886315/148‑1 (MQ=255)
                                                                                      tttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACt                                                          <  2:3488560/148‑1 (MQ=255)
                                                                                      tttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACt                                                          <  1:2717466/148‑1 (MQ=255)
                                                                                      tttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGGGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACt                                                          <  1:1561788/148‑1 (MQ=255)
                                                                                       ttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACtt                                                         >  2:954603/1‑148 (MQ=255)
                                                                                        tGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACttt                                                        >  1:3095547/1‑148 (MQ=255)
                                                                                          aTCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTCAACTTGATACTTTTa                                                      >  2:1814129/1‑148 (MQ=255)
                                                                                            cTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATc                                                    <  2:1380932/148‑1 (MQ=255)
                                                                                              gcgcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCa                                                  <  1:2437562/148‑1 (MQ=255)
                                                                                                gcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAAc                                                <  2:1528488/148‑1 (MQ=255)
                                                                                                gcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAAc                                                <  1:3024290/148‑1 (MQ=255)
                                                                                                gcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAAc                                                <  2:1434503/148‑1 (MQ=255)
                                                                                                 cAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACa                                               <  2:2456927/148‑1 (MQ=255)
                                                                                                 cAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACa                                               >  2:140124/1‑148 (MQ=255)
                                                                                                     atTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCt                                           <  2:2899793/148‑1 (MQ=255)
                                                                                                      ttgttgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTg                                          >  2:541830/1‑148 (MQ=255)
                                                                                                       tgttgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGt                                         >  2:1843221/1‑148 (MQ=255)
                                                                                                        gttgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTa                                        >  1:2644161/1‑148 (MQ=255)
                                                                                                          tgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTaaa                                      <  1:1083728/148‑1 (MQ=255)
                                                                                                            gTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACg                                    <  2:2609491/148‑1 (MQ=255)
                                                                                                            gTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACg                                    <  2:851109/148‑1 (MQ=255)
                                                                                                             tGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGc                                   <  1:1392536/148‑1 (MQ=255)
                                                                                                                cTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGt                                >  2:154940/1‑148 (MQ=255)
                                                                                                                cTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGt                                >  2:3115301/1‑148 (MQ=255)
                                                                                                                 tGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGtt                               <  1:2314175/148‑1 (MQ=255)
                                                                                                                    cGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCg                            <  2:105523/148‑1 (MQ=255)
                                                                                                                        tCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCa                        >  2:3097712/1‑148 (MQ=255)
                                                                                                                          gTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGt                      <  2:2865014/148‑1 (MQ=255)
                                                                                                                           tCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTg                     >  2:1386014/1‑148 (MQ=255)
                                                                                                                            cTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTgg                    >  1:2780890/1‑148 (MQ=255)
                                                                                                                              ggCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGca                  >  2:1893407/1‑148 (MQ=255)
                                                                                                                               gCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGcac                 >  2:311995/1‑148 (MQ=255)
                                                                                                                               gCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGcac                 >  2:1274701/1‑148 (MQ=255)
                                                                                                                                 ccTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACAt               <  2:2647504/148‑1 (MQ=255)
                                                                                                                                      tggGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTg          <  1:101374/148‑1 (MQ=255)
                                                                                                                                        gggTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAt        <  2:1928137/148‑1 (MQ=255)
                                                                                                                                           tAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAtctc     <  2:2081629/148‑1 (MQ=255)
                                                                                                                                             aaGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATctcct   <  1:3362029/148‑1 (MQ=255)
                                                                                                                                              aGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATctcctc  >  1:1495611/1‑148 (MQ=255)
                                                                                                                                              aGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATctcctc  <  2:2205485/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: