Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,451,422 A→T 48.8% 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,4220AT48.8% 12.0 / inf 291T246S (ACG→TCG) paaEring 1,2‑phenylacetyl‑CoA epoxidase, NAD(P)H oxidoreductase component
Reads supporting (aligned to +/- strand):  ref base A (91/58);  new base T (79/63);  total (170/121)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

Only 100 of 293 total aligned reads displayed.

AATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGATCT  >  CP009273/1451273‑1451570
                                                                                                                                                     |                                                                                                                                                    
    aCAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGcc                                                                                                                                                  >  2:2395872/1‑150 (MQ=255)
     cAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATTCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCt                                                                                                                                                 <  2:1640744/150‑1 (MQ=255)
     cAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCt                                                                                                                                                 >  1:1726788/1‑150 (MQ=255)
      aGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTg                                                                                                                                                >  1:214616/1‑150 (MQ=255)
       gTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTgg                                                                                                                                               <  2:469605/150‑1 (MQ=255)
        tCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGc                                                                                                                                              >  1:3178697/1‑150 (MQ=255)
          aCTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCAc                                                                                                                                            >  1:4275577/1‑150 (MQ=255)
             tGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCAcgcg                                                                                                                                         >  2:4025508/1‑150 (MQ=255)
                ggCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGt                                                                                                                                      >  2:182751/1‑150 (MQ=255)
                ggCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGt                                                                                                                                      >  1:2809743/1‑150 (MQ=255)
                ggCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGt                                                                                                                                      >  1:417218/1‑150 (MQ=255)
                  ccTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCa                                                                                                                                    <  2:2029540/150‑1 (MQ=255)
                    tCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCaaa                                                                                                                                  >  1:554395/1‑150 (MQ=255)
                     cGCTCATTACTTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAAc                                                                                                                                 <  1:942514/150‑1 (MQ=255)
                       cTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGTTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTTGAGCGGTTTAATTCGCCAGGCACGAGAGTCAAACGt                                                                                                                               >  1:2410553/1‑150 (MQ=255)
                        tCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTa                                                                                                                              <  2:56013/150‑1 (MQ=255)
                          aTTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGc                                                                                                                            >  1:312002/1‑150 (MQ=255)
                           ttAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCg                                                                                                                           <  1:1970984/150‑1 (MQ=255)
                              aTTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTa                                                                                                                        <  2:847157/150‑1 (MQ=255)
                              aTTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTa                                                                                                                        >  1:1399337/1‑150 (MQ=255)
                               ttttTGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTaa                                                                                                                       >  2:4143575/1‑150 (MQ=255)
                                    gTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGc                                                                                                                  <  1:1955472/150‑1 (MQ=255)
                                     tCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCa                                                                                                                 <  1:2067837/150‑1 (MQ=255)
                                      cTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCaa                                                                                                                >  2:4371282/1‑150 (MQ=255)
                                        ttATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAg                                                                                                              >  1:1078179/1‑150 (MQ=255)
                                        ttATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAg                                                                                                              >  1:1173718/1‑150 (MQ=255)
                                          atgatgAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTg                                                                                                            >  2:742293/1‑150 (MQ=255)
                                           tgatgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGa                                                                                                           <  1:4359056/150‑1 (MQ=255)
                                            gatgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGAc                                                                                                          >  1:4087616/1‑150 (MQ=255)
                                              tgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACgg                                                                                                        <  1:2971345/150‑1 (MQ=255)
                                                 ggCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACa                                                                                                     <  2:3791280/150‑1 (MQ=255)
                                                  gCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaa                                                                                                    >  1:2211458/1‑150 (MQ=255)
                                                      ttatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAg                                                                                                <  2:3288078/150‑1 (MQ=255)
                                                        atttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTg                                                                                              >  1:472273/1‑150 (MQ=255)
                                                           tgtgGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACt                                                                                           >  1:2659104/1‑150 (MQ=255)
                                                             tgGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGt                                                                                         <  2:2157529/150‑1 (MQ=255)
                                                             tgGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGt                                                                                         <  1:1397489/150‑1 (MQ=255)
                                                              ggTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTa                                                                                        >  1:664707/1‑150 (MQ=255)
                                                                  cggcggCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTc                                                                                    <  1:2728764/150‑1 (MQ=255)
                                                                   ggcggcGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCa                                                                                   <  2:2289517/150‑1 (MQ=255)
                                                                     cggcgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAGTGACGGACAAAAAGTGACTGTACGTCAGGa                                                                                <  1:3979949/150‑1 (MQ=255)
                                                                       gcgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGAt                                                                               <  1:3732637/150‑1 (MQ=255)
                                                                         gatgatGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATgg                                                                             >  2:3849473/1‑150 (MQ=255)
                                                                           tgatgGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGc                                                                           >  2:3718083/1‑150 (MQ=255)
                                                                           tgatgGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGc                                                                           >  2:3814050/1‑150 (MQ=255)
                                                                           tgatgGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGc                                                                           >  1:999824/1‑150 (MQ=255)
                                                                              tggatgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCggg                                                                        <  1:1106475/150‑1 (MQ=255)
                                                                                 atgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATc                                                                     >  2:3687034/1‑150 (MQ=255)
                                                                                  tgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCg                                                                    <  1:4227267/150‑1 (MQ=255)
                                                                                     cgcgGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGa                                                                 >  1:1572168/1‑150 (MQ=255)
                                                                                     cgcgGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGa                                                                 <  2:903870/150‑1 (MQ=255)
                                                                                      gcgGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGaa                                                                <  2:3357339/150‑1 (MQ=255)
                                                                                        ggAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAAt                                                              >  1:3645277/1‑150 (MQ=255)
                                                                                          aaaCCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCg                                                            >  2:4351863/1‑150 (MQ=255)
                                                                                          aaaCCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGACATCg                                                            >  2:3490690/1‑150 (MQ=255)
                                                                                           aaCCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGt                                                           >  1:1078809/1‑150 (MQ=255)
                                                                                             ccGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGc                                                         <  2:2217598/150‑1 (MQ=255)
                                                                                               gCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCtt                                                       <  2:4312485/150‑1 (MQ=255)
                                                                                                ccTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTa                                                      >  2:1236062/1‑150 (MQ=255)
                                                                                                   tAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATg                                                   <  2:2924842/150‑1 (MQ=255)
                                                                                                    aaaaGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGc                                                  >  2:3960470/1‑150 (MQ=255)
                                                                                                       aGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcga                                               >  1:3376719/1‑150 (MQ=255)
                                                                                                       aGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcga                                               >  1:3045418/1‑150 (MQ=255)
                                                                                                         cACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcgacg                                             >  1:4037584/1‑150 (MQ=255)
                                                                                                         cACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcgacg                                             >  2:3398378/1‑150 (MQ=255)
                                                                                                          aCTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcgacga                                            >  1:1554166/1‑150 (MQ=255)
                                                                                                            tGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATg                                          >  2:636049/1‑150 (MQ=255)
                                                                                                             gggAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGa                                         >  1:2597860/1‑150 (MQ=255)
                                                                                                              ggAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGaa                                        <  2:1870373/150‑1 (MQ=255)
                                                                                                               gAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGaaa                                       <  2:3098379/150‑1 (MQ=255)
                                                                                                                aaTGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGaaat                                      >  1:3217491/1‑149 (MQ=255)
                                                                                                                   gCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATTCCGACGATGAAAGCAt                                   <  1:1346910/150‑1 (MQ=255)
                                                                                                                   gCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAt                                   <  2:312002/150‑1 (MQ=255)
                                                                                                                      aGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATtct                                >  2:4223326/1‑150 (MQ=255)
                                                                                                                        aTAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCg                              <  1:3612965/150‑1 (MQ=255)
                                                                                                                          aaaaCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGAt                            >  2:3221308/1‑150 (MQ=255)
                                                                                                                          aaaaCCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAGTCTCGAt                            >  2:836412/1‑150 (MQ=255)
                                                                                                                           aaaCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATg                           >  2:536512/1‑150 (MQ=255)
                                                                                                                           aaaCCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATg                           <  1:3639832/150‑1 (MQ=255)
                                                                                                                           aaaCCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATg                           >  1:114713/1‑150 (MQ=255)
                                                                                                                                aTTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGTATGGGCGGGATCGGGAAAACGTGCTTAATGCCGACGATGAAACCATACACGATGCGGCa                      >  1:1117585/1‑150 (MQ=255)
                                                                                                                                  tCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCAtt                    >  2:690482/1‑150 (MQ=255)
                                                                                                                                  tCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTATCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGACGACGATGAAAGCATTCTCGATGCGGCAtt                    >  1:1997994/1‑150 (MQ=255)
                                                                                                                                   cATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTg                   <  1:4423806/150‑1 (MQ=255)
                                                                                                                                    aTCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgc                  >  2:1452245/1‑150 (MQ=255)
                                                                                                                                    aTCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgc                  >  2:1311552/1‑150 (MQ=255)
                                                                                                                                    aTCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgc                  <  1:2295038/150‑1 (MQ=255)
                                                                                                                                    aTCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgc                  <  1:4658129/150‑1 (MQ=255)
                                                                                                                                     tCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgcg                 <  1:1822667/150‑1 (MQ=255)
                                                                                                                                     tCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgcg                 <  2:4529402/150‑1 (MQ=255)
                                                                                                                                      cTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgcgc                >  2:2323663/1‑150 (MQ=255)
                                                                                                                                       tGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGcc               >  2:1712642/1‑150 (MQ=255)
                                                                                                                                        ggAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCa              <  1:2716424/150‑1 (MQ=255)
                                                                                                                                        ggAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCa              >  2:113677/1‑150 (MQ=255)
                                                                                                                                               tttAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCg       >  2:2664374/1‑150 (MQ=255)
                                                                                                                                                ttAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCgg      >  2:1609959/1‑150 (MQ=255)
                                                                                                                                                 tAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGa     <  1:2350049/150‑1 (MQ=255)
                                                                                                                                                  aaTACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGAt    >  1:1237143/1‑150 (MQ=255)
                                                                                                                                                  aaTACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGAt    >  2:2904168/1‑150 (MQ=255)
                                                                                                                                                    tACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGATCt  <  1:692120/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGATCT  >  CP009273/1451273‑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: