Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,451,422 A→T 59.1% T246S (ACG→TCG)  paaE → ring 1,2‑phenylacetyl‑CoA epoxidase, NAD(P)H oxidoreductase component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,451,4220AT59.1% 192.9 / inf 269T246S (ACG→TCG) paaEring 1,2‑phenylacetyl‑CoA epoxidase, NAD(P)H oxidoreductase component
Reads supporting (aligned to +/- strand):  ref base A (54/56);  new base T (81/78);  total (135/134)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 269 total aligned reads displayed.

TTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGATCT  >  CP009273/1451275‑1451570
                                                                                                                                                   |                                                                                                                                                    
ttACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCg                                                                                                                                                    >  1:2972992/1‑150 (MQ=255)
 tACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGc                                                                                                                                                   <  2:2841588/150‑1 (MQ=255)
  aCAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGcc                                                                                                                                                  >  1:385079/1‑150 (MQ=255)
   cAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCt                                                                                                                                                 >  1:3283295/1‑150 (MQ=255)
   cAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCt                                                                                                                                                 >  1:615967/1‑150 (MQ=255)
   cAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCt                                                                                                                                                 >  1:2274156/1‑150 (MQ=255)
     gTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTgg                                                                                                                                               >  1:3292071/1‑150 (MQ=255)
       cACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCa                                                                                                                                             >  2:3941229/1‑150 (MQ=255)
       cACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAACCACTGGGAATGCCACATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCa                                                                                                                                             >  2:184698/1‑150 (MQ=255)
         cTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCAcg                                                                                                                                           >  2:2284480/1‑150 (MQ=255)
           tGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCAcgcg                                                                                                                                         >  1:1578077/1‑150 (MQ=255)
            ggggCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCAcgcgc                                                                                                                                        <  2:926519/150‑1 (MQ=255)
              ggCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGt                                                                                                                                      >  1:1405547/1‑150 (MQ=255)
                   cGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAAc                                                                                                                                 >  2:4145034/1‑150 (MQ=255)
                     cTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGt                                                                                                                               >  2:1902915/1‑150 (MQ=255)
                     cTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGt                                                                                                                               <  2:3316851/150‑1 (MQ=255)
                     cTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGt                                                                                                                               <  2:2597966/150‑1 (MQ=255)
                          tAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGt                                                                                                                          >  2:954005/1‑150 (MQ=255)
                           aaTTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGtt                                                                                                                         >  2:2908194/1‑150 (MQ=255)
                            aTTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTa                                                                                                                        <  2:3178585/150‑1 (MQ=255)
                             ttttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTaa                                                                                                                       <  1:1700392/150‑1 (MQ=255)
                             ttttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTaa                                                                                                                       <  2:545311/150‑1 (MQ=255)
                             ttttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTaa                                                                                                                       >  2:2686666/1‑150 (MQ=255)
                              tttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAAc                                                                                                                      >  1:4103553/1‑150 (MQ=255)
                               ttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACg                                                                                                                     >  2:1618536/1‑150 (MQ=255)
                               ttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACg                                                                                                                     <  1:1795841/150‑1 (MQ=255)
                                  gTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGc                                                                                                                  <  1:1717703/150‑1 (MQ=255)
                                   tCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCa                                                                                                                 >  1:2362077/1‑150 (MQ=255)
                                    cTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCaa                                                                                                                >  1:3502014/1‑150 (MQ=255)
                                       tATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGt                                                                                                             <  2:460147/150‑1 (MQ=255)
                                            tgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACgg                                                                                                        >  2:3861017/1‑150 (MQ=255)
                                              aGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGAc                                                                                                      <  1:1029471/150‑1 (MQ=255)
                                               ggCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACa                                                                                                     <  2:733737/150‑1 (MQ=255)
                                                  atttatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaaaa                                                                                                  <  2:1339089/150‑1 (MQ=255)
                                                  atttatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaaaa                                                                                                  <  2:643301/150‑1 (MQ=255)
                                                  atttatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaaaa                                                                                                  >  2:2974000/1‑150 (MQ=255)
                                                      atttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTg                                                                                              >  1:524402/1‑150 (MQ=255)
                                                        ttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGAc                                                                                            <  1:222552/150‑1 (MQ=255)
                                                            ggTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTa                                                                                        >  1:903437/1‑150 (MQ=255)
                                                            ggTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTa                                                                                        >  2:3495705/1‑150 (MQ=255)
                                                              tCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACg                                                                                      >  1:3110492/1‑150 (MQ=255)
                                                              tCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACg                                                                                      <  1:323594/150‑1 (MQ=255)
                                                                cggcggCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTc                                                                                    >  1:566268/1‑150 (MQ=255)
                                                                cggcggCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTc                                                                                    >  1:3792137/1‑150 (MQ=255)
                                                                    ggcgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGa                                                                                >  2:3060010/1‑150 (MQ=255)
                                                                       gatgatGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATgg                                                                             >  2:1050338/1‑150 (MQ=255)
                                                                       gatgatGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATgg                                                                             >  2:2249333/1‑150 (MQ=255)
                                                                         tgatgGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAATCCATTCATCTGGAGCGGTTTAATTCGCCTGGAACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGTTGGGc                                                                           <  1:1928208/150‑1 (MQ=255)
                                                                         tgatgGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGc                                                                           >  1:2036218/1‑150 (MQ=255)
                                                                            tggatgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCggg                                                                        <  2:1695543/150‑1 (MQ=255)
                                                                             ggatgACGCGGACACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGa                                                                       <  2:415718/150‑1 (MQ=255)
                                                                             ggatgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGa                                                                       <  2:2472630/150‑1 (MQ=255)
                                                                               atgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATc                                                                     <  1:1504668/150‑1 (MQ=255)
                                                                                tgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCg                                                                    <  2:579965/150‑1 (MQ=255)
                                                                                   cgcgGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGa                                                                 <  1:871186/150‑1 (MQ=255)
                                                                                    gcgGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGaa                                                                <  2:3571338/150‑1 (MQ=255)
                                                                                      ggAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGaact                                                              >  2:2382677/1‑148 (MQ=255)
                                                                                       gAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATc                                                             >  1:126506/1‑150 (MQ=255)
                                                                                        aaaCCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCg                                                            <  2:1470416/150‑1 (MQ=255)
                                                                                        aaaCCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCg                                                            >  2:3552039/1‑150 (MQ=255)
                                                                                           ccGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGc                                                         >  2:2429411/1‑150 (MQ=255)
                                                                                           ccGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGc                                                         >  1:2077143/1‑150 (MQ=255)
                                                                                            cGCCTTAAAAGCACTGGGAAGGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCt                                                        <  2:3357151/150‑1 (MQ=255)
                                                                                                 tAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATg                                                   >  1:1865984/1‑150 (MQ=255)
                                                                                                   aaaGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGcc                                                 <  1:1612780/150‑1 (MQ=255)
                                                                                                   aaaGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGcc                                                 <  1:712891/150‑1 (MQ=255)
                                                                                                    aaGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcg                                                >  2:1925601/1‑150 (MQ=255)
                                                                                                     aGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcga                                               <  1:923511/150‑1 (MQ=255)
                                                                                                           gggAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGa                                         <  1:2306907/150‑1 (MQ=255)
                                                                                                             gAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGaaa                                       >  1:1216051/1‑150 (MQ=255)
                                                                                                             gAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGaaa                                       >  1:2455644/1‑150 (MQ=255)
                                                                                                              aaTGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAg                                      <  2:1136434/150‑1 (MQ=255)
                                                                                                              aaTGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGTACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAg                                      <  1:2798026/150‑1 (MQ=255)
                                                                                                                 gCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAt                                   <  1:1450449/150‑1 (MQ=255)
                                                                                                                 gCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAt                                   >  2:1260621/1‑150 (MQ=255)
                                                                                                                   cAGATAAAACCTTTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGCCGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATtc                                 <  2:857413/150‑1 (MQ=255)
                                                                                                                    aGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATtct                                <  1:393631/150‑1 (MQ=255)
                                                                                                                    aGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATtct                                <  1:189361/150‑1 (MQ=255)
                                                                                                                        aaaaCCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGAt                            >  1:2727914/1‑150 (MQ=255)
                                                                                                                         aaaCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATg                           <  1:2852223/150‑1 (MQ=255)
                                                                                                                          aaCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGc                          <  1:852287/150‑1 (MQ=255)
                                                                                                                          aaCCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGc                          >  2:3888683/1‑150 (MQ=255)
                                                                                                                           aCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCg                         >  2:882614/1‑150 (MQ=255)
                                                                                                                            ccATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCgg                        >  2:2740346/1‑150 (MQ=255)
                                                                                                                             cATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGc                       >  2:1228870/1‑150 (MQ=255)
                                                                                                                               ttCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCAt                     >  2:4138703/1‑150 (MQ=255)
                                                                                                                               ttCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCAt                     >  1:2603901/1‑150 (MQ=255)
                                                                                                                               ttCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCAt                     >  1:240437/1‑150 (MQ=255)
                                                                                                                                 cATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTg                   >  2:1518617/1‑150 (MQ=255)
                                                                                                                                   tCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgcg                 <  2:3755635/150‑1 (MQ=255)
                                                                                                                                     tGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGcc               >  1:3082912/1‑150 (MQ=255)
                                                                                                                                      ggAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCa              <  1:255970/150‑1 (MQ=255)
                                                                                                                                      ggAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCa              >  2:3615894/1‑150 (MQ=255)
                                                                                                                                      ggAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCa              <  2:385079/150‑1 (MQ=255)
                                                                                                                                       gAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAg             >  1:3373019/1‑150 (MQ=255)
                                                                                                                                        aGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGACATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAgg            >  2:754835/1‑150 (MQ=255)
                                                                                                                                        aGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAgg            >  2:754789/1‑150 (MQ=255)
                                                                                                                                           ggTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAggggg         >  1:2348440/1‑150 (MQ=255)
                                                                                                                                              ttAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCgg      <  1:3637123/150‑1 (MQ=255)
                                                                                                                                                 aTACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGATc   <  2:2274156/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGATCT  >  CP009273/1451275‑1451570

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: