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,7250GT99.0% 955.3 / ‑0.2 293intergenic (+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 (2/1);  new base T (132/158);  total (134/159)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 293 total aligned reads displayed.

GATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTC  >  CP009273/1457578‑1457871
                                                                                                                                                   |                                                                                                                                                  
gATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTcc                                                                                                                                                  <  1:1170809/150‑1 (MQ=255)
gATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTcc                                                                                                                                                  >  2:4002436/1‑150 (MQ=255)
   tcaGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCACTATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAACGCTCCAgg                                                                                                                                               <  2:1142681/148‑1 (MQ=255)
   tAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAgg                                                                                                                                               <  1:1864396/150‑1 (MQ=255)
    aaGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAggg                                                                                                                                              >  1:3215707/1‑150 (MQ=255)
     aGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGc                                                                                                                                             >  1:4198810/1‑150 (MQ=255)
        cGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAg                                                                                                                                          >  2:3925500/1‑150 (MQ=255)
           tggtggGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAg                                                                                                                                       <  1:1917020/150‑1 (MQ=255)
            ggtggGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGt                                                                                                                                      <  1:79622/150‑1 (MQ=255)
              tggGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCg                                                                                                                                    >  1:2923803/1‑150 (MQ=255)
                  gATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATAc                                                                                                                                <  2:2638246/150‑1 (MQ=255)
                  gATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATAc                                                                                                                                <  1:582310/150‑1 (MQ=255)
                  gATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATAc                                                                                                                                <  2:4340341/150‑1 (MQ=255)
                  gATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATAc                                                                                                                                >  1:286240/1‑150 (MQ=255)
                       ctACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGTGTGTTTGATCTGCTCTATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAATTCGATACGACCt                                                                                                                           <  2:1019903/150‑1 (MQ=255)
                        tACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCtg                                                                                                                          >  2:5184578/1‑150 (MQ=255)
                        tACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCtg                                                                                                                          <  1:4802764/150‑1 (MQ=255)
                         aCCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCtgt                                                                                                                         <  1:24741/150‑1 (MQ=255)
                            gATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCt                                                                                                                      <  2:1448070/150‑1 (MQ=255)
                            gATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCt                                                                                                                      >  2:212905/1‑150 (MQ=255)
                            gATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCt                                                                                                                      <  2:4155137/150‑1 (MQ=255)
                               gtgATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATg                                                                                                                   >  2:4535129/1‑150 (MQ=255)
                               gtgATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATg                                                                                                                   >  2:4042039/1‑150 (MQ=255)
                                 gatgatCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGAt                                                                                                                 <  2:2407084/150‑1 (MQ=255)
                                  atgatCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGAtt                                                                                                                <  1:1687607/150‑1 (MQ=255)
                                         ttAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAAt                                                                                                         <  1:5018506/150‑1 (MQ=255)
                                         ttAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAAt                                                                                                         >  1:3700042/1‑150 (MQ=255)
                                          tAACTGAGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGTGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATc                                                                                                        >  1:3497295/1‑150 (MQ=255)
                                            aCTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATcac                                                                                                      <  1:5231567/150‑1 (MQ=255)
                                             cTGTGGCAGTATCCCGCGTTCAGATGGCATGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATcaca                                                                                                     <  2:1485981/150‑1 (MQ=255)
                                             cTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATcaca                                                                                                     <  2:2889167/150‑1 (MQ=255)
                                              tgtgGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCacaa                                                                                                    <  2:3834005/150‑1 (MQ=255)
                                              tgtgGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGAGTCATAAATCacaa                                                                                                    >  2:1408470/1‑150 (MQ=255)
                                               gtgGCAGTATCCCGCGTTCAGACGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCacaac                                                                                                   <  2:4371307/150‑1 (MQ=255)
                                               gtgGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCacaac                                                                                                   >  2:1280177/1‑150 (MQ=255)
                                                 ggCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCAcaacaa                                                                                                 >  2:2317176/1‑150 (MQ=255)
                                                  gCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGAAGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAAt                                                                                                <  2:2214137/150‑1 (MQ=255)
                                                   cAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATa                                                                                               >  1:3848680/1‑150 (MQ=255)
                                                    aGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaa                                                                                              >  2:3088944/1‑150 (MQ=255)
                                                     gTATCCCGCGTTCAGTAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaac                                                                                             <  1:646955/150‑1 (MQ=255)
                                                       aTCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaacaa                                                                                           <  1:1067557/150‑1 (MQ=255)
                                                        tCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATaacaac                                                                                          <  1:2372250/150‑1 (MQ=255)
                                                           cgcgTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGa                                                                                       >  2:3599797/1‑150 (MQ=255)
                                                              gTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCacaaca                                                                                                  >  1:140805/1‑136 (MQ=255)
                                                               ttCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGa                                                                                   >  1:97090/1‑150 (MQ=255)
                                                                 cAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAg                                                                                 <  1:4484370/150‑1 (MQ=255)
                                                                  aGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGc                                                                                >  2:4586098/1‑150 (MQ=255)
                                                                  aGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGc                                                                                >  2:163523/1‑150 (MQ=255)
                                                                  aGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGc                                                                                <  2:862540/150‑1 (MQ=255)
                                                                      ggcaaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAAt                                                                            <  2:4667154/150‑1 (MQ=255)
                                                                      ggcaaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAAt                                                                            <  2:3417905/150‑1 (MQ=255)
                                                                         aaggcGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAg                                                                         <  2:4514241/150‑1 (MQ=255)
                                                                               tgtCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACt                                                                   >  1:5089531/1‑150 (MQ=255)
                                                                                gtCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACtt                                                                  >  2:1878283/1‑150 (MQ=255)
                                                                                  cGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGa                                                                >  1:3539676/1‑150 (MQ=255)
                                                                                    ggTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGTCCCTGGTGTGGTAAATCTCCAGGGACAGAAGTCGATACTACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCTAATGAGTAAACTTGATa                                                              <  2:4968050/150‑1 (MQ=255)
                                                                                    ggTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATa                                                              <  1:2152372/150‑1 (MQ=255)
                                                                                        tttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACttt                                                          >  1:4093298/1‑150 (MQ=255)
                                                                                        tttGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACttt                                                          <  1:4645388/150‑1 (MQ=255)
                                                                                          tGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTa                                                        >  2:626243/1‑150 (MQ=255)
                                                                                            ctCTGCGCAATATTGTTGGGGCCGGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATc                                                      <  1:2253326/149‑1 (MQ=255)
                                                                                             tCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATcc                                                     <  1:2035151/150‑1 (MQ=255)
                                                                                              cTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCa                                                    <  1:744091/150‑1 (MQ=255)
                                                                                                 cgcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACa                                                 >  1:5248947/1‑150 (MQ=255)
                                                                                                 cgcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACa                                                 >  2:2645677/1‑150 (MQ=255)
                                                                                                  gcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACAt                                                <  2:1387040/150‑1 (MQ=255)
                                                                                                  gcAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACAt                                                <  1:248532/150‑1 (MQ=255)
                                                                                                   cAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGACTGAGTAAACTTGATACTTTTATCCAACATg                                               >  2:4449532/1‑150 (MQ=255)
                                                                                                   cAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATg                                               >  1:912673/1‑150 (MQ=255)
                                                                                                   cAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATg                                               <  1:1606623/150‑1 (MQ=255)
                                                                                                   cAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATg                                               <  1:4934966/150‑1 (MQ=255)
                                                                                                    aaTATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGc                                              >  1:4257045/1‑150 (MQ=255)
                                                                                                    aaTATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGc                                              >  1:2810774/1‑150 (MQ=255)
                                                                                                    aaTATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGc                                              <  2:3190444/150‑1 (MQ=255)
                                                                                                    aaTATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGc                                              >  2:3223295/1‑150 (MQ=255)
                                                                                                     atatTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTCAACTTGATACTTTTATCCAACATGCt                                             >  1:4865732/1‑150 (MQ=255)
                                                                                                     atatTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCt                                             >  2:5272494/1‑150 (MQ=255)
                                                                                                     atatTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCt                                             >  2:4752585/1‑150 (MQ=255)
                                                                                                       atTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGt                                           >  1:3249364/1‑150 (MQ=255)
                                                                                                        ttgttgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTa                                          <  2:3700042/150‑1 (MQ=255)
                                                                                                         tgttgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTaa                                         >  1:4887943/1‑150 (MQ=255)
                                                                                                          gttgGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTaaa                                        <  2:1761880/150‑1 (MQ=255)
                                                                                                              gTGCCTGACGTATCGTCTGGCCCTGGTTAGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCAGAGCTATTATTAATAAATCACAACAAAAACAACAGACTGAATCGAATGAGTAAACTTGATACTGTTATCCAACATGCTGTAAACGCt                                    >  2:796847/1‑150 (MQ=255)
                                                                                                               tGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTg                                   >  2:2851102/1‑150 (MQ=255)
                                                                                                                  cTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTc                                >  1:1385843/1‑150 (MQ=255)
                                                                                                                   tGACGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTcc                               >  1:1765837/1‑150 (MQ=255)
                                                                                                                      cGTATCGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGt                            <  2:1691971/150‑1 (MQ=255)
                                                                                                                           cGTCTGGCCCTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTg                       >  2:2561403/1‑150 (MQ=255)
                                                                                                                                   ccTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCt               >  1:386226/1‑150 (MQ=255)
                                                                                                                                   ccTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCt               <  1:3925500/150‑1 (MQ=255)
                                                                                                                                   ccTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCt               >  2:905704/1‑150 (MQ=255)
                                                                                                                                    cTGGTGGGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCtt              <  1:1280177/150‑1 (MQ=255)
                                                                                                                                     tggtggGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCttt             <  1:3599797/150‑1 (MQ=255)
                                                                                                                                     tggtggGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCttt             <  2:2327385/150‑1 (MQ=255)
                                                                                                                                        tggGGTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAt          >  2:3760222/1‑150 (MQ=255)
                                                                                                                                           ggTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAtctc       >  1:3126930/1‑150 (MQ=255)
                                                                                                                                           ggTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGAtctc       >  1:1491474/1‑150 (MQ=255)
                                                                                                                                            gTAAAGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATctcc      <  1:2270742/150‑1 (MQ=255)
                                                                                                                                               aaGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCctct   >  1:1683626/1‑150 (MQ=255)
                                                                                                                                               aaGCTCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCctct   >  1:5021313/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GATTAAGTCGATGGTGGGGATCTCTACCGATGTGATGATCGTTAACTGTGGCAGTATCCCGCGTTCAGAAGGCAAGGCGTGTCGGGTGTTTGATCTGCGCAATATTGTTGGTGCCTGACGTATCGTCTGGCCCTGGTGGGGTAAAGCGCCAGGGCCAGAAGTCGATACGACCTGTGCTATGATTCATAAATCACAACAATAACAACAGACTGAAGCGAATGAGTAAACTTGATACTTTTATCCAACATGCTGTAAACGCTGTTCCGGTCAGTGGCACATCTTTGATCTCCTCTC  >  CP009273/1457578‑1457871

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: