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,7250GT97.8% 693.0 / 4.3 221intergenic (+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 (3/2);  new base T (87/129);  total (90/131)
Fisher's exact test for biased strand distribution p-value = 3.99e-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 221 total aligned reads displayed.

TAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTC  >  CP009273/1457581‑1457869
                                                                                                                                                |                                                                                                                                                
tAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCa                                                                                                                                               >  2:3481504/1‑148 (MQ=255)
   gTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAggg                                                                                                                                            >  1:530915/1‑148 (MQ=255)
   gTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAggg                                                                                                                                            <  2:1712580/148‑1 (MQ=255)
   gTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAggg                                                                                                                                            >  1:2426693/1‑148 (MQ=255)
      gATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCa                                                                                                                                         >  1:225102/1‑148 (MQ=255)
       aTGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAg                                                                                                                                        >  1:3009715/1‑148 (MQ=255)
             gggATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCg                                                                                                                                  >  1:720633/1‑148 (MQ=255)
             gggATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCg                                                                                                                                  <  1:852608/148‑1 (MQ=255)
              ggATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGa                                                                                                                                 >  1:2832746/1‑148 (MQ=255)
               gATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGAt                                                                                                                                <  2:2808885/148‑1 (MQ=255)
               gATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGAt                                                                                                                                >  2:1661269/1‑148 (MQ=255)
                  ctctACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCCGGGCCAGAAGTCGATACg                                                                                                                             >  2:3039785/1‑148 (MQ=255)
                      aactaTGTGATGATCGTTAACTGTGGAAGTATCCCGCGTTCATTAGGCAAGGCGTGTTGTGTGTTTGATCTGCGCTATTTTGTTGGTGCGTTTCGTTTCTTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCt                                                                                                                         <  1:2888768/144‑1 (MQ=255)
                         gATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCATGGCGTGTCGGGTGTTTTATCTGCGCAATATTGGTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCtgtg                                                                                                                      <  1:1112029/148‑1 (MQ=255)
                         gATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCtgtg                                                                                                                      <  1:2578042/148‑1 (MQ=255)
                         gATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCtgtg                                                                                                                      >  1:100727/1‑148 (MQ=255)
                          aTGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCGGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGc                                                                                                                     >  2:1309270/1‑148 (MQ=255)
                            gtgATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTa                                                                                                                   <  1:2031407/148‑1 (MQ=255)
                            gtgATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTa                                                                                                                   >  2:2813820/1‑148 (MQ=255)
                            gtgATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTa                                                                                                                   <  1:2498417/148‑1 (MQ=255)
                            gtgATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTa                                                                                                                   <  2:1502452/148‑1 (MQ=255)
                            gtgATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTa                                                                                                                   <  2:3339217/148‑1 (MQ=255)
                              gatgatCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATg                                                                                                                 <  2:2911976/148‑1 (MQ=255)
                              gatgatCGTTAACTGGGGCAGTATCCCGCGTTCAGACGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATg                                                                                                                 <  1:604152/148‑1 (MQ=255)
                                 gatCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGAtt                                                                                                              <  2:1338439/148‑1 (MQ=255)
                                  atCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTTTGCTATGATTc                                                                                                             >  1:1546190/1‑148 (MQ=255)
                                        aaCTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAAt                                                                                                       <  1:2168967/148‑1 (MQ=255)
                                         aCTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCTATACGACCTGTGCTATGATTAATAAATc                                                                                                      >  1:2271126/1‑148 (MQ=255)
                                          cTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATca                                                                                                     <  2:793687/148‑1 (MQ=255)
                                          cTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATca                                                                                                     >  2:620870/1‑148 (MQ=255)
                                          cTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATca                                                                                                     >  1:2203276/1‑148 (MQ=255)
                                          cTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGACAGAAGTCGATACGACCTGTGATATGATTCATAAATca                                                                                                     >  1:2160507/1‑148 (MQ=255)
                                               gCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCacaaca                                                                                                <  2:3133821/148‑1 (MQ=255)
                                               gCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCacaaca                                                                                                <  1:2428108/148‑1 (MQ=255)
                                                 aGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAAt                                                                                              <  1:172329/148‑1 (MQ=255)
                                                   tATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAACTCACAACAATaa                                                                                            >  2:1822067/1‑148 (MQ=255)
                                                   tATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaa                                                                                            <  1:1564385/148‑1 (MQ=255)
                                                     tCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaaca                                                                                          <  1:2183309/148‑1 (MQ=255)
                                                        cgcgTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAacaaca                                                                                       <  2:2712563/148‑1 (MQ=255)
                                                         gcgTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAg                                                                                      <  2:213167/148‑1 (MQ=255)
                                                          cgTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGa                                                                                     <  2:128258/148‑1 (MQ=255)
                                                             tCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTg                                                                                  <  1:1305636/148‑1 (MQ=255)
                                                             tCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTg                                                                                  >  1:1843552/1‑148 (MQ=255)
                                                                 aaggcaaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGc                                                                              <  2:2526088/148‑1 (MQ=255)
                                                                   ggcaaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGa                                                                            <  2:2397272/148‑1 (MQ=255)
                                                                     caaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAAt                                                                          >  2:1967088/1‑148 (MQ=255)
                                                                     caaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAAt                                                                          <  2:164343/148‑1 (MQ=255)
                                                                     caaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCGGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAAt                                                                          >  1:1366250/1‑148 (MQ=255)
                                                                       aggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGa                                                                        >  2:360071/1‑148 (MQ=255)
                                                                       aggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGa                                                                        >  1:60557/1‑148 (MQ=255)
                                                                        ggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCTGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAg                                                                       <  2:1674170/148‑1 (MQ=255)
                                                                        ggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAg                                                                       <  2:1104455/148‑1 (MQ=255)
                                                                         gcGTGTCGGGTGTTTGATCTGCGCAATATTGGGGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGt                                                                      <  1:1933500/148‑1 (MQ=255)
                                                                          cGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTa                                                                     <  1:163611/148‑1 (MQ=255)
                                                                          cGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTa                                                                     <  1:1961640/148‑1 (MQ=255)
                                                                            tgtCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTaaa                                                                   <  1:1829588/148‑1 (MQ=255)
                                                                             gtCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAAc                                                                  >  2:2888636/1‑148 (MQ=255)
                                                                              tCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGa                                                                                     >  1:2812579/1‑128 (MQ=255)
                                                                               cGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACtt                                                                >  2:356365/1‑148 (MQ=255)
                                                                               cGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACtt                                                                <  2:2991341/148‑1 (MQ=255)
                                                                                gggTGTTGGAACTGCGCAATAGTGTTGGTGCCTGACTTATCGTCTGGCCCTGGAGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTg                                                               <  1:1446787/148‑1 (MQ=255)
                                                                                  gtgtTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGAt                                                             >  2:2154729/1‑148 (MQ=255)
                                                                                  gtgtTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGAt                                                             <  2:2198891/148‑1 (MQ=255)
                                                                                    gtTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATAc                                                           >  1:2477316/1‑148 (MQ=255)
                                                                                    gtTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATAc                                                           <  1:2716367/148‑1 (MQ=255)
                                                                                      ttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACtt                                                         <  2:1853367/148‑1 (MQ=255)
                                                                                           cTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATc                                                    <  2:2066900/148‑1 (MQ=255)
                                                                                             gcgcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCa                                                  <  2:2961998/148‑1 (MQ=255)
                                                                                              cgcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCaa                                                 <  2:797275/148‑1 (MQ=255)
                                                                                              cgcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCaa                                                 <  2:354783/148‑1 (MQ=255)
                                                                                                cAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACa                                               <  2:250395/148‑1 (MQ=255)
                                                                                                 aaTATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGGAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACAt                                              <  1:484104/148‑1 (MQ=255)
                                                                                                      tgttgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGt                                         >  2:725607/1‑148 (MQ=255)
                                                                                                       gttgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTa                                        <  2:1730934/148‑1 (MQ=255)
                                                                                                           gTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACg                                    >  1:1852145/1‑148 (MQ=255)
                                                                                                             gCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCt                                  <  1:2928368/148‑1 (MQ=255)
                                                                                                                  aCGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTcg                             >  2:894498/1‑147 (MQ=255)
                                                                                                                  aCGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTcc                             >  1:2103676/1‑148 (MQ=255)
                                                                                                                  aCGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTcc                             >  1:3210528/1‑148 (MQ=255)
                                                                                                                  aCGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTcc                             >  1:506589/1‑148 (MQ=255)
                                                                                                                   cGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGACTGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCg                            >  1:1885822/1‑148 (MQ=255)
                                                                                                                   cGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCg                            <  2:1136285/148‑1 (MQ=255)
                                                                                                                    gTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCgg                           >  2:3102671/1‑148 (MQ=255)
                                                                                                                     tATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGt                          >  1:3556118/1‑148 (MQ=255)
                                                                                                                     tATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGt                          <  2:486757/148‑1 (MQ=255)
                                                                                                                       tCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCa                        >  2:1556934/1‑148 (MQ=255)
                                                                                                                          tCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTg                     >  2:449400/1‑148 (MQ=255)
                                                                                                                           cTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTgg                    >  2:3048476/1‑148 (MQ=255)
                                                                                                                            tGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGc                   <  1:664095/148‑1 (MQ=255)
                                                                                                                             ggCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCTATACGACCTGTTCTATGATTCATAAATCACAACAATAACAACAGACTTAAGCGAATGAGTAAAATTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGca                  <  2:1026475/148‑1 (MQ=255)
                                                                                                                             ggCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGca                  >  1:2910789/1‑148 (MQ=255)
                                                                                                                              gCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGcac                 >  2:2541909/1‑148 (MQ=255)
                                                                                                                              gCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGcac                 >  2:1280047/1‑148 (MQ=255)
                                                                                                                              gCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGcac                 >  1:3351469/1‑148 (MQ=255)
                                                                                                                                ccTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACAt               <  2:225102/148‑1 (MQ=255)
                                                                                                                                  tggtggGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCt             <  1:3031609/148‑1 (MQ=255)
                                                                                                                                    gtggGGTAAAGCTCCATGGCCAGTAGTCGATACGACCTGTGGTATGATTCATATATCACAACAATAACACCAGACTGTCGCGAATGAGTAAACTTGTTACTTTTATCCAACATGCTGTAAACGCTGTTCGGGTCAGTGGCACATCttt           <  2:2356453/148‑1 (MQ=255)
                                                                                                                                       gggTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAt        >  1:966095/1‑148 (MQ=255)
                                                                                                                                          tAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAtctc     <  1:2704235/148‑1 (MQ=255)
                                                                                                                                             aGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATctcctc  >  1:2385152/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTC  >  CP009273/1457581‑1457869

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